function f(c)
      implicit none
      real*8 c, f
      f = sqrt(c)
      return
      end