program Trajectory use Trueacceleration implicit none integer:: i real*8 :: x0,y0,ispeed,t,a,vx1,vy1,vx2,vy2,px1,py1,px2,py2 100 format (2F12.3) 101 format (A12,A12,A12,A12,A12) t=0 x0=0 y0=0 print*, "please enter initial speed: " read (*,*) ispeed print*, "please enter initial angle: " read (*,*) a vx1=ispeed*cos(a*(3.14159265359/180)) vy1=ispeed*sin(a*(3.14159265359/180)) print 101, "time", "x", "y", "x_dot", "y_dot" print 101, "(sec)", "(m)","(m)","(m/s)","(m/s)" do while (a>181) call velocity_position(vx1,vx2,vy1,vy2,px1,px2,py1,py2) print 100, t, vx2 t=t+0.01 a=a+1 end do print *, x0, y0 end program Trajectory
We use cookies to provide and improve our services. By using our site, you consent to our Cookies Policy. Accept Learn more