(см. скан)
Затем программа выполняется 
раз, причем на каждом шаге по времени вызываются подпрограммы 
и 
(именно в указанном порядке). Детали каждого обращения к подпрограммам указаны ниже. Вся последовательность действий аналогична блок-схеме на рис. 2.3, только операция взвешивания перенесена в подпрограммы 
и 
(см. скан)
После прохождения 
шагов по времени выходим из основного цикла через процедуры 
и 
(см. скан)
Остальные подробности, относящиеся к каждому шагу всех подпрограмм, приведены в последующих параграфах, в основном в порядке их использования в программе.
Есть несколько особенностей в 
которые обусловлены спецификой развитых в Ливерморе операционных систем для ЭВМ CDC-7600 и CRAY-1. Большинство из них относится к подпрограммам 
и подпрограммам, рисующим графики. В программу интегрирования уравнений движения включены специальные операторы, позволяющие использовать возможность векторных вычислений на компьютере CRAY-1. Кроме того, в этой программе использован предкомпилятор для удобства изменения размерности всех массивов и для повторения COMMON блоков в подпрограммах.
Операторы 
определяют величины, которые являются фиксированными во время компилирования и сохраняются постоянными (по сравнению с операторами 
которые присваивают начальные значения переменным, изменяющимся в процессе выполнения). CLICHE и ENDCLICHE ограничивают операторы, которые должны вставляться в любой источник, цитируемый оператором 
Таким образом, изменения в размерностях массивов и COMMON блоков необходимо делать только в одном месте.
Далее следует текст основной программы:
(см. скан)
Следом вызываются подпрограммы 
и 
(которая в свою очередь вызывает процедуру прямого и обратного преобразований Фурье) и 
вызывающая 
Обращение к подпрограммам происходит в указанном здесь порядке. Эти процедуры используются только 1 раз для задания соответствующих начальных условий, т. е. величин 
для всех частиц, а также 
в качестве начальных условий для основного вычислительного цикла. Все действия выполняются в следующем порядке:
