program dragged_fall_plot
      real k,m
      write(*,*)'Vx=?'
      read(*,*)vx
      write(*,*)'Vy=?'
      read(*,*)vy
      write(*,*)'k=?'
      read(*,*)k
      write(*,*)'m=?'
      read(*,*)m
      x=0.0
      y=1.0
      dt=0.0001
      i_dummy=pgopen('/xwin')
      call pgenv(0.0,10.0,0.0,10.0,1,0)
      do while (y>=0)
        v=sqrt(vx**2+vy**2)
        fd=(-1)*k*v
        ax=fd*(vx/v)/m
        ay=-9.8+fd*(vy/v)/m
        vx=vx+ax*dt
        vy=vy+ay*dt
        x=x+vx*dt
        y=y+vy*dt
        call pgpt(1,x,y,20)
      end do
      write(*,*)'Distance=',x,';Hight=',y
      call pgclos
      end