Subrutina POLN

Esta subrutina es la más sencilla y su papel consiste en calcular un polinomio de cuarto grado para la variable $\xi$ ó para $\eta$ según que el valor del índice IP sea 1 ó 2; en AZ se guarda el valor del polinomio, en $\xi$ y en BZ el de $\eta$.

C   ************************************************************************
C   *                                                                      *
C   *   LA SUBRUTINA "POLN" CALCULA UN POLINOMIO DE CUARTO GRADO           *
C   *   Y HACE UN INDICADOR IGUAL A UNO CUANDO ES NEGATIVO                 *
C   *                                                                      *
C   ************************************************************************
       SUBROUTINE  POLN  (IP,Y,A,B,C,D,E,)
       COMMON  AZ,BZ,X,IXIN,IETN,III,(600),JJJ(600),H,ALFA,PFI,G1,G2,Z1,Z2,
     @ EPP,EPM,GMI,GPL,JII(120)
       DIMENSION AA(5)
       AA(1)= A
       AA(2)= B
       AA(3)= C
       AA(4)= D
       AA(5)= E
       POL= AA(1)
       DO 3 I=2,5
       POL=POL*Y+AA(I)
  3    CONTINUE
       GO TO (4,5),IP
  4    AZ=POL
       IF (AZ.LT.0.) IXIN=1
       AZ= ABS(AZ)
       AZ=SQRT(AZ)
       GO TO 6
  5    BZ=POL
       IF (BZ.LT.0.) IETN=1
       BZ=ABS (BZ)
       BZ=SQRT(BZ)
  6    RETURN
       END

Como las expresiones de los momentos son en términos de la raíz del polinomio entonces preguntamos por el signo de éste para saber si la expresión de $p_\xi \;$ o $\; p_{\eta}$ es imaginaria. En cualquier caso se toma el valor de la raíz para ser utilizado en otro lado. La manera de calcular un polinomio se discutió ya cuando vimos la subrutina TURN del programa TWOC y no insistiremos en ese método. No hay más que explicar de esta subrutina.