In [1]:
using DifferentialEquations
f = @ode_def_nohes LorenzExample begin
dx = σ*(y-x)
dy = x*(ρ-z) - y
dz = x*y - β*z
end σ ρ β
u0 = big.([0.1;0.0;0.0])
tspan = (big(0.0),big(100.0))
prob = ODEProblem(f,u0,tspan,(10.0,28.0,2.6666))
sol = solve(prob);
In [2]:
using Plots; gr(); plot(sol)