C CWWFUD SOURCE CHAT 05/01/12 22:33:52 5004 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 complex*16 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/(4.D0*cost**2-1.D0) aux1=1.D0-e**2 aux2=4.D0*aux1*cost**2 aux3=2.D0*e-1.D0 aux4=sqrt(aux2+5.D0*e**2-4.D0*e) n =aux3*aux3+aux2 d =2.D0*aux1*cost+aux3*aux4 n_d =-8.D0*aux1*cost*aux0 d_d =-(2.D0+4.D0*aux3*cost/aux4)*aux1*aux0 c_wwfuncder=n_d/d-n*d_d/d**2 return end