I have written the Julia example for solving the Lorenz equation as in the page:
Ordinary Differential Equations · DifferentialEquations.jl :
du = 10.0*(u-u)
du = u*(28.0-u) - u
du = u*u - (8/3)*u
u0 = [1.0;0.0;0.0]
tspan = (0.0,100.0)
prob = ODEProblem(lorenz!,u0,tspan)
sol = solve(prob)
The problem appears when plotting because I obtain two different results when I was expecting the same result. Below I have attached the piece of code together with the plots.
In the code, x[1,:],x[2,:] and x[3,:] are respectively the x,y,z components of the solution points.
Thank you very much in advance.