% Loese das Roessler-System mit RK2

T  = 1000;           % Stop-Zeitpunkt
dt = 0.1;            % Zeitschritt
N  = floor(T/dt);    % Dann brauchen wir N Schritte

x = zeros(3,N);      % Fuer bessere Effizienz preallokiere Loesungsvektor
x(:,1) = [0; 0; 0];  % Anfangsdaten

for j = 1:N
  F1 = dt * f(x(:,j));
  F2 = dt * f(x(:,j)+F1);
  x(:,j+1)  = x(:,j) + (F1+F2)/2; 
end

gset parametric
gsplot (x(:,2000:N))'