The Euler method uses the first order Taylor series to find the next value of the ordinary differential equation. Therefore, this method is not very accurate. To achieve higher accuracy, Runge-Kutta method employs higher order terms of the Taylor series in its approximation. For example the second order Runge-Kutta method uses the Taylor series up to he second order term. where, Using the mathematical relations, it can be demonstrated that the above equation can be written in the following form: where,
Aand_{1}, A_{2}, P_{1} Q are some constant which can be obtained from following system of equations:_{11}It is clear that the number of unknowns are more than the number of equations. Therefore, a series of solution can be found for these constants. For
A ½ and _{1} = A_{2} =P the general Runge-Kutta method is transformed to the following form:_{1} = Q_{11} = 1This method is known as Heun’s method with single corrector. For
A _{1} =0, A_{2} =1 and P the general Runge-Kutta method can be written as follows:_{1} = Q_{11} = ½ This method is known as improved polygon method. For
A _{1} =1/3, A_{2} =2/3 and Pthe Ralston’s method can be obtained:_{1} = Q_{11} = ¾, The order of all above methods is
x_{n}.
(Source Code in C++) |