Descripción del programa

Una vez expuesta la naturaleza del problema y del método de integración para las ecuaciones del movimiento, el siguiente paso es discutir la parte computacional; todos los procesos de cálculo que requiere el problema han sido elaborados mediante el uso de un programa en FORTRAN para una computadora IBM 1130 (actualmente se tiene también en gran parte dicho programa en ALGOL para una máquina Burroughs-5500) existe también una versión para la computadora PDP-10. La manera más comoda y práctica de elaborar un programa con una serie de opciones y diversidad de cálculos, es la de dividirlo en subrutinas para realizar las diferentes tareas; algunas de las subrutinas usan variables comunes a todas ellas; para eso, hay que declararlas en común, por medio de la instrucción COMMON. En seguida daremos una lista de dichas variables:

T	Representa al tiempo
DT	Representa al incremento de tiempo
H	Contiene a la hamiltoniana (energía)
ALFA	Es la constante $\alpha$, que ya discutimos
IPR	Indica el número de pasos de integración que hay que
	efectuar antes de imprimir un renglón
IPO	Es el número de renglones deseados para cada gráfica
G1, G2	Representan las cargas magnéticas que hemos colocado en
	los dos centros
Z1, Z2	Son las cargas eléctricas. En las fórmulas que usa el
	programa, las cargas no aparecen aisladas sino que
	formando diferentes combinaciones tales como:
GPL = G1 + G2	correspondiente a $g_+$
GMI = G1 - G2	correspondiente a $g_-$
EPP = Z1 + Z2	correspondiente a $\varepsilon_+$
EPM = Z1 - Z2	correspondiente a $\varepsilon_-$
XK = $2\cdot G1\cdot G2$	correspondiente a $\kappa$. Hemos visto que $\kappa=\frac{1}{2}(g_+^2-g_-^2)$
MO	Ya no se usa en este programa y puede omitirse
IC, IP, IT, IG	Son indicadores que tienen la función de comunicar al
	programa si se desea o no realizar determinadas gráficas

El programa utiliza, además, en varias subrutinas, los arreglos III, y X; el primero contiene la imagen de una línea y el segundo contiene las coordenadas (posición y momento) de la partícula en el espacio fase. Antes de ver cómo funciona el programa, es conveniente discutir por separado cada una de las diferentes subrutinas, y después estudiar cómo van ensambladas en lo que es el programa mismo.