Beispiel - SAM - ETH Zürich

3 while ((t(end) < T) (h > hmin)) %
4 yh = Psihigh(h,y0); % ESV hoher Ordnung
5 yH = Psilow(h,y0); % ESV niedriger Ordnung
6 est = norm(yH-yh); %
Numerische
Mathemtik
7
8 i f (est < max(reltol*norm(y0),abstol)) %
9 y0 = yh; y = [y,y0]; h = 1.1*h; % Schritt akzeptiert
10 h = min(h,T-t(end)); t=[t,t+h]; %
11 else, h = h/2; end % Schritt verworfen
12 end
R. Hiptmair
Kommentare zu Code 2.4:
rev 35327,
25. April
2011
• Argumente vonodeintadapt:
–Psilow,Psihigh: Funktionshandles auf diskrete Evolutionen (für autonome ODE)
unterschiedlicher Ordnungen, Typ@(y,h), Zustandsvektor als erstes Argument, Schrittweite
als zweites,
–T: EndzeitpunktT > 0,
2.6
p. 296