function dX = rossler(X,a,b,c)
dX = [    -X(2) - X(3);          % dx/dt
     X(1) + a*X(2);        % dy/dt
     b + X(3)*(X(1) - c)]; % dz/dt
end