Subrutina GRAD

Aquí es donde se calculan las diferentes expresiones para las ecuaciones de movimiento dadas por (I-4-9) a (I-4-14); se encuentra el valor del hamiltoniano (con potencial repulsivo o sin él) así como la constante de separación. En las expresiones que aparecen en el listado los cálculos están hechos sin considerar el potencial repulsivo. Cuando se quieran ver los efectos del potencial repulsivo habrá que modificar ligeramente las expresiones pero esencialmente la subrutina conserva su estructura.

      SUBROUTINE  GRAD (Y,DY)
      DIMENSION   X(6),DX(6),Y(6),DY(6)
      COMMON      T,DT,H,ALFA,IPR,IPO,G1,G2,Z1,Z2,GPL,GMI,EPP,EPM,XK,MO
      COMMON      IC,IT,IG,IP
      EQUIVALENCE (XI,X(1)),(ETA,X(2)),(PHI,X(3))
      EQUIVALENCE (PXI,X(4)),(PETA,X(5)),(PPHI,X(6))
      EQUIVALENCE (DXI,DX(1)),(DETA,DX(2)),(DPHI,DX(3))
      EQUIVALENCE (DPXI,DX(4)),(DPETA,DX(5)),(DPPHI,DX(6))
C
      CALL CPYV (X,Y)
      D   = XI*XI-ETA*ETA
      XXI =XI*XI-1
      EETA=1.-ETA*ETA
      PPHX=(PPHI+GMI*XI)/XXI
      PPHE=(PPHI+GPL*ETA)/EETA
      HX  =PXI*PXI+PPHX*PPHX
      HE  =PETA*PETA+PPHE*PPHE
      H   =(XXI*HX+EETA*HE+EPP*XI-EPM*ETA+XK)/D
      ALFA=XI*XI*H-XXI*HX-EPP*XI
      DXI =PXI*XXI/D
      DETA=PETA*EETA/D
      DPHI=(PPHX+PPHE)/D
      DPXI=(XI*H-XI*PXI*PXI-PPHX*(GMI-PPHX*XI)-EPP/2.)/D
      DPETA=(ETA*PETA*PETA-PPHE*(GPL+PPHE*ETA)+EPM/2.-ETA*H)/D
      DPPHI=0.0
C
      CALL CPYV (DY,DX)
      RETURN
      END

Los argumentos de GRAD son (Y, DX) y su funcionamiento es muy simple después de haber visto la parte analítica del problema. Esta subrutina como hemos dicho, prepara las ecuaciones de movimiento que son resueltas numéricamente por otra subrutina que discutiremos enseguida.