Esope source files

   1 : C CWWFUD    SOURCE    CHAT      05/01/12    22:33:52     5004           
   2 : c
   3 :                  FUNCTION C_WWFUNCDER(T,E)
   4 : c n_d, d_d no estan multiplicadas por sin(3t) porq cancela 
   5 : c con las derivadas t_s
   6 :       IMPLICIT INTEGER(I-N)
   7 :       complex*16 c_wwfuncder
   8 :       complex*16 t,aux0,aux1,aux2,aux3,aux4,cost,n,d,n_d,d_d
   9 :       real*8 e
  10 :       cost=cos(t)
  11 : c      aux0=sin(t)
  12 :       aux0=1.D0/(4.D0*cost**2-1.D0)
  13 :       aux1=1.D0-e**2
  14 :       aux2=4.D0*aux1*cost**2
  15 :       aux3=2.D0*e-1.D0
  16 :       aux4=sqrt(aux2+5.D0*e**2-4.D0*e)
  17 :       n   =aux3*aux3+aux2
  18 :       d   =2.D0*aux1*cost+aux3*aux4
  19 :       n_d =-8.D0*aux1*cost*aux0
  20 :       d_d =-(2.D0+4.D0*aux3*cost/aux4)*aux1*aux0
  21 :       c_wwfuncder=n_d/d-n*d_d/d**2
  22 :       return
  23 :       end
  24 :  
  25 :  
  26 :