Lorenz方程式

NDSolve[{
x'[t]==-10 (x[t]-y[t]),
y'[t]==-x[t] z[t]+28 x[t]-y[t],
z'[t]==x[t] y[t]-8/3 z[t],
x[0]==z[0]==0,y[0]==1},{x,y,z},{t,0,200},MaxSteps->Infinity]

ParametricPlot3D[Evaluate[{x[t],y[t],z[t]}/.%],{t,0,200},PlotPoints->10000,ColorFunction->(ColorData["Rainbow"][#4]&), AxesLabel->{x,y,z}]