MWORKS Syslab Online

# Lorenz System - Sensitive Dependence
using MWorks, Plots, DifferentialEquations

function lorenz!(du, u, p, t)
    σ, ρ, β = p
    du[1] = σ * (u[2] - u[1])
    du[2] = u[1] * (ρ - u[3]) - u[2]
    du[3] = u[1] * u[2] - β * u[3]
end

u0 = [1.0; 0.0; 0.0]
p = [10.0, 28.0, 8/3]
tspan = (0.0, 50.0)

prob = ODEProblem(lorenz!, u0, tspan, p)
sol = solve(prob)

plot(sol, vars=(1,2,3), lw=0.5)
title!("Lorenz Attractor")