s=NDSolve[{x'[t]==-(y[t]+z[t]),y'[t]==x[t] +0.2*y[t],z'[t]==.2+x[t]*z[t]-5.7*z[t],
x[0]==.3,z[0]==4.2,y[0]==1.3},{x,y,z},{t,0,400},MaxSteps->12000];
S=ParametricPlot3D[Evaluate[{x[t],y[t],z[t]}/.s],{t,0,400}]
