Sources Esope | Cast3M

cwwfun
C CWWFUN    SOURCE    CB215821  25/04/24    21:15:07     12248          
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=CMPLX(4.D0,0.D0)*aux1*cost**CMPLX(2)
      aux3=CMPLX(2.D0,0.D0)*e-CMPLX(1.D0,0.D0)
      c_wwfunc=(aux3**CMPLX(2)+aux2)/
     .   (CMPLX(2.D0,0.D0)*aux1*cost+aux3*sqrt(aux2+CMPLX(5.D0*e**2)-
     .    CMPLX(4.D0*e)))
      return
      end