# Copyright: Public domain.
# Filename: CONIC_SUBROUTINES.agc
# Purpose: Part of the source code for Luminary 1A build 099.
# It is part of the source code for the Lunar Module's (LM)
# Apollo Guidance Computer (AGC), for Apollo 11.
#
# Assembler: yaYUL
# Contact: Jim Lawton <jim DOT lawton AT gmail.com>
# Website: www.ibiblio.org/apollo.
# Pages: 1159-1204
# Mod history: 2009-05-28 JL Started updating from page images.
# 2009-06-02 JL Finished updating from page images.
# 2010-12-31 JL Fixed page number comments.
#
# This source code has been transcribed or otherwise adapted from digitized
# images of a hardcopy from the MIT Museum.  The digitization was performed
# by Paul Fjeld, and arranged for by Deborah Douglas of the Museum.  Many
# thanks to both.  The images (with suitable reduction in storage size and
# consequent reduction in image quality as well) are available online at
# www.ibiblio.org/apollo.  If for some reason you find that the images are
# illegible, contact me at info@sandroid.org about getting access to the
# (much) higher-quality images which Paul actually created.
#
# Notations on the hardcopy document read, in part:
#
#    Assemble revision 001 of AGC program LMY99 by NASA 2021112-061
#    16:27 JULY 14, 1969

# Page 1159
# PROGRAM DESCRIPTION -- ENTIRE CONIC SUBROUTINE LOG SECTION DATE - 1 SEPTEMBER 1967
# MOD NO. - 0 LOG SECTION - CONIC SUBROUTINES
# MOD BY KRAUSE ASSEMBLY - COLOSSUS REVISION 88
#
# FUNCTIONAL DESCRIPTION -
# THE FOLLOWING SET OF SUBROUTINES SOLVE VARIOUS PROBLEMS INVOLVING THE TRAJECTORY PRODUCED BY A CENTRAL
# INVERSE-SQUARE FORCE ACTING ON A POINT MASS, AS OUTLINED IN THE CMC AND LGC LUNAR LANDING MISSION GSOP, SECTION
# 5.5.1.2.  A GENERAL USAGE POINT-OF-VIEW WAS TAKEN IN FORMULATING, MECHANIZING, AND SCALING THE SUBROUTINES,
# RATHER THAN OPTIMIZING EACH FOR A PARTICULAR USE.  THEREFORE, MULTIPLE USAGE CAN BE MADE OF THE SUBROUTINES
# INVOLVING ANY REALISTIC SET OF CONSTRAINTS.  IT SHOULD BE NOTED THAT ONLY ONE SET OF CODING IS USED, WHETHER THE
# EARTH, MOON, OR ANY OTHER CELESTIAL BODY IS SPECIFIED AS THE CENTRAL BODY OF THE PROBLEM, PROVIDED ONE OBSERVES
# THE INHERENT SCALE CHANGE REQUIRED IN POSITION, VELOCITY, MU, AND TIME, AS OUTLINES IN MISSION PROGRAMMING
# DEFINITION MEMO NO. 10.  THIS CAN BE ACCOMPLISHED BY SIMPLY ADDING TO THE MUTABLE AND INITIALIZING THE SUBROUTINES
# APPROPRIATELY.
#
# DUE TO THE UNIFORMITY OF THE EQUATIONS INVOLVED, CODING WAS MINIMIZED BY TREATING INDIVIDUAL EQUATIONS AND
# BLOCKS OF EQUATIONS AS SUBROUTINES OF LOWER RANK WHENEVER POSSIBLE.  AS A RESULT, THREE BY-PRODUCTS SUBROUTINES,
# DIRECTLY USABLE AS INDEPENDENT SUBROUTINES, WERE GENERATED.
#
# RESTRICTIONS -
# THE ONLY LIMITATION IN THE SCOPE OF THE PROBLEM WHICH CAN BE SOLVED BY A PARTICULAR SUBROUTINE IS THE SCALING
# LIMIT OF EACH PARAMETER AS SPECIFIED IN THE GSOP.  THESE SCALING LIMITS WERE CHOSEN SO THAT ALL FEASIBLE TRAJECTORIES
# COULD BE HANDLED.
#
# SINCE THE SUBROUTINES (EXCEPT KEPLER) USE COMMON SUBROUTINES OF LOWER RANK WHICH USE ERASABLE OTHER THAN
# THE PUSHLIST (DUE TO ITS LIMITED SIZE) AND COMMON INTERPRETIVE SWITCHES, THE CONIC SUBROUTINES CANNOT BE ALLOWED
# TO INTERRUPT EACH OTHER.  IT IS UP TO THE USER TO GUARANTEE THIS CONDITION.

# Page 1160
# PROGRAM DESCRIPTION - KEPLER SUBROUTINE DATE - 11 OCTOBER 1967
# MOD NO. -1 LOG SECTION - CONIC SUBROUTINES
# MOD BY KRAUSE ASSEMBLY - COLOSSUS 103 AND SUNDANCE 222
# MOD NO. - 2 (AUGUST 1968) BY ROBERTSON: TO PERMIT BACKDATING BY MORE THAN ONE ORBITAL PERIOD.
# MOD NO. - 3 (DEC 1968) BY ROBERTSON:  SUPPRESSION OF X-MODULO-ING
# MOD NO. - 4 (JAN 1969) BY ROBERTSON:  CLEAR OVFIND AT KEPLER ENTRY
#
# FUNCTIONAL DESCRIPTION -
# THIS SUBROUTINE, GIVEN AN INITIAL STATE VECTOR AND THE DESIRED TRANSFER TIME THROUGH WHICH THE STATE IS TO
# BE UPDATED ALONG A CONIC TRAJECTORY, COMPUTES THE NEW, UPDATED STATE VECTOR.  THE TRAJECTORY MAY BE ANY CONIC
# SECTION - CIRCULAR, ELLIPTIC, PARABOLIC, HYPERPOLIC, OR RECTILINEAR WITH RESPECT TO THE EARTH OR THE MOON.  THE
# USE OF THE SUBROUTINE CAN BE EXTENDED USING OTHER PRIMARY BODIES BY SIMPLE ADDITIONS TO THE MUTABLE WITHOUT
# INTRODUCING ANY CODING CHANGES, ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY.  AN ITERATION
# TECHNIQUE IS UTILIZED IN THE COMPUTATION.
#
# IF A NEGATIVE TIME-OF-FLIGHT IS INPUT, THE PROGRAM WILL SOLVE FOR THE STATE WHICH WOULD BE PRODUCED BY
# EXTRAPOLATING THE POSITION BACKWARD IN TIME.
#
# IF THE ABSOLUTE VALUE OF THE DESIRED TRANSFER TIME EXCEEDS THE ORBITAL PERIOD, THE SUBROUTINE, THROUGH A
#  MODULAR TECHNIQUE, WILL COMPUTE THE STATE CORRESPONDING TO THE DESIRED TIME (WHETHER POSITIVE OR NEGATIVE).
#
# THE RESTRICTIONS ARE -
# 1. (PREVIOUS RESTRICTION ON THE NEGATIVE DESIRED TRANSFER TIME IS NOW DELETED.)
# 2. THE PARAMETERS IN THE PROBLEM CANNOT EXCEED THEIR SCALING LIMITS AS SPECIFIED IN THE GSOP.  IF
#    ANY OF THESE LIMITS ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS.
#
#  THE NUMBER OF ITERATIONS AND, THEREFORE, THE COMPUTATION SPEED IS DEPENDENT ON THE ACCURACY OF THE
# GUESS, XKFPNEW.  THE AGC COMPUTATION TIME IS APPROXIMATELY .061 SECONDS FOR INITIALIZATION, .065 SECONDS FOR THE
# FINAL COMPUTATIONS, PLUS .083 SECONDS FOR EACH ITERATION.
#
# REFERENCES -
# R-479, MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP, SECTION 5.5, SGA
# MEMO 67-4.
#
# INPUT - ERASABLE INITIALIZATION REQUIRED
# * SCALE FACTOR  *
# VARIABLE *IN POWERS OF 2 *  DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# RRECT * +29 FOR EARTH * DP INITIAL POSITION VECTOR IN METERS
# * +27 FOR MOON *
# Page 1161
# VRECT * +7 FOR EARTH * DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND
# * +5 FOR MOON *
# X1 (38D) * NONE * INDEX REGISTER SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON,
# * * RESPECTIVELY, IS THE CENTRAL BODY
# TAU * +28 * DESIRED TRANSFER TIME IN CENTISECONDS (DP)
# * * MAY BE POS OR NEG AND ABSOLUTE VALUE MAY BE GREATER OR LESS THAN ONE ORBITAL PERIOD.
# XKEPNEW * +17 FOR EARTH * DP GUESS OF ROOT X OF KEPLERS EQN IN SQRT(METERS). SIGN SHOULD AGREE WITH THAT OF TAU.
# * +16 FOR MOON * AND ABS VALUE SHOULD BE LESS THAN THAT CORRESPONDING TO A PERIOD, VIZ, 2PI SQRT(SEMI-
# * * MAJOR AXIS), FOR SPEED OF CONVERGENCE, BUT IF EITHER CONDITION FAILS, XKEPNEW IS RESET
# * * BY KEPLER TO A POOR BUT VALID GUESS.
# TC * +28 * DP PREV. VALUE OF TIME IN CENTISECS.  MUST BE LESS THAN ONE ORBITAL PERIOD.
# XPREV * +17 FOR EARTH * DP PREV. VALUE OF X IN SQRT(METERS).  MUST BE LESS THAN AN X CORRESPONDING TO ONE
# * +16 FOR MOON * ORBITAL PERIOD, VIZ, 2PI SQRT(SEMI-MAJOR AXIS)
#
# SUBROUTINES CALLED -
# DELTIME
#
# CALLING SEQUENCE AND NORMAL EXIT MODES -
# KEPRTN-2 GOTO # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY.
# KEPRTN-1 KEPLER # RETURNS WITH XPREV IN MPAC. PL IS AT 0.
# KEPRTN ... # CONTINUE
#
# KEPLER MUST NOT BE CALLED DIRECTLY SINCE AN INTERRUPTION OF IT WOULD DESTROY THE ERASABLES IT NEEDS TO COMPLETE
# THE INTERRUPTED JOB.  THEREFORE THE USER MUST CALL CSMCONIC OR LEMCONIC WHICH GUARANTEES NO INTERRUPTS AND WHICH
# ALSO CALLS KEPPREP TO COMPUTE A GUESS OF XKEPNEW.
#
# ABORT EXIT MODES -
# NONE
#
# OUTPUT -
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# RCV * +29 FOR EARTH * DP TERMINAL POSITION VECTOR IN METERS
# * +27 FOR MOON *
# VCV * +7 FOR EARTH * DP TERMINAL VELOCITY VECTOR IN METERS/CENTISEC
# * +5 FOR MOON *
# TC * +28 * DP TRANSFER TIME IN CENTISECS TO WHICH KEPLER CONVERGED. ALWAYS LESS THAN ONE PERIOD.
# XPREV * +17 FOR EARTH * DP VALUE OF X IN SQRT(METERS) TO WHICH KEPLER CONVERGED. ALWAYS LESS THAN THE X
# * +16 FOR MOON * CORRESPONDING TO ONE PERIOD.
# Page 1162
# FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS.
#
# DEBRIS -
# PARAMETERS WHICH MAY BE OF USE -
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# URRECT * +1 * DP UNIT VECTOR OF INITIAL POSITION
# R1 * +29 FOR EARTH * DP MAGNITUDE OF INITIAL POSITION IN METERS
# * +27 FOR MOON *
#  ALPHA * -22 FOR EARTH * DP INVERSE OF SEMI-MAJOR AXIS IN 1/METERS
# * -20 FOR MOON *
# TMODULO * +28 * DP INTEGRAL NUMBER OF PERIODS IN CENTISECS, WHICH WAS SUBTRACTED FROM TAU, TO PRODUCE A
# * * TAU. OF LESS THAN ONE PERIOD.
#
# PARAMETERS OF NO USE -
# DP PARAMETERS - EPSILONT, DELX, DELT, RCNORM, XMODULO, PLUS PUSHLIST REGISTERS 0 THROUGH 39D.

# Page 1163

# PROGRAM DESCRIPTION - LAMBERT SUBROUTINE DATE - 1 SEPTEMBER 1967
# MOD NO. - 0 LOG SECTION - CONIC SUBROUTINES
# MOD BY KRAUSE ASSEMBLY - COLOSSUS REVISION 88
#
# FUNCTIONAL DESCRIPTION -
# THIS SUBROUTINE CALCULATES THE INITIAL VELOCITY REQUIRED TO TRANSFER A POINT-MASS ALONG A CONIC TRAJECTORY
# FROM AN INITIAL POSITION TO A TERMINAL POSITION IN A PRESCRIBED TIME INTERVAL.  THE RESULTING TRAJECTORY MAY BE
# A SECTION OF A CIRCLE, ELLIPSE, PARABOLA, OR HYPERBOLA WITH RESPECT TO TEH EARTH OR THE MOON.  THE USE OF THE
# SUBROUTINE CAN BE EXTEDED USING OTHER PRIMARY BODIES BY SIMPLE ADDITIONS TO THE MUTABLE WITHOUT INTRODUCING ANY
# CODING CHANGES, ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY.  AN ITERATION TECHNIQUE IS
# UTILIZED IN THE COMPUTATION.
#
# THE RESTRICTIONS ARE: -
# 1. RECTILINEAR TRAJECTORIES CANNOT BE COMPUTED.
# 2. AN ACCURACY DEGRADATION OCCURS AS THE COSINE OF THE TRUE ANOMALY DIFFERENCE APPROACHES +1.0.
# 3. THE ANGLE BETWEEN ANY POSITION VECTOR AND ITS VELOCITY VECTOR MUST BE GREATER THAN 1 DEGREE 47.5 MINUTES
#    AND LESS THAN 178 DEGREES 12.5 MINUTES.
# 4. NEGATIVE TRANSFER TIME IS AMBIGUOUS AND WILL RESULT IN NO SOLUTION.
# 5. THE PARAMETERS IN THE PROBLEM MUST NOT EXCEED THEIR SCALING LIMITS SPECIFIED IN THE GSOP.  IF THE
#    LIMITS ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS.
#
# THE NUMBER OF ITERATIONS AND, THEREFORE, THE COMPUTATIONS SPEED IS DEPENDENT ON THE ACCURACY OF THE FIRST
# GUESS OF THE INDEPENDENT VARIABLE, COGA.  THE AGC COMPUTATION TIME IS APPROXIMATELY
# .105 SECONDS FOR INITIALIZATION, .069 SECONDS FOR FINAL COMPUTATIONS, PLUS .205 SECONDS FOR EACH ITERATION.
#
# REFERENCES -
# R-479, MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP - SECTION 5.5, SGA MEMO 67-8,
# SGA MEMO 67-4.
#
# INPUT - ERASABLE INITIALIZATION REQUIRED
#
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# R1VEC * +29 FOR EARTH * DP INITIAL POSITION VECTOR IN METERS
# * +27 FOR MOON *
# R2VEC * +29 FOR EARTH * DP TARGET OR TERMINAL POSITION VECTOR IN METERS
# * +27 FOR MOON *
# TDESIRED * +28 * DP DESIRED TRANSFER TIME IN CENTISECONDS
# X1 (38D) * NONE * INDEX REGISTER SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON,
# * * RESPECTIVELY, IS THE CENTRAL BODY
# GEOMSGN * NONE * SP +.5 IF DESIRED TRANSFER ANGLE IS LESS THAN 180 DEGREES, -.5 IF GREATER THAN 180 DEG.
# GUESSW * NONE * AN INTERPRETER SWITCH TO BE SET IF NO GUESS OF COGA IS AVAILABLE, CLEAR IF A GUESS OF
# Page 1164
# * * COGA IS TO BE USED BY LAMBERT
# COGA * +5 * DP GUESS OF COTANGENT OF FLIGHT PATH ANGLE (MEASURED FROM VERTICAL).  THIS WILL BE
# * * IGNORED IF GUESSW IS SET.
# NORMSW * NONE * AN INTERPRETER SWITCH TO BE SET IF UN IS TO BE AN INPUT TO THE SUBROUTINE, CLEAR IF
# * * LAMBERT IS TO COMPUTE ITS OWN NORMAL (UN).
# UN * +1 * DP UNIT NORMAL TO THE DESIRED ORBIT PLANE IN THE DIRECTION OF THE RESULTING ANGULAR
# * * MOMENTUM VECTOR.  THIS WILL BE IGNORED IF NORMSW IS CLEAR.
# VTARGTAG * NONE * A S.P. TAG TO BE SET TO ZERO IF LAMBERT IS TO COMPUTE THE VELOCITY AT R2VEC AS WELL AS
# * * AT R1VEC.
# ITERCTR * NONE * A S.P. COUNTER WHICH SPECIFIES THE MAXIMUM NUMBER OF ITERATIONS ALLOWABLE.
# * * (AN ITERATION MEANS A PASS THRU KEPLER EQN (DELTIME). AT LEAST ONE OF THESE MUST
# * * ALWAYS OCCUR, EVEN IF COGA CORRESPONDING TO SOLUTION WERE INPUT AS A GUESS.)
# * * TWENTY ITERATIONS ARE SUFFICIENT TO SOLVE ALL PROBLEMS INCLUDING THOSE WITHOUT GUESS.
#
# SUBROUTINES CALLED -
# GEOM, GETX, DELTIME, ITERATOR, LAMENTER (PART OF NEWSTATE)
#
# CALLING SEQUENCE AND NORMAL EXIT MODES -
# L  CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY.
# L+1 LAMBERT # RETURNS WITH PL AT 0 AND WITH VVEC IN MPAC IF VTARGTAG WAS WAS NON-ZERO OR VTARGET
# # IN MPAC IF VTARGTAG WAS ZERO
# L+2 BON # CONTINUE IF SOLNSW CLEAR SINCE SOLUTION IS ACCEPTABLE
# L+3 SOLNSW
# L+4 LAMABORT
#
# IF A LAMBERT RESULT IS TO BE A FIRST GUESS FOR THE NEXT LAMBERT CALCULATION, COGA MUST BE PRESERVED AND
# GUESSW MUST BE CLEAR FOR EACH SUCCEEDING LAMBERT CALL.
#
# ABORT EXIT MODES -
# IF SOLNSW WAS SET UPON EXITING, EITHER LAMBERT WAS ASKED TO COMPUTE A TRANSFER TOO NEAR 0 OR 360 DEG, OR T
# WAS TOO SMALL TO PRODUCE A REALISTIC TRANSFER BETWEEN R1VEC AND R2VEC.  IN EITHER CASE THE FIX MUST BE MADE
# ACCORDING TO THE NEEDS OF THE PARTICULAR USER.  THE ABORT EXIT MODE MAY BE CODED AS ...
# LAMBABORT DLOAD ABS # A MEASURE OF THE PROXIMITY TO 0 OR
# 1-CSTH # 360 DEGREES.
# DSU BMN
# ONEBIT
# CHANGER2 # CHANGE R2VEC DIRECTION SLIGHTLY.
# DLOAD DAD
# TDESIRED
# SOMETIME
# STCALL TDESIRED # INCRESE TDESIRED
# LAMBERT
#
# Page 1165
# OUTPUT -
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# VVEC * +7 FOR EARTH * DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND REQUIRED TO SATISFY THE BOUNDARY VALUE
# * +5 FOR MOON * PROBLEM.
# VTARGET * +7 FOR EARTH * DP RESULTANT VELOCITY VECTOR AT R2VEC IN METERS/CENTISECOND.
# * +5 FOR MOON *
# SOLNSW * NONE * INTERPRETER SWITCH WHICH IS SET IF THE SUBROUTINE CANNOT SOLVE THE PROBLEM, CLEAR IF THE
# *  * SOLUTION EXISTS.
#
# FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS.
#
# DEBRIS -
# PARAMETERS WHICH MAY BE OF USE -
#
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# SNTH * +1 * DP SIN OF ANGLE BETWEEN R1VEC AND R2VEC
# CSTH * +1 * DP COSINE OF ANGLE
# 1-CSTH * +2 * DP 1-CSTH
# COGA * +5 * DP COTAN OF INITIAL REQUIRED FLIGHT PATH ANGLE MEASURED FROM VERTICAL
# P * +4 * DP RATIO OF SEMILATUS RECTUM TO INITIAL RADIUS
# R1A * +6 * DP RATIO OF INITIAL RADIUS TO SEMI-MAJOR AXIS
# R1 (32D) * +29 FOR EARTH * DP INITIAL RADIUS IN METERS
# * +27 FOR MOON *
# UR1 * +1 * DP UNIT VECTOR OF R1VEC
# U2 * +1 * DP UNIT VECTOR OF R2VEC
#
# PARAMETERS OF NO USE -
# DP PARAMETERS - EPSILONL, CSTH-RHO, TPREV, TERRLAMB, R2, RTNLAMB (SP), PLUS PUSHLIST REGISTER 0 THROUGH 41D
# ADDITIONAL INTERPRETIVE SWITCHES USED - INFINFLG, 360SW, SLOPESW, ORDERSW

# Page 1166
# PROGRAM DESCRIPTION - TIME-THETA SUBROUTINE DATE - 1 SEPTEMBER 1967
# MOD NO. - 0 LOG SECTION - CONIC SUBROUTINES
# MOD BY KRAUSE ASSEMBLY - COLOSSUS REVISION 88
#
# FUNCTIONAL DESCRIPTION -
# THIS SUBROUTINE, GIVEN AN INITIAL STATE VECTOR AND A DESIRED TRUE-ANOMALY-DIFFERENCE THROUGH WHICH THE
# STATE IS TO BE UPDATED ALONG A CONIC TRAJECTORY, CALCULATES THE CORRESPONDING TIME-OF-FLIGHT AND, IN ADDITION,
# PROVIDES THE OPTION OF COMPUTING THE NEW UPDATED STATE VECTOR.  THE RESULTING TRAJECTORY MAY BE A SECTION OF A
# CIRCLE, ELLIPSE, PARABOLA, OR HYPERBOLA WITH RESPECT TO THE EARTH OR THE MOON.  THE USE OF TEH SUBROUTINE CAN BE
# EXTENDED USING OTHER PRIMARY BODIES BY SIMPLE ADDITIONS TO THE MUTABLE WTIHOUT INTRODUCING ANY CODING CHANGES,
# ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY.
#
# THE RESTRICTIONS ARE -
# 1. THE ANGLE BETWEEN ANY POSITION VECTOR AND ITS VELOCITY VECTOR MUST BE GREATER THAN 1 DEGREE 47.5 MINUTES
#    AND LESS THAN 178 DEGREES 12.5 MINUTES.
# 2. THE PARAMETERS IN THE PROBLEM MUST NOT EXCEED THEIR SCALING LIMITS SPECIFIED IN THE GSOP.  IF THE LIMITS
#    ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLY.
#
# THE AGC COMPUTATION TIME IS APPROXIMATELY .292 SECONDS.
#
# REFERENCES -
# R-479, MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP-SECTION 5.5, SGA MEMO 67-8.
#
# INPUT - ERASABLE INITIALIZATION REQUIRED
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# RVEC * +29 FOR EARTH * DP INITIAL POSITION VECTOR IN METERS
# * +27 FOR MOON *
# VVEC * +7 FOR EARTH * DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND
# * +5 FOR MOON *
# SNTH * +1 * DP SINE OF TRUE-ANOMALY-DIFFERENCE THROUGH WHICH THE STATE IS TO BE UPDATED
# CSTH * +1 * DP COSINE OF THE ANGLE
# RVSW * NONE * AN INTERPRETIVE SWITCH TO BE SET IF ONLY TIME IS TO BE AN OUTPUT, CLEAR IF THE NEW STATE
# * * IS TO BE COMPUTED ALSO.
# X1 (38D) * NONE * INDEX REGISTER TO BE SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON,
# * * RESPECTIVELY, IS THE CENTRAL BODY.
#
# SUBROUTINES CALLED -
# Page 1167
# PARAM, GEOM, GETX, DELTIME, NEWSTATE
#
# CALLING SEQUENCE AND NORMAL EXIT MODES -
# IF ONLY TIME IS DESIRED AS OUTPUT -
# L SET CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY.
# L+1 RVSW
# L+2 TIMETHET # RETURN WITH PL AT 0 AND T IN MPAC
# L+3 ... # CONTINUE
#
# IF THE UPDATE STATE VECTOR IS DESIRED AS WELL -
# L CLEAR CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY.
# L+1 RVSW
# L+2 TIMETHET # RETURNS WITH PL AT 6.  THE INITIAL POSITION VECTOR IS IN 0D OF THE PUSHLIST AND
# # THE INITIAL VELOCITY VECTOR IN MPAC.
# L+3 STOVL NEWVVEC
# L+4 STADR
# L+5 STORE NEWRVEC # NEWVVEC AND NEWRVEC ARE SYMBOLIC REPRESENTATIONS OF THE USERS LOCATIONS.
# L+6 ... # CONTINUE.
#
# ABORT EXIT MODES -
# IF COGAFLAG AND/OR INFINFLG IS SET AT THE EXIT TO TIME-THETA, TIME-THETA WILL TRANSFER TO POODOO WITH
# AN ALARM CODE (ORIGINALLY 00607), AND NOT RETURN TO THE CALLING PROGRAM.  (PCR 692 AND 721).
#
# OUTPUT -
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# T (30D) * +28 * DP TRANSFER TIME IN CENTISECONDS
# INFINFLG * NONE * AN INTERPRETIVE SWITCH WHICH IS SET IF THE TRANSFER ANGLE REQUIRES CLOSURE THROUGH
# * * INFINITY (NO SOLUTION), CLEAR IF A PHYSICAL SOLUTION IS POSSIBLE.
# COGAFLAG * NONE * AN INTERPRETIVE SWITCH WHICH IS SET IF RESTRICTION 1 HAS BEEN VIOLATED (NO SOLUTION),
# *  * CLEAR IF NOT.
#
# IN ADDITION, IF RVSW IS CLEAR, THE FOLLOWING ARE OUTPUT -
# MPAC -  * +7 FOR EARTH * DP TERMINAL VELOCITY VECTOR IN METERS/CENTISEC.
#  MPAC +5 * +5 FOR MOON *
# 0D - 5D * +29 FOR EARTH * DP TERMINAL POSITION VECTOR IN METERS (PL AT 6D)
# * +27 FOR MOON *
#
# FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS.
#
# Page 1168
# DEBRIS -
# PARAMETERS WHICH MAY BE OF USE -
#
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# R1 (32D) * +29 FOR EARTH * DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS
# * +27 FOR MOON *
# R1A * +6 * DP RATIO OF R1 TO SEMIMAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES)
# P * +4 * DP RATIO OF SEMILATUS RECTUM TO R1
# COGA * +5 * DP COTAN OF ANGLE BETWEEN RVEC AND VVEC
# UR1 * +1 * DP UNIT VECTOR OF RVEC
# U2 * +1 * DP UNIT VECTOR OF VVEC
# UN * +1 * DP UNIT VECTOR OF UR1*U2
#
# PARAMETERS OF NO USE -
# SP PARAMETERS -- RTNTT, GEOMSGN, RTNPRM, MAGVEC2=R2 (DP), PLUS PUSHLIST LOCATIONS 0-11D, 14D-21D, 24D-39D, 41D
# ADDITIONAL INTERPRETIVE SWITCHES USED -- NORMSW, 360SW

# Page 1169
# PROGRAM DESCRIPTION - TIME-RADIUS SUBROUTINE DATE - 11 OCTOBER 1967
# MOD NO. -1 LOG SECTION - CONIC SUBROUTINES
# MOD BY KRAUSE ASSEMBLY - COLOSSUS REVISION 88
#
# FUNCTIONAL DESCRIPTION -
#
# THIS SUBROUTINE, GIVEN AN INITIAL STATE VECTOR AND A DESIRED RADIUS TO WHICH THE
# STATE IS TO BE UPDATED ALONG A CONIC TRAJECTORY, CALCULATES THE CORRESPONDING TIME-OF-FLIGHT AND, IN ADDITION,
# PROVIDES THE OPTION OF COMPUTING THE NEW UPDATED STATE VECTOR.  THE RESULTING TRAJECTORY MAY BE A SECTION OF A
# CIRCLE, ELLIPSE, PARABOLA, OR HYPERBOLA WITH RESPECT TO THE EARTH OR THE MOON.  THE USE OF THE SUBROUTINE CAN BE
# EXTENDED USING OTHER PRIMARY BODIES BY SIMMPE ADDITIONS TO THE MUTABLE WITHOUT INTRODUCING ANY CODING CHANGES,
# ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY.
#
# IF THE DESIRED RADIUS IS BEYOND THE RADIUS OF APOCENTER OF THE CONIC OR BELOW THE RADIUS OF PERICENTER,
# APSESW WILL BE SET AND THE SUBROUTINE WILL RETURN THE APOCENTER OR PERICENTER SOLUTION, RESPECTIVELY.
#
# THE RESTRICTIONS ARE -
# 1. THE ANGLE BETWEEN ANY POSITION VECTOR AND ITS VELOCITY VECTOR MUST BE GREATER THAN 1 DEGREE 47.5 MINUTES
#    AND LESS THAN 178 DEGREES 12.5 MINUTES.
# 2. THE PARAMETERS IN THE PROBLEM MUST NOT EXCEED THEIR SCALING LIMITS SPECIFIED IN THE GSOP.  IF THE LIMITS
#    ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS.
# 3. AN ACCURACY DEGRADATION OCCURS AS THE SENSITIVITIES OF TIME AND UPDATED STATE VECTOR TO CHANGES IN
#    RDESIRED INCREASE.  THIS WILL OCCUR NEAR EITHER APSIS OF THE CONIC AND WHEN THE CONIC IS NEARLY CIRCULAR.  IN
#    PARTICULAR, IF THE CONIC IS AN EXACT CIRCLE, THE PROBLEM IS UNDEFINED AND THE SUBROUTINE WILL ABORT.
#
# THE AGC COMPUTATION TIME IS APPROXIMATELY .363 SECONDS.
#
# REFERENCES -
# R-479, MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP-SECTION 5.5, SGA MEMO 67-8.
#
# INPUT - ERASABLE INITIALIZATION REQUIRED.
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# RVEC * +29 FOR EARTH * DP INITIAL POSITION VECTOR IN METERS
# * +27 FOR MOON *
# VVEC * +7 FOR EARTH * DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND
# * +5 FOR MOON *
# RDESIRED * +29 FOR EARTH * DP TERMINAL RADIAL DISTANCE ON CONIC TRAJECTORY FOR WHICH TRANSFER TIME IS TO BE
# * +27 FOR MOON * COMPUTED
# SGNRDOT * NONE * SP TAG SET TO +.5 OR -.5 ACCORDING TO WHETHER THE RADIAL VELOCITY AT RDESIRED IS TO BE
# * * POSITIVE OR NEGATIVE, RESPECTIVELY.   THIS TAG REDUCES THE DOUBLE-VALUED PROBLEM TO A
# Page 1170
# * * SINGLE-VALUED PROBLEM.
# X1 (38D) * NONE * INDEX REGISTER TO BE SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON,
# *  * RESPECTIVELY, IS THE CENTRAL BODY.
# RVSW * NONE * AN INTERPRETIVE SWITCH TO BE SET IF ONLY TIME IS TO BE AN OUTPUT, CLEAR IF THE NEW STATE
# * * IS TO BE COMPUTED ALSO.
#
# SUBROUTINES CALLED -
# PARAM, GEOM, GETX, DELTIME, NEWSTATE
#
# CALLING SEQUENCE AND NORMAL EXIT MODES -
#
# IF ONLY TIME IS DESIRED AS OUTPUT -
# L SET CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY.
# L+1 RVSW
# L+2 TIMERAD # RETURN WITH PL AT 0 AND T IN MPAC
# L+3 ... # CONTINUE
#
# IF THE UPDATE STATE VECTOR IS DESIRED AS WELL -
# L CLEAR CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY.
# L+1 RVSW
# L+2 TIMERAD # RETURNS WITH PL AT 6.  THE INITIAL POSITION VECTOR IS IN 0D OF THE PUSHLIST AND
# # THE INITIAL VELOCITY VECTOR IN MPAC.
# L+3 STOVL NEWVVEC
# L+4 STADR
# L+5 STORE NEWRVEC # NEWVVEC AND NEWRVEC ARE SYMBOLIC REPRESENTATIONS OF THE USERS LOCATIONS.
# ... # CONTINUE
#
# ABORT EXIT MODES -
# IF SOLNSW AND/OR COGAFLAG AND/OR INFINFLG IS SET AT THE EXIT TO TIME-RADIUS, TIME-RADIUS WILL TRANSFER
# TO POODOO WITH AN ALARM CODE (ORIGINALLY 00607), AND NOT RETURN TO THE CALLING PROGRAM. (PCR 692 & 721)
#
# OUTPUT -
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# T (30D) * +28 * DP TRANSFER TIME IN CENTISECONDS.
# INFINFLG * NONE * AN INTERPRETIVE SWITCH WHICH IS SET IF RDESIRED AND SGNRDOT REQUIRE CLOSURE THROUGH
# * * INFINITY (NO SOLUTION), CLEAR IF A PHYSICAL SOLUTION IS POSSIBLE.
# COGAFLAG * NONE * AN INTERPRETIVE SWITCH WHICH IS SET IF RESTRICTION 1 HAS BEEN VIOLATED (NO SOLUTION),
# * * CLEAR IF NOT.
# APSESW * NONE * AN INTERPRETIVE SWITCH WHICH IS SET IF RDESIRED WAS GREATER THAN RADIUS OF APOCENTER OR
# Page 1171
# * * LESS THAN RADIUS OF PERICENTER.  THE APOCENTER OR PERICENTER SOLUTION, RESPECTIVELY,
# * * WILL THEN BE RETURNED.  THE SWITCH IS CLEAR IF RDESIRED WAS BETWEEN PERICENTER AND
# * * APOCENTER.
# SOLNSW * NONE * AN INTERPRETIVE SWITCH WHICH IS SET IF THE CONIC IS SO CLOSE TO A CIRCLE THAT THE TERMIN
# * * POINT IS AMBIGUOUS, VIOLATING RESTRICTION 3.  IF ECCENTRICITY IS GREATER THAN 2-TO-THE-
# * * MINUS-18, THE SWITCH IS CLEAR.
#
# IN ADDITION, IF RVSW IS CLEAR, THE FOLLOWING ARE OUTPUT -
#
# MPAC -  * +7 FOR EARTH * DP TERMINAL VELOCITY VECTOR IN METERS/CENTISEC.
#  MPAC +5 * +5 FOR MOON *
# 0D - 5D * +29 FOR EARTH * DP TERMINAL POSITION VECTOR IN METERS (PL AT 6D)
# * +27 FOR MOON *
#
# FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS.
#
# DEBRIS -
# PARAMETERS WHICH MAY BE OF USE -
#
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# R1 (32D) * +29 FOR EARTH * DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS
# * +27 FOR MOON *
# R1A * +6 * DP RATIO OF R1 TO SEMIMAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES)
# P * +4 * DP RATIO OF SEMILATUS RECTUM TO R1
# COGA * +5 * DP COTAN OF ANGLE BETWEEN RVEC AND VVEC
# UR1 * +1 * DP UNIT VECTOR OF RVEC
# U2 * +1 * DP UNIT VECTOR OF VVEC
# UN * +1 * DP UNIT VECTOR OF UR1*U2
# CSTH * +1 * DP COSINE OF TRUE ANOMALY DIFFERENCE BETWEEN RVEC AND RDESIRED.
# SNTH * +1 * DP SINE OF TRUE ANOMALY DIFFERENCE.
#
#  PARAMETERS OF NO USE -
# SP PARAMETERS -- RTNTT, GEOMSGN, RTNPRM, MAGVEC2=R2 (DP), PLUS PUSHLIST LOCATIONS 0-11D, 14D-21D, 24D-39D, 41D
# ADDITIONAL INTERPRETIVE SWITCHES USED -- NORMSW, 360SW
#

# Page 1172
# PROGRAM DESCRIPTION - APSIDES SUBROUTINE DATE - 1 SEPTEMBER 1967
# MOD NO. - 0 LOG SECTION - CONIC SUBROUTINES
# MOD BY KRAUSE ASSEMBLY - COLOSSUS REVISION 88
#
# FUNCTIONAL DESCRIPTION -
#
# THIS SUBROUTINE, GIVEN AN INITIAL STATE VECTOR CALCULATES THE RADIUS OF PERICENTER AND OF APOCENTER AND THE
# ECCENTRICITY OF THE RESULTING CONIC TRAJECTORY, WHICH MAY BE A STRAIGHT LINE,
# CIRCLE, ELLIPSE, PARABOLA, OR HYPERBOLA WITH RESPECT TO THE EARTH OR THE MOON.  THE USE OF THE SUBROUTINE CAN
# BE EXTENDED USING OTHER PRIMARY BODIES BY SIMPLE ADDITIONS TO THE MUTABLE WITHOUT INTRODUCING ANY CODING CHANGES,
# ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY.
#
# THE RESTRICTIONS ARE -
# 1. IF APOCENTER IS BEYOND THE SCALING OF POSITION, THE SCALE FACTOR LIMIT (536,870,910 METERS WITH RESPECT
#    TO THE EARTH OR 134,217,727.5 METERS WITH RESPECT TO THE MOON) WILL BE RETURNED.
# 2. THE PARAMETERS IN THE PROBLEM MUST NOT EXCEED THEIR SCALING LIMITS SPECIFIED IN THE GSOP.  IF THE LIMITS
#    ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS.
#
# THE AGC COMPUTATION TIME IS APPROXIMATELY .103 SECONDS.
#
# REFERENCES -
# MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP-SECTION 5.5
#
# INPUT - ERASABLE INITIALIZATION REQUIRED
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# RVEC * +29 FOR EARTH * DP INITIAL POSITION VECTOR IN METERS
# * +27 FOR MOON *
# VVEC * +7 FOR EARTH * DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND
# * +5 FOR MOON *
# X1 (38D) * NONE * INDEX REGISTER TO BE SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON,
# * * RESPECTIVELY, IS THE CENTRAL BODY.
#
# SUBROUTINES CALLED -
# PARAM, GEOM
#
# CALLING SEQUENCE AND NORMAL EXIT MODES -
# Page 1173
# IF ONLY TIME IS DESIRED AS OUTPUT -
# L CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY.
# L+1 APSIDES # RETURNS WITH PL AT 0, RADIUS OF APOCENTER IN MPAC AND RADIUS OF PERICENTER IN 0D
# L+2 STODL APOAPSE
# L+3 0D
# L+4 STORE PERIAPSE # APOAPSE AND PERIAPSE ARE SYMBOLIC REPRESENTATIONS OF THE USERS LOCATIONS
# L+5 ... # CONTINUE
#
# OUTPUT -
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# MPAC * +29 FOR EARTH * DP RADIUS OF APOCENTER IN METERS
# * +27 FOR MOON *
# 0D-1D * +29 FOR EARTH * DP RADIUS OF PERICENTER IN METERS
# * +27 FOR MOON *
# ECC * +3 * DP ECCENTRICITY OF CONIC TRAJECTORY.
#
# FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS.
#
# DEBRIS -
#
# PARAMETERS WHICH MAY BE OF USE -
#
# * SCALE FACTOR *
# VARIABLE *IN POWERS OF 2 * DESCRIPTION AND REMARKS
# -------- *-------------- * -----------------------
# R1 (32D) * +29 FOR EARTH * DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS
# * +27 FOR MOON *
# R1A * +6 * DP RATIO OF R1 TO SEMI-MAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES)
# P * +4 * DP RATIO OF SEMILATUS RECTUM TO R1
# COGA * +5 * DP COTAN OF ANGLE BETWEEN RVEC AND VVEC
# UR1 * +1 * DP UNIT VECTOR OF RVEC
# U2 * +1 * DP UNIT VECTOR OF VVEC
# UN * +1 * DP UNIT VECTOR OF UR1*U2
# MAGVEC2 * +7 FOR EARTH * DP MAGNITUDE OF VVEC
# * +5 FOR MOON *
#
# PARAMETERS OF NO USE -
# SP PARAMETERS - RTNAPSE, GEOMSGN, RTNPRM, PLUS PUSHLIST LOCATIONS 0-5, 10D-11D, 14D-21D, 31D-38D.
# ADDITIONAL INTERPRETIVE SWITCHES USED - NORMSW

  SETLOC CONICS
# Page 1174
  BANK

  COUNT* $$/CONIC
  EBANK= UR1
KEPLERN  SETPD BOV
   0
   +1
  VLOAD*
   MUTABLE,1
  STOVL 14D
   RRECT
  UNIT SSP
   ITERCTR
   20D
  STODL URRECT
   36D
  STOVL R1
   RRECT
  DOT SL1R
   VRECT
  DMP SL1R
   1/ROOTMU # 1/ROOTMU (-17 OR -14)
  STOVL KEPC1  # C1=R.V/ROOTMU (+17 OR +16)
   VRECT
  VSQ DMPR
   1/MU  # 1/MU (-34 OR -28)
  DMP SL3
   R1
  DSU ROUND
   D1/64
  STORE KEPC2  # C2=RV.V/MU -1 (+6)

  BDSU SR1R
   D1/64
  DDV
   R1
  STORE ALPHA  # ALPHA=(1-C2)/R1 (-22 OR -20)

  BPL DLOAD  # MAXIMUM X DEPENDS ON TYPE OF CONIC
   1REV
   -50SC  # -50SC (+12)
  DDV BOV
   ALPHA
   STOREMAX
  SQRT GOTO
   STOREMAX

1REV  SQRT BDDV
# Page 1175
   2PISC  # 2PISC (+6)
  BOV
   STOREMAX
STOREMAX STORE XMAX
  DMP PDDL
   1/ROOTMU
   ALPHA
  NORM PDDL
   X1
  SL* DDV
   0 -6,1
  BOV BMN
   MODDONE
   MODDONE  # MPAC=PERIOD
PERIODCH PDDL ABS  # 0D=PERIOD
   TAU.
  DSU BMN
   0D
   MODDONE
  SIGN
   TAU.
  STODL TAU.
  GOTO
   PERIODCH
MODDONE  SETPD DLOAD
   0
   XKEPNEW
  STORE X
  SIGN BZE
   TAU.
   BADX
  BMN ABS
   BADX
  DSU BPL
   XMAX
   BADX
STORBNDS DLOAD BPL
   TAU.
   STOREMIN
  DLOAD DCOMP
   XMAX
  STODL XMIN
   KEPZERO
  STORE XMAX
  GOTO
   DXCOMP
STOREMIN DLOAD
   KEPZERO
  STORE XMIN
DXCOMP  DLOAD DMPR
# Page 1176
   TAU.
   BEE22
  ABS
  STODL EPSILONT
   XPREV
XDIFF  BDSU
   X
  STORE DELX

KEPLOOP  DLOAD DSQ
   X  # X=XKEP
  NORM PUSH  # 0D=XSQ (+34 OR +32 -N1) PL AT 2
   X1
  DMP SRR*
   ALPHA
   0 -6,1
  STCALL XI  # XI=ALPHA XSQ (+6)
   DELTIME
  BOV BDSU
   TIMEOVFL # UNLIKELY
   TAU.
  STORE DELT  # DELT=DELINDEP
  ABS BDSU
   EPSILONT
  BPL DLOAD
   KEPCONVG
   T
  DSU NORM
   TC
   X1
  PDDL NORM
   DELX
   X2
  XSU,1 DMP
   X2
   DELT
  SLR* DDV
   1,1
  SR1 PUSH  # 0D=TRIAL DELX PL AT 2
  BPL DLOAD
   POSDELX
   X
  STORE XMAX  # MOVE MAX BOUND IN
  BDSU DSU  # PL AT 0
   XMIN
  BOV BPL
   NDXCHNGE
   NDXCHNGE
  DLOAD GOTO
# Page 1177
   0D
   NEWDELX

NDXCHNGE DLOAD DSU
   XMIN
   X
  DMPR GOTO  # TO FORCE MPAC +2 TO ZERO
   DP9/10
   NEWDELX

POSDELX  DLOAD
   X
  STORE XMIN  # MOVE MIN BOUND IN
  BDSU DSU  # PL AT 0
   XMAX
  BOV BMN
   PDXCHNGE
   PDXCHNGE
  DLOAD
   0D
NEWDELX  STORE DELX
  BZE DAD
   KEPCONVG
   X
  STODL X
   T
  STORE TC
BRNCHCTR RTB BHIZ
   CHECKCTR
   KEPCONVG
  GOTO
   KEPLOOP  # ITERATE

PDXCHNGE DLOAD DSU
   XMAX
   X
  DMPR GOTO  # TO FORCE MPAC +2 TO ZERO
   DP9/10
   NEWDELX

BADX  DLOAD SR1
   XMAX
  SIGN
   TAU.
  STORE X
  GOTO
# Page 1178
   STORBNDS
TIMEOVFL DLOAD BMN  # X WAS TOO BIG
   X
   NEGTOVFL
  STORE XMAX
CMNTOVFL DLOAD SR1
   DELX
  STORE DELX
  BZE BDSU
   KEPRTN
   X
  STODL X
   TC
  STORE T
  GOTO
   BRNCHCTR
NEGTOVFL STORE XMIN
  GOTO
   CMNTOVFL
KEPCONVG DLOAD SR4R
   R1
  DSU VXSC
   XSQC(XI)
   URRECT
  VSL1 PDDL  # 0D=(R1-XSQC(XI))URRECT (+33 OR +31)
   X
  DSQ NORM
   X1
  DMPR DMPR
   1/ROOTMU
   X
  DMP SRR*
   S(XI)
   0 -7,1
  BDSU
   T
  SL1 VXSC
   VRECT
  VSL1 VAD  #  PL AT 0
  VSL4
  STORE RCV  # RCV (+29 OR +27)

  ABVAL NORM
   X2
  STODL RCNORM
   XI
  DMPR DSU
   S(XI)
   D1/128
# Page 1179
  DMP SL1R
   ROOTMU
  DMP SLR*
   X
   0 -3,2
  DDV VXSC
   RCNORM
   URRECT
  VSL1 PDDL  # 0D=URRECT(XI S(XI)-1)X ROOTMU/RCV (+15
   XSQC(XI) # OR +13) PL AT 6
  SLR* DDV
   0 -4,2
   RCNORM
  BDSU VXSC
   D1/256
   VRECT
  VAD VSL8
  STADR   # PL AT 0
  STODL VCV  # VCV (+7 OR +5)
   T
  STODL TC
   X
  STORE XPREV
  GOTO
   KEPRTN

# Page 1180
DELTIME  EXIT   # MPAC=XI (+6), 0D=XSQ (+34 OR +32 -N1)
  TC POLY
  DEC 8
  2DEC .083333334
  2DEC -.266666684
  2DEC .406349155
  2DEC -.361198675
  2DEC .210153242
  2DEC -.086221951
  2DEC .026268812
  2DEC -.006163316
  2DEC .001177342
  2DEC -.000199055

  TC INTPRET
  STODL S(XI)
   XI
  EXIT
  TC POLY
  DEC 8
  2DEC .031250001
  2DEC -.166666719
  2DEC .355555413
  2DEC -.406347410
  2DEC .288962094
  2DEC -.140117894
  2DEC .049247387
  2DEC -.013081923
  2DEC .002806389
  2DEC -.000529414
  TC INTPRET
# Page 1181
  DMP SRR*  # PL AT 0
   0D
   0 -5,1
  STORE XSQC(XI) # XSQC(XI) (+33 OR +31)
  DMP SL1
   KEPC1
  RTB PDDL  # XCH WITH PL. 0D=C1 XSQ C(XI) (+49 OR +46)
   TPMODE  # PL AT 0,3
  DMP SRR*
   S(XI)
   0 -5,1
  DMP SL1
   KEPC2
  RTB PDDL  # 3D=C2 XSQ S(XI) (+35 OR +33)  PL AT 6
   TPMODE
   R1
  SR TAD  # PL AT 3
   6
  NORM DMP  # TO PRESERVE SIGNIF.
   X1
   X
  SR* TAD  # X(C2 XSQ S(XI) +R1) (+49 OR +46)  PL AT 0
   0 -3,1
  SL4R DMPR
   1/ROOTMU
  STORE T
  RVQ

# Page 1182
ITERATOR BONCLR DLOAD
   SLOPESW
   FIRSTIME
   DEP
  DSU NORM
   DEPREV
   X1
  PDDL NORM
   DELINDEP
   X2
  XSU,1 DMP
   X2
   DELDEP
  SLR* DDV  # PL UP 2
   1,1
  SR1 BOFF
   ORDERSW
   SGNCHECK
  ABS SIGN  # IN CASE 2ND DERIV. CHANGED SIGN, MUST
   DELDEP  # DISREGARD IT TO FIND MIN.

SGNCHECK PUSH BPL  # TRIAL DELINDEP PL DOWN 2
   POSDEL
  DLOAD BON
   INDEP
   ORDERSW
   MINCHECK
  STORE MAX  # IF NOT 2ND ORDER, CAN MOVE MAX BOUND IN.

MINCHECK BDSU DSU
   MIN
  BOV BPL
   MODNGDEL
   MODNGDEL
  GOTO
   DELOK

MODNGDEL DLOAD DSU  # TRIAL DELINDEP WOULD EXCEED MIN BOUND
   MIN
   INDEP
  DMP GOTO
   DP9/10
   NEWDEL

FIRSTIME DLOAD DMP
   MIN
   TWEEKIT  # DLOAD TWEEKIT(40D) SENSITIVE TO CHANGE.
  PDDL DMP  # S2(41D) SHOULDNT CONTAIN HI ORDER ONES
# Page 1183
   MAX
   TWEEKIT
  DSU
  SIGN GOTO
   DELDEP
   SGNCHECK

POSDEL  DLOAD BON
   INDEP
   ORDERSW
   MAXCHECK
  STORE MIN  # IF NOT 2ND ORDER, CAN MOVE MIN BOUND IN.

MAXCHECK BDSU DSU
   MAX
  BOV BMN
   MODPSDEL
   MODPSDEL
DELOK  DLOAD
   0D
NEWDEL  STORE DELINDEP
  RVQ

MODPSDEL DLOAD DSU
   MAX
   INDEP
  DMP GOTO
   DP9/10
   NEWDEL

CHECKCTR CS ONE
  INDEX FIXLOC
  AD ITERCTR
  INDEX FIXLOC
  TS ITERCTR
  TS MPAC
  TC DANZIG

# Page 1184
NEWSTATE DLOAD SR4R
   R1
  DSU VXSC
   XSQC(XI)
   UR1
  VSL1 PDDL  # 0D=(R1-XSQC(XI))UR1 (+33 OR 31) PL AT 6
   X
  DSQ NORM
   X1
  DMPR DMPR
   1/ROOTMU
   X
  DMP SRR*
   S(XI)
   0 -7,1
  BDSU
   T
  SL1 VXSC
   VVEC
  VSL1 VAD  # PL AT 0
  VSL4 PUSH
  ABVAL
LAMENTER NORM
   X1
  STODL R2
   XI
  DMP DSU
   S(XI)
   D1/128
  DMP SL1R
   ROOTMU
  DMP SLR*
   X
   0 -3,1
  DDV VXSC
   R2
   UR1
  VSL1 PDDL  # 6D=V2VEC PART (+15 OR 13) PL AT 12
   XSQC(XI)
  SLR* DDV
   0 -4,1
   R2
  BDSU
   D1/256
  VXSC VAD  # PL AT 6
   VVEC
  VSL8 RVQ

# Page 1185
  SETLOC CONICS1
  BANK

  COUNT* $$/CONIC
# DO NOT DISTURB THE ORDER OF THESE CDS, OVERLAYS HAVE BEEN MADE.
BEE17  DEC 0  # KEEP WITH D1/8 2DEC 1.0B-17 (0000004000)
D1/8  2DEC 1.0 B-3
D1/128  2DEC 1.0 B-7
D1/64  2DEC 1.0 B-6
D1/4  2DEC 1.0 B-2
D1/16  2DEC 1.0 B-4
D1/32  2DEC 1.0 B-5
D1/1024  2DEC 1.0 B-10
D1/256  2DEC 1.0 B-8
DP9/10  2DEC .9
KEPZERO  EQUALS LO6ZEROS
-50SC  2DEC -50.0 B-12
2PISC  2DEC 6.28318530 B-6
BEE19  EQUALS D1/32 -1 # 2DEC 1.0 B-19 (00000 01000)
BEE22  EQUALS D1/256 -1 # 2DEC 1.0 B-22 (00000 00100)
ONEBIT  2DEC 1.0 B-28
COGUPLIM 2DEC .999511597
COGLOLIM 2DEC -.999511597

# Page 1186
  SETLOC CONICS
  BANK

  COUNT* $$/CONIC
TIMETHET STQ SETPD  # PL AT 0
   RTNTT
   0
  BOV
   +1
  VLOAD PDVL  # SETUP FOR PARAM CALL PL AT 6
   RVEC
   VVEC
  CALL
   PARAM
  BOV CALL  # PL AT 0
   COGAOVFL
   GETX
COMMNOUT DLOAD BON
   XI
   INFINFLG
   ABTCONIC
  CLEAR CALL
   COGAFLAG
   DELTIME
  BON CALL
   RVSW
   RTNTT
   NEWSTATE
  GOTO
   RTNTT

COGAOVFL SETGO
   COGAFLAG
   ABTCONIC

  BANK 4
  SETLOC CONICS1
  BANK
  COUNT* $$/CONIC
PARAM  STQ CLEAR  # MPAC=V1VEC, 0D=R1VEC PL AT 6
   RTNPRM
   NORMSW
  CLEAR
   COGAFLAG
  SSP CALL
   GEOMSGN
   37777  # GAMMA ALWAYS LESS THAN 180DEG
   GEOM  # MPAC=SNGA (+1), 0D=CSGA (+1) PL AT 2
  STODL 36D  # 36D=SIN GAMMA (+1) PL AT 0
  SR DDV
# Page 1187
   5
   36D
  STOVL* COGA
   MUTABLE,1
  STODL 1/MU
   MAGVEC2
  DSQ NORM
   X1
  DMPR DMP
   1/MU
   R1
  SRR*
   0 -3,1
  PUSH BDSU  # 0D=R1 V1SQ/MU (+6) PL AT 2
   D1/32
  STODL R1A  # R1A (+6) PL AT 0

  DMP NORM
   36D
   X1
  DMP SR*
   36D
   0 -4,1
  STORE P  # P (+4)
  GOTO
   RTNPRM

# Page 1188
GEOM  UNIT   # MPAC=V2VEC, 0D=R1VEC PL AT 6
  STODL U2  # U2 (+1)
   36D
  STOVL MAGVEC2  # PL AT 0
  UNIT
  STORE UR1  # UR1 (+1)
  DOT SL1
   U2
  PDDL   # 0D=CSTH (+1) PL AT 2
   36D
  STOVL R1  # R1 (+29 OR +27)
   UR1
  VXV VSL1
   U2
  BON SIGN
   NORMSW
   HAVENORM
   GEOMSGN
  UNIT BOV
   COLINEAR
UNITNORM STODL UN  # UN (+1)
   36D
  SIGN RVQ  # MPAC=SNTH (+1), 34D=SNTH.SNTH (+2)
   GEOMSGN

COLINEAR VSR1 GOTO
   UNITNORM

HAVENORM ABVAL SIGN
   GEOMSGN
  RVQ   # MPAC=SNTH (+1), 34D=SNTH.SNTH (+2)

# Page 1189
  BANK 12
  SETLOC CONICS
  BANK

  COUNT* $$/CONIC
GETX  AXT,2 SSP  # ASSUMES P (+4) IN MPAC
   3
   S2
   1
  CLEAR
   360SW
  SQRT PDDL  # 0D=SQRT(P) PL AT 2
   CSTH
  SR1 BDSU
   D1/4
  PDDL SRR  # PL AT 4D
   SNTH
   6
  DDV   # PL AT 2
  BOV
   360CHECK
  DSU DMP
   COGA  # PL AT 0
  SL2R BOV
   360CHECK
WLOOP  PUSH DSQ  # 0D=W (+5) PL AT 2
  TLOAD PDDL  # 2D=WSQ (+10) PL AT 5
   MPAC
   R1A
  SR4 TAD  # PL AT 2
  BMN SQRT
   INFINITY
  ROUND DAD  # PL AT 0D
  BOV TIX,2
   RESETX2
   WLOOP

  BDDV BOV
   D1/128
   INFINITY
POLYCOEF BMN PUSH  # 0D=1/W (+2) OR 16/W (+6) PL AT 2
   INFINITY
  DSQ
  NORM DMP
   X1
   R1A
  SRR* EXIT
   0 -10D,1
  TC POLY
# Page 1190
  DEC 5
  2DEC .5
  2DEC -.166666770
  2DEC .100000392
  2DEC -.071401086
  2DEC .055503292
  2DEC -.047264098
  2DEC .040694204

  TC INTPRET
  DMP SL1R  # PL AT 0D
  PUSH BON
   360SW
   TRUE360X
XCOMMON  DSQ NORM
   X1
  DMP SRR*
   R1A
   0 -12D,1
  STODL XI  # XI (+6)
   R1
  SR1 SQRT
  ROUND DMP
  SL4R   # PL AT 0
  STORE X  # X (+17 OR +16)

  DSQ NORM
   X1
  PDDL DMP  # 0D=XSQ (+34 OR +32 -N1) PL AT 2
   P
   R1
  SL3 SQRT
  DMP SL3R
   COGA
  STODL KEPC1
   R1A
  BDSU CLEAR
   D1/64
   INFINFLG
  STORE KEPC2
  RVQ

# Page 1191
RESETX2  AXT,2
   3

360CHECK SETPD BPL
   0D
   INVRSEQN
  SET
   360SW

INVRSEQN DLOAD SQRT
   P
  PDDL DMP  # 0D=SQRT(P) (+2) PL AT 2
   SNTH
   COGA
  SL1 PDDL  # 2D=SNTH COGA (+5) PL AT 4
   CSTH
  SR4 DAD
   D1/32
  DSU DMP  # PL AT 2,0
  NORM BDDV
   X1
   SNTH
  SLR* ABS  # NOTE: NEAR 360 CASE TREATED DIFFERENTLY
   0 -5,1
  PUSH DSQ  # 0D=1/W (-1) PL AT 2
  STODL 34D
   D1/16
1/WLOOP  PUSH DSQ  # 2D=G (+4) PL AT 4
  RTB PDDL  # PL AT 7
   TPMODE
   R1A
  DMP SR4
   34D
  TAD   # PL AT 4
  BMN SQRT
   INFINITY
  DAD   # PL AT 2
  TIX,2 NORM
   1/WLOOP
   X1
  BDDV
  SLR* GOTO  # PL AT 0
   0 -7,1
   POLYCOEF

TRUE360X DLOAD BMN
   R1A
# Page 1192
   INFINITY
  SQRT NORM
   X1
  BDDV SL*
   2PISC
   0 -3,1
  DSU PUSH  # 0D=2PI/SQRT(R1A) -X PL AT 0,2
  GOTO
   XCOMMON
INFINITY SETPD BOV  # NO SOLUTION EXISTS SINCE CLOSURE THROUGH
   0  # INFINITY IS REQUIRED
   OVFLCLR
OVFLCLR  SET RVQ
   INFINFLG

# Page 1193
LAMBERT  STQ SETPD
   RTNLAMB
   0D
  BOV
   +1
  CLEAR VLOAD*
   SOLNSW
   MUTABLE,1
  STODL 1/MU
   TDESIRED
  DMPR
   BEE19
  STORE EPSILONL
  SET VLOAD
   SLOPESW
   R1VEC
  PDVL CALL  # 0D=R1VEC (+29 OR +27) PL AT 6
   R2VEC  # MPAC=R2VEC (+29 OR +27)
   GEOM
  STODL SNTH  # 0D=CSTH (+1) PL AT 2
   MAGVEC2
  NORM PDDL
   X1
   R1
  SR1 DDV  # PL AT 2
  SL* PDDL  # DXCH WITH 0D, 0D=R1/R2 (+7) PL AT 0,2
   0 -6,1
  STADR
  STORE CSTH  # CSTH (+1)
  SR1 BDSU
   D1/4
  STORE 1-CSTH  # 1-CSTH (+2)

  ROUND BZE
   360LAMB
  NORM PDDL  # PL AT 4
   X1
   0D
  SR1 DDV  # PL AT 2
  SL* SQRT
   0 -3,1
  PDDL SR  # 2D=SQRT(2R1/R2(1-CSTH)) (+5)  PL AT 4
   SNTH
   6
  DDV DAD  # PL AT 2
   1-CSTH
  STADR
  STORE COGAMAX
  BOV BMN  # IF OVFL, COGAMAX=COGUPLIM
   UPLIM  # IF NEG, USE EVEN IF LT COGLOLIM, SINCE
# Page 1194
   MAXCOGA  #  THIS WOULD BE RESET IN LAMBLOOP
  DSU BMN  # IF COGAMAX GT COGUPLIM, COGAMAX=COGUPLIM
   COGUPLIM
   MAXCOGA  # OTHERWISE OK, SO GO TO MAXCOGA
UPLIM  DLOAD
   COGUPLIM # COGUPLIM=.999511597 = MAX VALUE OF COGA
  STORE COGAMAX  # NOT CAUSING OVFL IN R1A CALCULATION
MAXCOGA  DLOAD
   CSTH
  SR DSU  # PL AT 0
   6
  STADR
  STODL CSTH-RHO
   GEOMSGN
  BMN DLOAD
   LOLIM
   CSTH-RHO
  SL1 DDV
   SNTH
  BOV
   LOLIM
MINCOGA  STORE COGAMIN  # COGAMIN (+5)
  BON SSP
   GUESSW
   NOGUESS
   TWEEKIT
   00001
  DLOAD
   COGA

LAMBLOOP DMP
   SNTH
  SR1 DSU
   CSTH-RHO
  NORM PDDL  # 0D=SNTH COGA-(CSTH-RHO) (+7+C(X1)) PL=2
   X1
   1-CSTH
  SL* DDV  # 1-CSTH (+2) PL AT 0
   0 -9D,1
  BMN BZE
   NEGP
   NEGP
  STODL P  # P=(1-CSTH)/(SNTH COGA-(CSTH-RHO)) (+4)
   COGA
  DSQ DAD
   D1/1024
  NORM DMP
   X1
   P
# Page 1195
  SR* BDSU
   0 -8D,1
   D1/32
  STODL R1A  # R1A=2-P(1+COGA COGA) (+6)
   P
  BOV CALL
   HIENERGY
   GETX
  DLOAD
   T
  STODL TPREV
   XI
  BON CALL
   INFINFLG
   NEGP  # HAVE EXCEEDED THEORETICAL BOUNDS
   DELTIME
  BOV BDSU
   BIGTIME
   TDESIRED
  STORE TERRLAMB
  ABS BDSU
   EPSILONL
  BPL RTB
   INITV
   CHECKCTR
  BHIZ CALL
   SUFFCHEK
   ITERATOR
  DLOAD BZE
   MPAC
   SUFFCHEK
  DAD
   COGA
  STORE COGA
  GOTO
   LAMBLOOP

NEGP  DLOAD BPL  # IMPOSSIBLE TRAJECTORY DUE TO INACCURATE
   DCOGA  # BOUND CALCULATION. TRY NEW COGA.
   LOENERGY

HIENERGY SETPD DLOAD  # HIGH ENERGY TRAJECTORY RESULTED
   0
   COGA  # IN OVFL OF P OR R1A, OR XI EXCEEDING 50.
  STORE COGAMIN  # THIS IS THE NEW BOUND.
COMMONLM DLOAD SR1
   DCOGA
# Page 1196
  STORE DCOGA  # USE DCOGA/2 AS DECREMENT
  BZE BDSU
   SUFFCHEK
   COGA
  STORE  COGA
  GOTO   # RESTART THIS LOOP
   LAMBLOOP

BIGTIME  DLOAD
   TPREV
  STORE T

LOENERGY SETPD DLOAD  # LOW ENERGY TRAJECTORY RESULTED
   0
   COGA  # IN OVERFLOW OF TIME.
  STORE COGAMAX  # THIS IS THE NEW BOUND.
  GOTO
   COMMONLM

SUFFCHEK DLOAD ABS
   TERRLAMB
  PDDL DMP  # PL AT 2D
   TDESIRED
   BEE17
  DAD DSU  # PL AT 0D
   ONEBIT
  BPL SETGO
   INITV
   SOLNSW
   INITV
360LAMB  SETPD SETGO  # LAMBERT CANNOT HANDLE CSTH=1
   0
   SOLNSW
   RTNLAMB

NOGUESS  SSP DLOAD
   TWEEKIT
   20000
   COGAMIN
  SR1 PDDL  # PL AT 2
   COGAMAX
  SR1 DAD
  STADR   # PL AT 0
  STORE COGA
  STORE DCOGA
  GOTO
# Page 1197
   LAMBLOOP

LOLIM  DLOAD GOTO
   COGLOLIM # COGLOLIM=-.999511597
   MINCOGA

INITV  DLOAD NORM
   R1
   X1
  PDDL SR1  # PL AT 2
   P
  DDV   # PL AT 0
  SL* SQRT
   0 -4,1
  DMP SL1
   ROOTMU
  PUSH DMP  # 0D=VTAN (+7) PL AT 2
   COGA
  SL VXSC
   5
   UR1
  PDDL   # XCH WITH 0D PL AT 0,6
  VXSC VSL1
   UN
  VXV VAD  # PL AT 0
   UR1
  VSL1
  STORE VVEC
  SLOAD BZE
   VTARGTAG
   TARGETV
  GOTO
   RTNLAMB

TARGETV  DLOAD CALL
   MAGVEC2
   LAMENTER
  STORE VTARGET
  GOTO
   RTNLAMB

# Page 1198
TIMERAD  STQ SETPD  # PL AT 0
   RTNTR
   0
  BOV
   +1
  VLOAD PDVL  # PL AT 6
   RVEC
   VVEC
  CALL
   PARAM
  BOV DLOAD  # PL AT 0
   COGAOVFL
   D1/32
  DSU DMP
   R1A
   P
  SQRT DMP
   COGA
  SL4 VXSC
   U2
  PDDL DSU  # PL AT 6
   D1/64
   R1A
  VXSC VSU  # PL AT 0
   UR1
  VSL4 UNIT
  BOV
   CIRCULAR
  PDDL NORM  # 0D=UNIT(ECC) (+3) PL AT 6
   RDESIRED # 35D=ECC (+3)
   X1
  PDDL DMP  # PL AT 8
   R1
   P
  SL* DDV  # PL AT 6
   0,1
  DSU DDV
   D1/16
   36D  # 36D=ECC (+3)
  STORE COSF
  BOV DSQ
   BADR2
  BDSU BMN
   D1/4
   BADR2
  SQRT SIGN
   SGNRDOT
  CLEAR
   APSESW
# Page 1199
TERMNVEC VXSC VSL1
   UN
  VXV PDVL  # VXCH WITH 0D PL AT 0,6
   0D
  VXSC VAD  # PL AT 0
   COSF
  VSL1 PUSH  # 0D=U2 PL AT 6

  DOT DDV  # LIMITS RESULT TO POSMAX OR NEGMAX
   UR1
   DP1/4
  SR1 BOV  # SCALE BACK DOWN TO NORMAL
   +1  # CLEAR OVFIND IF SET
  STOVL CSTH  # CSTH (+1)
   UR1
  VXV VSL1
  DOT SL1
   UN
  STODL SNTH  # SNTH (+1)
   P
  CALL
   GETX
  CLRGO
   SOLNSW
   COMMNOUT

CIRCULAR SETPD SETGO
   0
   SOLNSW
   ABTCONIC

BADR2  DLOAD SIGN
   LODPHALF
   COSF
  STODL COSF
   KEPZERO
  SETGO
   APSESW
   TERMNVEC

# Page 1200
APSIDES  STQ SETPD  # PL AT 0
   RTNAPSE
   0D
  BOV
   +1
  VLOAD PDVL  # PL AT 6
   RVEC
   VVEC
  CALL
   PARAM
  BOV   # PL AT 0
   GETECC
GETECC  DMP SL4
   R1A
  BDSU SQRT
   D1/64
  STORE ECC
  DAD PDDL  # PL AT 2
   D1/8
   R1
  DMP SL1
   P
  DDV   # PL AT 0
  PDDL NORM  # 0D=RP (+29 OR +27) PL AT 2
   R1A
   X1
  PDDL SL*  # PL AT 4
   R1
   0 -5,1
  DDV DSU  # PL AT 2,0
  BOV BMN
   INFINAPO
   INFINAPO
  GOTO
   RTNAPSE
INFINAPO DLOAD GOTO  # RETURNS WITH APOAPSIS IN MPAC, PERIAPSIS
   LDPOSMAX
   RTNAPSE  # THAT PL IS AT 0.

# Page 1201
ABTCONIC EXIT
  TC POODOO
  OCT 00607

# Page 1202
LDPOSMAX EQUALS  LODPMAX  # DPPOSMAX IN LOW MEMORY.

# ERASABLE ASSIGNEMENTS

# KEPLER SUBROUTINE

#  INPUT -
# RRECT ERASE  +5
# VRECT ERASE +5
# TAU. ERASE +1
# XKEP ERASE +1
# TC ERASE +1
# XPREV ERASE +1
1/MU  EQUALS 14D
ROOTMU  EQUALS 16D
1/ROOTMU EQUALS 18D

#  OUTPUT -
# RCV ERASE +5
# VCV ERASE +5
# RC ERASE +1
# XPREV ERASE +1

#  DEBRIS -
ALPHA  EQUALS 8D
XMAX  EQUALS 10D
XMIN  EQUALS 12D
X  EQUALS 20D
XI  EQUALS 24D
S(XI)  EQUALS 26D
XSQC(XI) EQUALS 28D
T  EQUALS 30D
R1  EQUALS 32D
KEPC1  EQUALS 34D
KEPC2  EQUALS 36D
# DELX ERASE +1
# DELT ERASE +1
# URRECT ERASE +5
# RCNORM ERASE +1
# XPREV EQUALS XKEP


# LAMBERT SUBROUTINE
#
#  INPUT -
# R1VEC  ERASE +5
# R2VEC ERASE  +5
# TDESIRED ERASE +1
# GEOMSGN ERASE +0
# GUESSW # 0 IF COGA GUESS AVIABLE, 1 IF NOT
# Page 1203
# COGA ERASE +1 # INPUT ONLY IF GUESS IS ZERO.
# NORMSW # 0 IF UN TO BE COMPUTED, 1 IF UN INPUT
# UN ERASE +5 # ONLY USED IF NORMSW IS 1
# VTARGTAG ERASE +0
# TWEEKIT EQUALS 40D # ONLY USED IF GUESSW IS 0

#  OUTPUT -
# VTARGET ERASE +5 # AVAILABLE ONLY IF VTARGTAG IS ZERO.
# V1VEC EQUALS MPAC

#  DEBRIS -
# RTNLAMB ERASE +0
# U2 ERASE +5
# MAGVEC2 ERASE +1
# UR1 ERASE +5
# R1 EQUALS 31D
# UN ERASE +5
# SNTH ERASE +1
# CSTH ERASE +1
# 1-CSTH ERASE +1
# CSTH-RHO ERASE +1
COGAMAX  EQUALS 14D # CLOBBERS 1/MU
COGAMIN  EQUALS 8D
DCOGA  EQUALS 12D
# TWEEKIT EQUALS 40D
# P ERASE +1
# COGA ERASE +1
# R1A ERASE +1
# X EQUALS 20D
# XSQ EQUALS 22D
# XI EQUALS 24D
# S(XI) EQUALS 26D
# XSQC(XI) EQUALS 28D
# T EQUALS 30D
# KEPC1 EQUALS 34D
# KEPC2 EQUALS 36D
# SLOPESW
# SOLNSW

#  OTHERS -
# RVEC EQUALS R1VEC
# VVEC ERASE +5
# COGAFLAG
# RVSW
# INFINFLG
# APSESW
# 360SW
# RTNTT EQUALS RTNLAMB
# ECC ERASE +1
# RTNTR EQUALS RTNLAMB
# Page 1204
# RTNAPSE EQUALS RTNLAMB
# R2 EQUALS MAGVEC2
COSF  EQUALS 24D
# RTNPRM ERASE +0
# SCNRDOT ERASE +0
# RDESIRED ERASE +1


# ITERATOR SUBROUTINE
# ORDERSW
MAX  EQUALS 14D  # CLOBBERS 1/MU
MIN  EQUALS 8D
# INDEP ERASE +1
DELINDEP EQUALS 12D
ITERCTR  EQUALS 22D
DEP  EQUALS 30D
# DELDEP ERASE +1
# DEPREV ERASE +1
TWEEKIT  EQUALS 40D


# MORE KEPLER
# EPSILONT ERASE +1


# MORE LAMBERT
# TERRLAMB EQUALS DELDEP
# TPREV EQUALS DEPREV

# EPSILONL EQUALS EPSILONT +2 # DOUBLE PRECISION WORD