% "Experimentelle" Bestimmung der Ordnung eines Verfahrens
%
% Hier am Beispiel RK2 fuer die Logistische Differentialgleichung

T = 1;            % Stop-Zeitpunkt
NN = 2.^(1:10);   % Vektor mit allen Schrittanzahlen

for i = 1:10
  
  N  = NN(i);     % Anzahl Schritte
  dt = T/N;       % Zeitschritt
  x  = 1/2;       % Anfangswert  

  for j = 1:N
    F1 = dt * x * (1-x);
    F2 = dt * (x+F1) * (1-(x+F1));
    x  = x + (F1+F2)/2;
  end
   
  error (i) = abs (x - e/(1+e));

end

loglog (NN, error, "*"), hold on
loglog (NN, NN.^(-2))  % Plotte eine Gerade der Steigung -2 zum Vergleich
grid, replot, hold off





