Sources Esope | Cast3M

cwwfun
C CWWFUN    SOURCE    KK2000    14/04/09    21:15:14     8027
c
                 FUNCTION C_WWFUNC(T,E)
      IMPLICIT INTEGER(I-N)
 
      complex*16 t,aux1,aux2,aux3,cost,c_wwfunc
 
      real*8 e
 
      cost=cos(t)
      aux1=CMPLX(1.D0-e**2)
      aux2=(4.D0,0.D0)*aux1*cost**CMPLX(2)
      aux3=(2.D0,0.D0)*e-(1.D0,0.D0)
      c_wwfunc=(aux3**CMPLX(2)+aux2)/
     .   ((2.D0,0.D0)*aux1*cost+aux3*sqrt(aux2+CMPLX(5.D0*e**2)-
     .    CMPLX(4.D0*e)))
      return
      end