numabj.nb |
Integroimalla differentiaaliyhtälö puolittain välin
yli saadaan
![[Graphics:Images/numabj_gr_3.gif]](Images/numabj_gr_3.gif)
Adamsin -- Bashforthin menetelmässä integraalille lasketaan approksimaatio korvaamalla funktio kolmannen asteen interpolaatiopolynomilla. Tämän tukipisteinä (so. pisteinä, joiden kautta polynomin kuvaaja kulkee) ovat neljä edellistä jo laskettua pistettä
,
.
![[Graphics:Images/numabj_gr_7.gif]](Images/numabj_gr_7.gif)
![[Graphics:Images/numabj_gr_9.gif]](Images/numabj_gr_9.gif)
Interpolaatiopolynomin muodostamisessa tarvittavat pisteet ovat
![[Graphics:Images/numabj_gr_11.gif]](Images/numabj_gr_11.gif)
Interpolaatiopolynomi saadaan yhdellä komennolla (muuttujana ):
![[Graphics:Images/numabj_gr_14.gif]](Images/numabj_gr_14.gif)
Integraalin approksimaatio saadaan integroimalla polynomi:
![[Graphics:Images/numabj_gr_16.gif]](Images/numabj_gr_16.gif)