Yhtälö y′ = f(x, y)
Kirjoita yhtälön oikea puoli muuttujilla x ja y. Käytä merkkejä + − * / ja ^ sekä funktioita sin, cos, exp, ln, sqrt. Esimerkiksi "y", "x + y" tai "x*y".
Ratkaise ensimmäisen kertaluvun differentiaaliyhtälön alkuarvotehtävä y′ = f(x, y) numeerisesti Eulerin tai Runge–Kutta-menetelmällä.
Kirjoita yhtälön oikea puoli muuttujilla x ja y. Käytä merkkejä + − * / ja ^ sekä funktioita sin, cos, exp, ln, sqrt. Esimerkiksi "y", "x + y" tai "x*y".
Moni ensimmäisen kertaluvun differentiaaliyhtälö ei ratkea suljetussa muodossa, joten ratkaisu lasketaan numeerisesti. Tämä laskuri ratkaisee alkuarvotehtävän y′ = f(x, y), y(x0) = y0 etenemällä alkupisteestä pienin askelin ja arvioimalla y:n arvon yhä uusissa pisteissä. Voit valita yksinkertaisen Eulerin menetelmän tai tarkemman neljännen kertaluvun Runge–Kutta-menetelmän.
Tehtävä koostuu yhtälöstä ja yhdestä tunnetusta pisteestä:
y′ = f(x, y), y(x0) = y0
Etsitään ratkaisukäyrä y(x), joka kulkee pisteen (x0, y0) kautta. Numeerinen menetelmä korvaa jatkuvan käyrän pisteiden jonolla, jossa peräkkäisten x-arvojen väli on askelpituus h.
Yksinkertaisin menetelmä ottaa askeleen kulmakertoimen suuntaan välin alussa:
yₙ₊₁ = yₙ + h · f(xₙ, yₙ)
Menetelmä on helppo ymmärtää, mutta virhe kasvaa askelpituuden myötä.
RK4 laskee neljä kulmakertoimen arviota askeleen sisällä ja yhdistää ne:
k₁ = f(xₙ, yₙ)
k₂ = f(xₙ + h/2, yₙ + h·k₁/2)
k₃ = f(xₙ + h/2, yₙ + h·k₂/2)
k₄ = f(xₙ + h, yₙ + h·k₃)
yₙ₊₁ = yₙ + h/6 · (k₁ + 2k₂ + 2k₃ + k₄)
Tulos on selvästi tarkempi kuin Eulerin, koska se ottaa huomioon kulmakertoimen muutoksen askeleen aikana.
Ratkaistaan y′ = y alkuarvolla y(0) = 1. Tarkka ratkaisu on y = eˣ.
Esimerkki näyttää, kuinka RK4 saavuttaa paljon paremman tarkkuuden samalla askelpituudella.