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))