Beispiel - SAM - ETH Zürich

AWP aus Bsp. 2.4.19
ode45 mit verschiedenen absoluten/relativen Toleranzen
Im Gegensatz zu Bsp. 2.6.8:
!
Toleranzen sagen nichts über globalen Fehler
Erklärung: wie in Bsp. 2.6.9 liegt ein schlecht konditioniertes
AWP vor, was den Einfluss von Einschrittfehlern
auf den Diskretisierungsfehler unkalkulierbar
macht.
Fehler |y h
(4)−y(4)|
10 1
10 0
10 −1
10 −2
10 −3
rtol = atol^(1/2)
rtol=atol^(2/3)
rtol=10*atol
10 −4
10 −6 10 −5 10 −4 10 −3 10 −2 10 −1
10 2 (Absolute) Toleranz
Numerische
Mathemtik
R. Hiptmair
rev 35327,
25. April
2011
2.6
p. 318