## I have written the Julia example for solving the Lorenz equation as in the page:

Ordinary Differential Equations · DifferentialEquations.jl :

function lorenz!(du,u,p,t)

du[1] = 10.0*(u[2]-u[1])

du[2] = u[1]*(28.0-u[3]) - u[2]

du[3] = u[1]*u[2] - (8/3)*u[3]

end

## 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.