Sources Esope | Cast3M

fau_p5
C FAU_P5    SOURCE    BECC      10/05/05    21:15:06     6674
      function fau_p5(M,coef,alpha)
C
      real*8 fau_p5, M, coef, alpha, fau_m1, fau_m2
     &     , f1, f2, f2m
C
      f1 = fau_m1(M,coef)
      f2 = fau_m2(M,coef)
      f2m = fau_m2(M,(-1.0D0 * coef))
      if (abs(M) .gt. 1.0D0) then
         fau_p5 = f1 / M
      else
         fau_p5 = coef * 2.0D0 - M
         fau_p5 = fau_p5 - (coef*16.0D0*alpha*M*f2m)
         fau_p5 = fau_p5 * f2
      endif
C
      return
      end