Sources Esope | Cast3M

cwwfud
C CWWFUD    SOURCE    CB215821  25/04/24    21:15:06     12248          
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=CMPLX(1.D0,0.D0)/(CMPLX(4.D0,0.D0)*cost**CMPLX(2.D0,0.D0)-
     &     CMPLX(1.D0,0.D0))
      aux1=CMPLX(1.D0,0.D0)-e**CMPLX(2.D0,0.D0)
      aux2=CMPLX(4.D0,0.D0)*aux1*cost**CMPLX(2.D0,0.D0)
      aux3=CMPLX(2.D0,0.D0)*e-CMPLX(1.D0,0.D0)
      aux4=sqrt(aux2+CMPLX(5.D0,0.D0)*e**CMPLX(2.D0,0.D0)-
     &     CMPLX(4.D0,0.D0)*e)
 
      n   =aux3*aux3+aux2
      d   =CMPLX(2.D0,0.D0)*aux1*cost+aux3*aux4
      n_d =CMPLX(-8.D0,0.D0)*aux1*cost*aux0
      d_d =-(CMPLX(2.D0,0.D0)+CMPLX(4.D0,0.D0)*aux3*cost/aux4)*aux1*aux0
      c_wwfuncder=n_d/d-n*d_d/d**CMPLX(2.D0,0.D0)
      return
      end