Julia Scientific Computing

using DifferentialEquations
using Plots

function lorenz!(du, u, p, t)
    sigma, rho, beta = p
    du[1] = sigma * (u[2] - u[1])
    du[2] = u[1] * (rho - u[3]) - u[2]
    du[3] = u[1] * u[2] - beta * 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, Tsit5())

plot(sol, vars=(1,2,3),
     title="Lorenz Attractor",
     label="Trajectory")