part 2
write (unit=*,fmt="(2x'dist to target ? ')",advance="no")
write (unit=*,fmt="(2x'Init veloc ')",advance="no")
read (unit=*,fmt=*) Init_vel
write (unit=*,fmt="(2x'Init angle ')",advance="no")
read (unit=*,fmt=*) Angle
Ang_cos=cos(Angle*Rad_cnv)
Ang_sin=sin(Angle*Rad_cnv)
Time=Init_vel*Ang_sin/(0.5*Grav)
! Displ=2*Init_vel**2*Ang_sin*Ang_cos/Grav
Displ=Init_vel*Ang_cos*Time