(*Adapted by Rafael Luís: rafaelluis@netmadeira.com
2016-09-11: updated for Mathematica v.10 by Ross Moore  *)

(*Rossler attractor for (x,y,z)=(-y-z,x+ay,b+xz-cz)*)

f={X,Y,Z}/.NDSolve[{
X'[t]==-(Y[t]+Z[t]),
Y'[t]==X[t]+a Y[t],
Z'[t]==b+X[t]Z[t]-c Z[t],
X[0]==1,Y[0]==1,Z[0]==1}/.{a->.2,b->.2,c->8},
{X,Y,Z},{t,0,50}][[1]]

<<Graphics`Colors`

ParametricPlot3D[Evaluate[#[t]&/@f],{t,0,50},BoxRatios->{1,1,1},PlotRange->All,PlotPoints->2000,AxesLabel->TraditionalForm/@l]