The program

What follows is the entire program — about a thousand lines of fixed-form FORTRAN IV, the same deck that runs when you press SUBMIT on the front page. It is printed the way it would come off the printer, columns and all; the notes in the margin say what each part is doing and where it comes from in the 1964 paper. Begin with the flow of the calculation, then read the listing itself. The boxes of the flow chart are links into the code.

2–1The flow of the calculation

Fig. 7. The logical pattern of the calculation, after Fig. 2 of Henyey, Forbes & Gould (1964). The left column is one Newton–Raphson iteration; the right column advances the evolutionary sequence. Every box is a link into the listing below.

2–2The listing

Fixed form: columns 1–5 hold statement labels (printed red), column 6 marks a continuation, columns 7–72 carry the statement. Lines whose number is red carry a note in the margin; on a narrow screen the notes move inline. The deck compiles unchanged with the repository's 7094-dialect compiler and with any modern Fortran compiler in legacy mode.

fetching the deck…

← back to the machine the repository the 1964 paper