subroutine tiro_parabolico(neq,t,y,ydot)

  implicit none

!!$  Ecuacion del tiro parabolico con rozamiento
!!$  en dos dimensiones
!!$  Se hacen cosas raras con las variables para 
!!$  adaptarse al solver slsode

  integer,intent(in) :: neq
  real,intent(in) :: t
  real,dimension(4),intent(in) :: y
  real,dimension(4),intent(out) :: ydot

!!$  constantes del problema
!!$    aceleracion de la gravedad
  real,parameter :: g=9.8
!!$    coeficiente de rozamiento
  real,parameter :: c=0.01

  ydot(1)=y(3)
  ydot(2)=y(4)
  ydot(3)=-c*y(3)*sqrt(y(3)**2+y(4)**2)
  ydot(4)=-g-c*y(4)*sqrt(y(3)**2+y(4)**2)

end subroutine tiro_parabolico