Subrutina PRYPO

Esta subrutina consta de dos partes, a saber: lo que podría llamarse el programa principal y una subrutina auxiliar que hemos llamado CALCA la cual tiene como tarea calcular dos funciones y cargar los símbolos que las representan, en los arreglos IARR y JARR, de acuerdo con lo indicado por ciertos parámetros que calcula el programa principal. Las funciones en cuestión son las discutidas en la sección (I.5.1). El programa principal decide bajo qué condiciones ha de operar la subrutina, con qué parámetros, en qué arreglo hay que cargar un símbolo, preve que en la vecindad de los puntos singulares no se calcule una función, esto es para evitar un sobreflujo por división por una cantidad cercana a cero. PRYPO tiene como argumentos (IR, K, GPL, GMI, H, ALFA, EPP, EPM) y la subrutina auxiliar CALCA utiliza esos mismos parámetros pero se les da nombres más generales porque deseamos obtener las curvas de la variable $\xi$ y las de $\eta$ mediante el mismo procedimiento. CALCA tiene tres parámetros adicionales: el primero es un número entero que toma dos valores, el valor 1 indica que hay que almacenar los puntos de la variable y el valor 2 va a indicar que los puntos que se almacenarán corresponden a la variable ?. El segundo argumento indica la columna donde se va a situar el símbolo, el renglón corresponde a la parte entera del valor de la función, habiendo un punto al pie de la gráfica o en la parte superior cuando la curva sale de la hoja. El tercer parámetro es también un número entero que puede tomar el valor 1 ó el valor 2, decide si hay que omitir o realizar el cálculo de las funciones $f_{2}$ y $s_{2}$ ( página 15) en los puntos donde ocurren las singularidades para dichas funciones. Eso lo decide el programa PRYPO al comparar $p_{\phi}$ con los parámetros $-g_{+}$ y $-g_{-}$. La subrutina PRYPO tiene un ciclo controlado por el índice k que permite variar ligeramente alguno de los parámetros que intervienen en el cálculo de nuestras funciones y obtener nuevamente las gráficas para cada variación; con ello se obtendrán las familias de curvas que discutimos en la sección (I.5.1). Esta subrutina puede usarse independientemente del programa principal y se tiene una versión en fortran para la computadora 3150 y en algol para la B-5500, es más ventajoso usarla por separado por que hay mayor libertad en el manejo de los parámetros. El inconveniente es que al hacerlo así no se conocen los valores iniciales. Actualmente, sabiendo que los valores de los momentos estan tratados por las expresiones (I-4-19) y (I-4-20) se ha elaborado una subrutina auxiliar de PRYPO que calcula el valor de dichas expresiones en cualquier punto dado, eso permite invertir el proceso de derivación de los diferentes parámetros del problema y establecer así condiciones iniciales para valores de las constantes seleccionadas de antemano. La versión más completa de PRYPO está escrita en FORTRAN para la computadora PDP-10; en el Apéndice 1, se discute con más detalle.