Apollo-11/Comanche055/CONIC_SUBROUTINES.agc
2018-12-28 12:42:00 -06:00

1923 lines
46 KiB
Plaintext

# Copyright: Public domain.
# Filename: CONIC_SUBROUTINES.agc
# Purpose: Part of the source code for Colossus 2A, AKA Comanche 055.
# It is part of the source code for the Command Module's (CM)
# Apollo Guidance Computer (AGC), for Apollo 11.
# Assembler: yaYUL
# Contact: Ron Burkey <info@sandroid.org>.
# Website: www.ibiblio.org/apollo.
# Pages: 1262-1308
# Mod history: 2009-05-08 RSB Adapted from the Colossus249/ file of the
# same name, using Comanche055 page images.
# 2009-05-20 RSB Corrected: Fixed four interpreter
# instructions.
#
# 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 055 of AGC program Comanche by NASA
# 2021113-051. 10:28 APR. 1, 1969
#
# This AGC program shall also be referred to as
# Colossus 2A
# Page 1262
# 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 OUTLINED 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 1263
# 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, HYPERBOLIC, 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, XKEPNEW. 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 1264
# 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 1265
# 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 SEMIMAJOR 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 1266
# 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 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.
#
# 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 1267
# COGA IS TO BE USED BY LAMBERT
# COGA +5 DP GUESS OF COTANGNT 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.
#
# 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 ...
# LAMBERT 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
#
# OUTPUT --
# SCALE FACTOR
# VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS
# Page 1268
# -------- -------------- -----------------------
# 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 SEMIMAJOR 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 1269
# 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 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. 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.
#
# 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 1270
# 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 P00DOO 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 1271
#
# 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 1272
# 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 SIMPLE 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
# 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 1273
# 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.
# L+6 ... # 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 P00DOO 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 1274
# 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 1275
# 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 1276
# 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 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
# MAGVEC2 +7 FOR EARTH DP MAGNITUDE OF VVEC
# +5 FOR MOON
#
# PARAMETERS OF NO USE --
# SP PARAMETERS -- RTNAPSE, GOMSGN, RTNPRM, PLUS PUSHLIST LOCATIONS 0-5, 10D-11D, 14D-21D, 31D-38D.
# ADDITIONAL INTERPRETIVE SWITCHES USED -- NORMSW
SETLOC CONICS
# Page 1277
BANK
COUNT 12/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
# Page 1278
1REV SQRT BDDV
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
STCALL XMAX
DXCOMP
STOREMIN DLOAD
KEPZERO
STORE XMIN
DXCOMP DLOAD DMPR
# Page 1279
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 1280
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.
STCALL X
STORBNDS
# Page 1281
TIMEOVFL DLOAD BMN # X WAS TOO BIG
X
NEGTOVFL
STORE XMAX
CMNTOVFL DLOAD SR1
DELX
STORE DELX
BZE BDSU
KEPRTN
X
STODL X
TC
STCALL T
BRNCHCTR
NEGTOVFL STCALL XMIN
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
DMP SL1R
ROOTMU
DMP SLR*
# Page 1282
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
STCALL XPREV
KEPRTN
# Page 1283
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 1284
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 1285
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 1286
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 1287
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 1288
SETLOC CONICS1
BANK
COUNT 04/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 1289
SETLOC CONICS
BANK
COUNT 12/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
# Page 1290
SR DDV
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
STCALL P # P (+4)
RTNPRM
# Page 1291
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 1292
BANK 12
SETLOC CONICS
BANK
COUNT 12/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
# Page 1293
TC POLY
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
# Page 1294
RVQ
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 # PL AT 0
SLR* GOTO
0 -7,1
POLYCOEF
# Page 1295
TRUE360X DLOAD BMN
R1A
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 1296
LAMBERT STQ SETPD
RTNLAMB
0D
BOV
+1
SSP VLOAD*
ITERCTR
20D
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 # PL AT 4
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
# Page 1297
UPLIM # IF NEG, USE EVEN IF LT COGLOLIM, SINCE
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(XI)) 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
# Page 1298
P
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
STCALL COGA
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 1299
STORE DCOGA # USE DCOGA/2 AS DECREMENT
BZE BDSU
SUFFCHEK
COGA
STCALL COGA
LAMBLOOP
BIGTIME DLOAD
TPREV
STORE T
LOENERGY SETPD DLOAD # LOW ENERGY TRAJECTORY RESULTED
0
COGA # IN OVERFLOW OF TIME.
STCALL COGAMAX # THIS IS THE NEW BOUND.
COMMONLM
SUFFCHEK DLOAD ABS
TERRLAMB
PDDL DMP # PL AT 2D
TDESIRED
D1/4
DAD DSU # PL AT 0D
ONEBIT
BPL SETGO
INITV
SOLNSW
RTNLAMB
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
STCALL DCOGA
LAMBLOOP
# Page 1300
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 CLEAR
SOLNSW
STORE VVEC
SLOAD BZE
VTARGTAG
TARGETV
GOTO
RTNLAMB
TARGETV DLOAD CALL
MAGVEC2
LAMENTER
STCALL VTARGET
RTNLAMB
# Page 1301
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 # 36D=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 1302
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 1303
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 1304
ABTCONIC EXIT
TC P00DOO
OCT 00607
# Page 1305
SETLOC CONICS1
BANK
COUNT 04/CONIC
MUTABLE 2DEC* 3.986032 E10 B-36* # MUE
2DEC* .25087606 E-10 B+34* # 1/MUE
2DEC* 1.99650495 E5 B-18* # SQRT(MUE)
2DEC* .50087529 E-5 B+17* # 1/SQRT(MUE)
2DEC 4.902778 E8 B-30 # MUM
2DEC .203966 E-8 B+28 # 1/MUM
2DEC* 2.21422176 E4 B-15* # SQRT(MUM)
2DEC* .45162595 E-4 B+14* # 1/SQRT(MUM)
LDPOSMAX EQUALS LODPMAX # DPPOSMAX IN LOW MEMORY.
# ERASABLE ASSIGNMENTS
# 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
# Page 1306
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 AVAILABLE, 1 IF NOT
# COGA ERASE +1 # INPUT ONLY IF GUESSW 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
# Page 1307
# 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
# RTNAPSE EQUALS RTNLAMB
# R2 EQUALS MAGVEC2
COSF EQUALS 24D
# RTNPRM ERASE +0
# SGNRDOT 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
# Page 1308
# MORE LAMBERT
# TERRLAMB EQUALS DELDEP
# TPREV EQUALS DEPREV
# EPSILONL EQUALS EPSILONT +2 # DOUBLE PRECISION WORD