Card 3 (8F10.0 Format) - Only Necessary If |A=1. 
YAXISB Height of bow acceleration axis in inches 
YAXISC Height of CG acceleration axis in inches 
PROGRAM MAIN 
Program MAIN reads all necessary input data from cards, sets up initial values, computes 
constants, calls KUTMER to determine the state variables at TIME for the period from XA 
to XE in increments of HMAX. A table state variables is created for every PTIME-th value. 
The values for \/H and 0/27 H/d are calculated and printed. If the plot option is on, a 
printer plot will be produced. 
Subroutine COMPUT(X) 
This routine computes pitch moment NL and lift force FL, excluding added mass terms, 
using values of integrals computed in subroutine FUNCT. The argument X contains the state 
vector. 
Subroutine DAUX 
This subroutine is called from KUTMER or EULER. It determines the values of m,. b, 
and bl*, based on the following equations 
hed) = Zeg = E(I) sin 8 + (1) cos 6 - r(1) 
where r(I) = r, cos k ixee + €(1) cos 0 + €(I) sin @ + ct] 
Then for 
hy (1) 0) 
hy, 
OD) = ae at) Ga 
where V(I) = -rk sin 0 ixee + &(1) cos 6 + (1) sin 6 + ct] 
If 
d(1) > b,,() tan (BC) 2/n) 
43 
