DiffEq Solutions

Solution to the Lorenz Problem

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)