%Sample script file to compute and plot a sample
% three dimensional orbit.

k=[0.2 0.2 5];
x(1)=0.1; y(1)=0.1; z(1)=0.3;

for n=1:500
    x(n+1)=x(n)-(y(n)+z(n))*0.1;
    y(n+1)=y(n)+(x(n)+k(1)*y(n))*0.1;
    z(n+1)=z(n)+(k(2)+z(n)*(x(n)-k(3)))*0.1;
end

plot3(x,y,z);
title('The Rossler Equation');

%Other plot options:
%plot3(x,y,z,'ro-');