function X = rossler_gen(time)
a = 0.2;
b = 0.2;
c = 5.7;


% Initial condition
X0 = [1; 1; 1];
% Solve ODE
[t, X] = ode45(@(t,X) rossler(X,a,b,c), time, X0);

for j=1:3
X(:,j) = X(:,j) - mean(X(:,j));
X(:,j) = X(:,j)/std(X(:,j));
end
end