clear;

a = 0.2;
b = 0.2;
c = 5.6;

function dx = rossler(t,x)
  dx(1) = - x(2) - x(3);
  dx(2) = x(1) + a * x(2);
  dx(3) = b + x(1) * x(3) - c * x(3);
endfunction

T = linspace(0,100,10000);
x0 = 1;
y0 = 0;
z0 = 0;
X0 = [x0; y0; z0];

X = ode(X0,0,T,rossler);
plot3d3(X(1,:),X(2,:),X(3,:));