! Ecuatia-1 a sagetii: f=x*l
subroutine SUB(fun, dfun, x)
common B, C
real k


fun =B/x -phi(x) -C

dfun =-B/(x*x) -2*k(x)

END
 