Define substitutions in terms of fundamental variables; add list of constants
Polish code output, and finish precession equations