subroutine func(t, ft, dft)
c
c -- the solution to f(t) = a + b/t + c log(t) - log(.01)
c
      real t, ft, dft, a, b, c, logp
      data a/18.19/, b/-23180.0/, c/-0.8858/ logp/-4.60517/
c
      ft = a + b / t + c * alog(t) - logp
      dft = -b / (t*t) + c / t
      return
      end