% Der harmonische Oszillator - explizites Eulerverfahren
%
% Version 1: Als ein System von zwei skalaren Gleichungen
% geschrieben 

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

x = zeros(N,1);
y = zeros(N,1);
y(1) = 1;

for i=1:N-1
  x(i+1) = x(i) + dt*y(i);
  y(i+1) = y(i) - dt*x(i);
end

plot(x,y), grid, replot 
