Sources Esope | Cast3M

wwfunc
C WWFUNC    SOURCE    CB215821  16/04/21    21:18:42     8920
CCC
C **********************************************************************
CCC
      FUNCTION WWFUNC(T,E)
      IMPLICIT INTEGER(I-N)
      real*8 t,e,aux1,aux2,aux3,cost,wwfunc
      cost=COS(t)
      aux1=1.D0-e**2
      aux2=4.D0*aux1*cost**2
      aux3=2.D0*e-1.D0
      wwfunc=(aux3**2+aux2)/
     .       (2.D0*aux1*cost+aux3*SQRT(aux2+5.D0*e**2-4.D0*e))
      return
      end