EXP-7 4th ORDER RK METHOD

def f(x,y):

z=y-x

return z

x0=float(input("Enter the initial value of x = "))

y0=float(input("Enter the initial value of y = "))

h=float(input("Enter the step size h = "))

n=int(input("Enter the number of iterations n = "))

count=1

while (count<=n):

k1=f(x0,y0)

k2=f(x0+h/2,y0+(h*k1/2))

k3=f(x0+h/2,y0+(h*k2/2))

k4=f(x0+h,y0+h*k3)

y1=y0+(h/6)*(k1+2*k2+2*k3+k4)

print("The value of k1 =",k1,"k2=",k2,"k3= ",k3,"k4=",k4)

print("The new value of y is",y1)

x0=x0+h

y0=y1

count=count+1