Discusión del programa general

El programa consiste en llamados a las diferentes subrutinas que ya hemos discutido y además, realiza algunos cálculos auxiliares, por lo cual no hay ninguna dificultad en comprender su funcionamiento. Es en el programa principal donde se declaran los archivos de disco y las variables comunes. Ya se indicó el papel que desempeñan algunos índices como IT, IC, IG en el sentido de señalar al programa si deben realizarse determinados cálculos o de qué manera deben trabajar las subrutinas cuando existen varias opciones. En la siguiente hoja se podrá ver todo el proceso de cálculo mediante un diagrama de flujo y un listado del programa principal. No hemos incluido el programa PRYPO porque preferimos discutirlo en un apéndice aparte. No se ha atendido aquí nada con respecto a los detalles del lenguaje FORTRAN, se ha discutido lo que se consideró necesario para entender los algoritmos mediante los cuales obtenemos los diferentes aspectos del problema. Los detalles de FORTRAN pueden verse en los libros de Bork [12] y Louden [13]. Todas las subrutinas en FORTRAN discutidas en este capítulo están escritas en la versión de la 1130 para la cual fueron planeadas originalmente, aunque a últimas fechas se adaptó el programa para la computadora PDP-10. Más adelante se presentan otros programas para la PDP-10 y se aclara en cada caso como han sido planeados.

//DELETE                       TWOC
// FOR
* IOCS(CARD,1132 PRINTER,TYPEWRITER,KEYBOARD,DISK)
* ONE WORD INTEGERS
* LIST ALL
* NAME TWOC
C
C      PROGRAMA PRINCIPAL
C
       DIMENSION   X(6),DX(6)
       DIMENSION   III(100)
       COMMON      T,DT,H,ALFA,IPR,IPO,G1,G2,Z1,Z2,GPL,GMI,EPP,EPM,XK,MO
       COMMON      IC,IT,IG,IP
       EQUIVALENCE (PPHI,X(6))
       DEFINE FILE 10(50,120,U,K10)
       DEFINE FILE 20(50,120,U,K20)
       MO=1
     1 CONTINUE
       CALL INDI (X)
       IC=IC+1
       IT=IT+1
       IG=IG+1
       IP=IP+1
       GPL=G1+G2
       GMI=G1-G2
       EPP=Z1+Z2
       EPM=Z1-Z2
       XK=2.0*G1*G2
       CALL GRAD (X,DX)
       GO TO (2,3),IT
     2 CONTINUE
       CALL TURN (H,ALFA,PPHI,EPP,EPM,GPL,GMI,XK)
     3 CONTINUE
       GO TO (4,5),IC
     4 CONTINUE
       CALL CNTU (G1,G2,Z1,Z2)
     5 CONTINUE
       CALL INDO (X)
       CALL PAGE (X,1)
       DO 20 J=1,IPO
       CALL GRAF (III,X,1)
       DO 10 I=1,IPR
       CALL RUKU (X,DX)
       CALL AUGV (X,X,DT,DX)
       GO TO (6,7),IG
     6 CONTINUE
       CALL GRAF (III,X,2)
     7 CONTINUE
    10 CONTINUE
       CALL GRAF (III,X,3)
       CALL PAGE (X,3)
    20 CONTINUE
       CALL GRAF (III,X,4)
       CALL PAGE (X,5)
       GO TO 1
       END