To calculate result you have to disable your ad blocker first.
Runge-Kutta Method Calculator
To use Runge-Kutta method calculator, enter the equation, write initial condition, interval values, step size, and click calculate button
Runge-Kutta Method Calculator
Runge-Kutta Method Calculator is used to calculate ordinary differential equations (ODEs) numerically.
What is Runge-kutta method?
Runge-kutta method is a numerical technique to find the solution of the ordinary differential equation. The method provides us with the value of y as the corresponding value of x. It is also known as the RK method. Here we use the most common RK-4 which is known as runge-kutta of order fourth.
Formula of Runge-kutta method
The formula of the Runge-Kutta method of order fourth is given as:
yi+1=yi+h/6(k1+2k2+2k3+k4) where
xi+1 = xi+h
where i is the number of iterations.
For the value of k’s, we have
- k1= f(xi,yi)
- k2= f(xi+h/2, yi+hk1/2)
- k3=f(xi+h/2,yi+hk2/2)
- k4= f(xi+h,yi+hk3)
How to calculate Runge-kutta method problems?
Here is a numerical example to understand this iterative method.
Example:
If we have differential function x2+1 with the "2" number of iterations. The initial value of the x is 1 to 3 and the initial value of y is 2.
Solution:
Step 1: The step size is
h= (3-1)/2 = 1, x0= 1, y0 = 2 and n= 2
f(x, y) = x^2+1
Step 2: For 1st iteration
x1= x0+h = 1+1 = 2
k1= f(x0,y0)= f(1,2) = 2
k2= f(x0+h/2, y0+hk1/2) = f(1+1/2,2+(1)(2)/2) = f(3/2,3) = 3.25
k3= f(x0+h/2, y0+hk2/2) = f(1+1/2,2+(1)(3.25)/2) = f(2/2,3.625) = 3.25
k4 = f(x0+h,y0+hk3) = f(1+1,2+(1)(3.25)) = f(2,5.25) = 5
y1 = y(x1) = y(2) = y0+ h/6(k1+2k2+2k3+k4) = 2+1/6(2+2(3.25)+2(3.25)+5) = 5.333
Step 3: For the 2nd iteration we have x1 and y1 to calculate x2 and y2.
x2= x1+h = 2+1 = 3
k1= f(x1,y1) = f(2,5.333) = 5
k2= f(x1+h/2, y1+hk1/2) = f(2+1/2,5.333+(1)(5)/2) = f(5/2,7.83) = 7.25
k3= f(x1+h/2, y1+hk2/2) = f(2+1/2,5.333+(1)(7.25)/2= f(5/2,8.958) = 7.25
k4 = f(x1+h,y1+hk3) = f(2+1,5.333+(1)(7.25)) = f(3,12.583) = 10
y2 = y(x2) = y (3) = y1+ h/6(k1+2k2+2k3+k4) = 5.333+1/6(5+2(7.25) +2(7.25) +10) = 12.666
Hence the approximate value of y2 = 12.666.