terrame / sysdyn Goto Github PK

import ("sysdyn")

Lorenz = SysDynModel{
    x  = 1.0, y = 1.0, z = 1.0,
    ro   = 28.0, gamma = 10.0, beta = 2.6666667,
    yplot = 0.0,
    delta = 0.01,
    deltaTime = 1,
    finalTime = 10000,
    changes = function(model, time)
        local nx = model.x + model.delta*model.gamma*(model.y - model.x)
        local ny = model.y + model.delta*(model.ro*model.x - model.y - model.x*model.z)
        local nz = model.z + model.delta*(model.x*model.y - model.beta*model.z)

        model.x = nx
        model.y = ny
        model.z = nz

        if (time % 100) then model.yplot = model.y end
    end,
    graphics = {
        timeseries = {{"yplot"}},
        phasespace = {{"x","z"}, {"x","y"}}
    }
}

Lorenz:execute()