import math


# constants

#baseball
m = 0.15 #kg
D = 0.07  #m
gamma = 0.25 #  N s^2/m^2
c = gamma * D**2

g = 9.8

print(math.sqrt(m*g/c))

# Initial Conditions

x = 0.0
y = 0.0

vx = 19.3
vy = 23.0


def   dvdt(vxp,vyp):
  v = math.sqrt(vxp*vxp + vyp*vyp)
  dvxdt= - c/m*v*vxp
  dvydt = -g  - c/m*v*vyp
  return dvxdt,dvydt

dt = 0.1 
t = 0
print(t,vx,vy)

for i in range(100):
   dvxdt,dvydt = dvdt(vx,vy) 
   vx=vx+dvxdt*dt
   vy=vy+dvydt*dt
   t=t+dt
   x=x+vx*dt
   y=y+vy*dt
   print(t,vx,vy,x,y)