Sources Esope | Cast3M

cwwfud
C CWWFUD    SOURCE    KK2000    14/04/09    21:15:13     8027
c
                 FUNCTION C_WWFUNCDER(T,E)
c n_d, d_d no estan multiplicadas por sin(3t) porq cancela
c con las derivadas t_s
      IMPLICIT INTEGER(I-N)
 
      complex*16 c_wwfuncder,t,aux0,aux1,aux2,aux3,aux4,cost,n,d,n_d,d_d
 
      real*8 e
 
      cost=cos(t)
c      aux0=sin(t)
      aux0=(1.D0,0.D0)/((4.D0,0.D0)*cost**(2.D0,0.D0)-(1.D0,0.D0))
      aux1=(1.D0,0.D0)-e**(2.D0,0.D0)
      aux2=(4.D0,0.D0)*aux1*cost**(2.D0,0.D0)
      aux3=(2.D0,0.D0)*e-(1.D0,0.D0)
      aux4=sqrt(aux2+(5.D0,0.D0)*e**(2.D0,0.D0)-(4.D0,0.D0)*e)
 
      n   =aux3*aux3+aux2
      d   =(2.D0,0.D0)*aux1*cost+aux3*aux4
      n_d =(-8.D0,0.D0)*aux1*cost*aux0
      d_d =-((2.D0,0.D0)+(4.D0,0.D0)*aux3*cost/aux4)*aux1*aux0
      c_wwfuncder=n_d/d-n*d_d/d**(2.D0,0.D0)
      return
      end