I used the following commands:
s = NDSolve[{X'[t] == 10 (X[t] - Y[t]),
Y'[t] == X[t] (28 - Z[t]) - Y[t], Z'[t] == X[t] Y[t] - (8/3) Z[t],
X[0] == Y[0] == Z[0] == 0.001}, {X, Y, Z}, {t, 100}]
ParametricPlot3D[Evaluate[{X[t], Y[t], Z[t]} /. s], {t, 0, 100}]
but I am not sure if I am right, the results are strange. Am I right? How can I plot $X,Y,Z?$
PlotRange
:ParametricPlot3D[Evaluate[{X[t], Y[t], Z[t]} /. s], {t, 0, 100}, PlotRange -> 40]
-- looks like it then goes off to infinity instead of getting orbiting another critical point. $\endgroup$Plot
, look in the documentation forNDSolve
, or look at these examples: mathematica.stackexchange.com/questions/134222/… $\endgroup$ParametricPlot3D
withPlot
? -- Also, if you want the Lorenz system, I think the first equation has X and Y switched. $\endgroup$