Operation Manual

1036
De methode Runge-Kutta
De methode Runge-KuttaDe methode Runge-Kutta
De methode Runge-Kutta
Voor integraties van gewone differentiaalvergelijkingen met de Runge-Kutta methode,
gebruikt de
TI-89 Titanium / Voyage™ 200 de Bogacki-Shampine 3(2) formule, zoals die
beschreven stond in het tijdschrift Applied Math Letters, 2 (1989), pag. 1–9.
De Bogacki-Shampine 3(2) formule
De Bogacki-Shampine 3(2) formuleDe Bogacki-Shampine 3(2) formule
De Bogacki-Shampine 3(2) formule
De Bogacki-Shampine 3(2) formule biedt een resultaat met een nauwkeurigheid van de
derde orde en een foutschatting die gebaseerd is op een ingebedde formule van de
tweede orde. Voor een probleem in de vorm:
y' = ƒ(x, y)
en een gegeven stapgrootte h, kan de Bogacki-Shampine formule als volgt geschreven
worden:
F
1
= ƒ(x
n
, y
n
)
F
2
= ƒ (x
n
+ h , y
n
+ h F
1
)
F
3
= ƒ (x
n
+ h , y
n
+ h F
2
)
y
n+1
= y
n
+ h ( F
1
+ F
2
+ F
3
)
x
n+1
= x
n
+ h
F
4
= ƒ (x
n+1
, y
n+1
)
errest = h ( F
1
ì F
2
ì F
3
+ F
4
)
De foutschatting errest wordt gebruikt om de stapgrootte automatisch te controleren. Voor
een diepgaande bespreking van hoe dit gedaan kan worden, zie Numerical Solution of
Ordinary Differential Equations door L. F. Shampine (New York: Chapman & Hall,
1994).
De software van de TI-89 Titanium / Voyage™ 200 past de stapgrootte niet aan om uit te komen bij
specifieke uitvoerpunten. Hij neemt, eerder, de grootste stappen die hij kan nemen (gebaseerd op
de fouttolerantie diftol) en verkrijgt resultaten voor x
n
_{ x {_ x
n+1
gebruik makend van het
derdegraads interpolatie-polynoom dat gaat door het punt
(x
n
, y
n
) met richtingscoëfficiënt F
1
,
en door het punt
(x
n+1
, y
n+1
) met richtingscoëfficiënt F
4
. De interpolant is efficiënt en biedt
resultaten tijdens de stap die even betrouwbaar zijn als de resultaten aan het eind van de stap.
1
2
---
1
2
---
3
4
---
3
4
---
2
9
---
1
3
---
4
9
---
5
72
------
1
12
------
1
9
---
1
8
---