Sources Esope | Cast3M

wwfude
C WWFUDE    SOURCE    CB215821  16/04/21    21:18:41     8920
c
             FUNCTION 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)
      real*8 wwfuncder
      real*8 t,e,aux0,aux1,aux2,aux3,aux4,cost,n,d,n_d,d_d
      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**2+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
      wwfuncder=n_d/d-n*d_d/d**2
      return
      end