matlab

x=0.0
y=5.0
z=10.0

h=0.1

print("\n"*40)
print("x","y")
print (x,y)

a=[x]
b=[y]

while(x<2.0):
    k1=z
    l1=60-5*y-4*z
    k2=z+(h/2.)*l1
    x_2=x+(h/2.)
    y_2=y+(h/2.)*k1
    z_2=k2
    l2=60-5*y_2-4*z_2
    k3=z+(h/2.)*l2
    x_3=x+(h/2.)
    y_3=y+(h/2.)*k2
    z_3=k3
    l3=60-5*y_3-4*z_3
    k4=z+(h/2.)
    x_4=x+(h/2.)
    y_4=y+(h/2.)*k3
    z_4=k4
    l4=60-5*y_4-4*z_4
    x=x+h
    round(x,2)
    y=y+(h/6.)*(k1+(2.*k2)+(2.*k3)+k4)
    z=z+(h/6.)*(l1+(2.*l2)+(2.*l3)+k4)
    
    print(round(x,6),round (y,6))