matlab
12345678910111213141516171819202122232425262728293031323334353637
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))
Python
INFO