subroutine derivs(n, t, y, ff)
!use work
! Calculeaza functia din BDF
implicit real(8) (a-h,o-z)
real(8) y(n), ff(n)

! y =xa(1,:)
! Membrii 2:
!      Fcn(1) = 1.0d0/y(1) - x**2 - 2.0d0/x**3
!      Fcn(2) = y(1)/y(2)**2 - 1.0d0/x - 0.5d0/SQRT(x**3)

ff(1) =1.0d0/y(1) - t**2 - 2.0d0/t**3
ff(2) =y(1)/y(2)**2 - 1.0d0/t - 0.5d0/dSQRT(t**3)


END subroutine derivs

