Beispiel - SAM - ETH Zürich

Bemerkung 2.1.4 (Notation fuer Einschrittverfahren).
Numerische
Mathemtik
Oft spezifiziert man Einschrittverfahren durch Angabe des ersten Schritts
y 1 = Ausdruck iny 0 undf .
Dieser Gepflogenheit wird sich auch diese Vorlesung manchmal anschliessen.
△
Einschrittverfahren
function [T,Y] = esv(Psi,tspan,y0)
t = tspan(1); y = y0; Y = y; T = t;
while (t < tspan(2))
h = Aktuelle Zeitschrittweite
y = Psi(t,t+h,y); t = t+h;
Y = [Y,y]; T = [T,t];
end
Funktionshandle
Psi = @(t0,t1,y) ...
Beachte: Die aktuelle Zeitschrittweite
wird jeweils aus den Genauigkeitsanforderungen
und y k berechnet(→
Sect. 2.6).
R. Hiptmair
rev 35327,
25. April
2011
2.1
p. 120