Esope source files

   1 : C CWWFUN    SOURCE    CHAT      05/01/12    22:33:58     5004           
   2 : c
   3 :                  FUNCTION C_WWFUNC(T,E)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       complex*16 t,aux1,aux2,aux3,cost,c_wwfunc
   6 :       real*8 e
   7 :       cost=cos(t)
   8 :       aux1=1.D0-e**2
   9 :       aux2=4.D0*aux1*cost**2
  10 :       aux3=2.D0*e-1.D0
  11 :       c_wwfunc=(aux3**2+aux2)/
  12 :      .       (2.D0*aux1*cost+aux3*sqrt(aux2+5.D0*e**2-4.D0*e))
  13 :       return
  14 :       end
  15 :  
  16 :  
  17 :