Wolfram / ローレンツアトラクター

sol=NDSolve[{x'[t]==10(y[t]-x[t]),y'[t]==x[t](28-z[t])-y[t],z'[t]==x[t]y[t]-(8/3)z[t],x[0]==1,y[0]==0,z[0]==0},{x[t],y[t],z[t]},{t,0,50}];
p[t_]={x[t],y[t],z[t]}/.sol[[1]];
lorenz=ParametricPlot3D[p[t],{t,0,50},PlotRange->All]

参考: レクチャーズオンMathematica, pp.141-142