accurate method for a system that may oscillate more rapidly than the initial integration 
interval. A minimum step size prevents the routine from subdividing the interval indefinitely. 
The input arguments are: 
NEQS Number of dependent variables in the x array 
TIME Actual time (independent variable) 
HMAX Increment for which the solution is to be returned 
X Vector of dependent variables 
EPSE Relative error criteria specified for each component of x and used for the 
components of x less than the absolute value of A 
A Absolute error criteria 
HMIN Minimum step size allowed 
FIRST Set to zero on first call; a value of | is assigned by KUTMER on subsequent 
calls for which the error criteria are satisfied, otherwise a value of 2 is 
assigned 
Subroutine PLOT2 (F, FMIN, FMAX, NVAR, NFUN, N1, N, XO, DELX) 
Data stored in the two-dimensional array F are plotted, using the printer by subroutine 
PLOT2. As many as 26 different functions, having evenly spaced abscissa values, can be 
plotted. The output is written on Unit 6. A description of variables follows. 
F Array containing data to be plotted; the Jth point of the Ith function is 
stored in F(I,J) 
FMIN An array of minimum functional values; the minimum of the Ith function 
is stored in FMIN(1) 
FMAX Same as FMIN only for maximum values 
NVAR An array of titles for each function to be plotted 
NFUN Number of functions to be plotted 
Nl First dimension of array F 
N Number of points to be plotted 
XO First abscissa value 
DELX Abscissa increment 
Subroutine PLOTER (FX, XA, HMAX, LAMBDA, IB, NWAVE) 
The routine initializes various values required to generate printer plots and computes 
pitch-and-heave ratios. The printer plots that are generated consists of pitch-and-heave time 
histories. A description of input variables follows. 
47 
