Runge-Kutta method: Difference between revisions

This equation has an exact solution:

:<math>y(t) = \tfrac{1}{16}(t^2 +4)^2</math>

;Task

Demonstrate the commonly used explicit &nbsp; [[wp:Runge–Kutta_methods#Common_fourth-order_Runge.E2.80.93Kutta_method|fourth-order Runge–Kutta method]] &nbsp; to solve the above differential equation.

* Solve the given differential equation over the range <math>t = 0 \ldots 10</math> with a step value of <math>\delta t=0.1</math> (101 total points, the first being given)

* Print the calculated values of <math>y</math> at whole numbered <math>t</math>'s (<math>0.0, 1.0, \ldots 10.0</math>) along with error as compared to the exact solution.

;Method summary

Starting with a given <math>y_n</math> and <math>t_n</math> calculate:

:<math>y_{n+1} = y_n + \tfrac{1}{6} (\delta y_1 + 2\delta y_2 + 2\delta y_3 + \delta y_4)</math>

:<math>t_{n+1} = t_n + \delta t\quad</math>

=={{header|Ada}}==