program fly
      real m

      write(*,*)'number of gun powder?'
      read(*,*) num_gun_powder
      write(*,*)'shooting angle?'
      read(*,*) angle
      angle = angle/180 * 3.14
      write(*,*)'mass?'
      read(*,*) m

c     Ek = (1/2)m*V*V
      Ek = 1.0 * num_gun_powder
      v_ini = sqrt(2*Ek/m)

      dt = 0.001
      vx = v_ini*cos(angle)
      vy = v_ini*sin(angle)
      x = 0.0
      y = 1.0
      g = -9.8
      do while(y >= 0.0)
        ax = 0
        ay = g
c       ax = a_drag_x
c       ay = g + a_drag_y
        vx = vx + ax*dt
        vy = vy + ay*dt
        x = x + vx*dt
        y = y + vy*dt
      end do
      write(*,*)'final hight is:',y
      write(*,*)'final distance is:',x



      end