Lorenz attractor

Implementation

NDSolve[ {x'[t] == -3 (x[t] - y[t]),
          y'[t] == -x[t] z[t] + 26.5 x[t] - y[t],
          z'[t] == x[t] y[t] - z[t],
          x[0] == z[0] == 0, y[0] == 1},
         {x, y, z}, {t, 0, 20}, MaxSteps->3000 ]

ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. %],
                          {t, 0, 20}, PlotPoints -> 1000]