******************************************
**** Scientific Subroutine Library ****
**** for Turbo Pascal ****
******************************************
September 27, 1985

The following programs were written by Allen Miller
and appear in the book entitled "Pascal Programs For Scientists
And Engineers" which is published by Sybex, 1981.
They were originally typed and submitted to MTPUG
in Oct. 1982 Juergen Loewner,
D-4400 Muenster, Hoher Heckenweg 3, West Germany.
They have had minor corrections and adaptations for
Turbo Pascal by Jeff Weiss, 1572 Peacock Ave., Sunnyvale,
CA 94087.

Full programs have .PAS extension; library procedures and
functions have .LIB extension and are generally used with
{$I ...} INCLUDE statements.


MEANS.PAS: Program to test MEANSTD.LIB. (p.26)

MEANSTD.LIB: This routine computes the mean and the standard
deviation of a set of numbers. (p.26)

RANDOM.LIB: A random generator (0..1). Note the internal routine
with Turbo Pascal is superior and used subsequently.
(p. 29)

RANTST.PAS: Program to test RANDOM.LIB,RANDG.LIB. (p. 32)

RANDG.LIB: A random generator with gaussian distribution. (p.35)

MATR1.PAS: A matrix multiplication program. (p. 50)

DETERM.PAS: A program to calculate the determinant of a 3x3
matrix. (p. 55)

SIMQ1.PAS: A program to solve three simultaneous equations
by Cramer's Rule. (p. 67)

GAUSS.PAS: A program to perform simultaneous solution by
Gaussian elimination. (p. 75)

SOLVGJ.PAS: same as above but Gauss-Jordan algorithm. (p. 84)

GAUSSJ.LIB: Gauss-Jordan matrix inversion and solution. (p. 87)

SOLVGV.PAS: A program to perform simultan. sol. by Gauss-
Jordan elimination with (mult.) const. vect. (p.96)

SOLVIT.PAS: Solve Hilbert matrix by Gauss-Jordan elimination,
example of ill-conditioning. (p. 106)

SOLVGJ2.PAS: A program to perform simult. solution when using
more equations than unknowns (by Gauss-Jordan
elimination). (p. 111)

SOLVEC.PAS: A program to perform simultaneous solution for
complex coefficients. (p. 119)

GAUSID.PAS: A program to perform simult. sol. by Gauss-Seidel.
(p. 129)

CFIT1.PAS: A program to perform a linear least-squares curve-
fit. (p. 139)

CFIT1a.PAS: same as above but with a random generator. (p.142)

PLOT.LIB: A (printer-) plotter subroutine. (p. 147)

CFIT2.PAS: A plotting program using PLOT.LIB. (p. 153)

LINFIT1.LIB: A program which fits a straight line through n
sets of x and y pairs of points. (p. 159)

LINFIT2.LIB: Another example of the above prog. (p. 167)

CFIT4.PAS: A linear least-squares fit program. (p. 164)

TSTSORT.PAS: Program to test various sort routines. (p. 172)

SORT-B1.LIB A bubble-sort. (p. 173)

SORT-B.LIB: A variation of the above prog. (p. 176)

SORT-S.LIB: A Shell-sort procedure. (p. 178)

SORT-Q-R.LIB: A recursive Quick-sort. (p. 180)

SORT-Q-N.LIB: A nonrecursive of Quick-sort. (p. 183)

LEAST1.PAS: A program to perform a linear least-squares fit. (p. 191)

LEAST2.PAS: As above but now with Gauss-Jordan procedure. (p.203)

LEAST3.PAS: A variation of LEAST2. (p. 209)

LEAST4.PAS: Fit to heat capacity Equation. (p. 216)

LEAST5.PAS: Fit to vapor pressure equation. (p. 220)

LEAST6.PAS: Variations with different problems. (p. 226)

NEWDR.PAS: A program to solve equations by Newtons method. (p. 243)

NEWDR2.PAS: Another version of the program above. (p. 249)

NEWTON.LIB: The Newton program for a library. (p. 252)

NEWTON-L.LIB: The Newton program with an iteration counter.

TRAP1.PAS: A program for integrations by the trapezoidal
rule. (p. 264)

TRAP2.PAS: Another better version of the above. (p. 266)

TRAP3.PAS: Final improved version calling trapez.lib (p. 270)

TRAPEZ.LIB: Improved Trapezoidal with end-correction. (p. 270)

SIMP1.PAS: Another integration program now by Simpson's rule. (p. 273)

SIMP2.PAS: Program calls SIMPS.LIB. (p. 278)

SIMPS.LIB: Simpson procedure with end-correction. (p. 278)

ROMB1.PAS: Last not least a integration program by the
Romberg method. (p. 281)

ROMB3.PAS: As above but now with adjustable panels. (p. 287)

FITPOL.PAS: A program to perform linear least-squares fit to
the ratio of 2 polynomals. (p. 295)

DIFFUS.PAS: An example of the above: diffuson of Zn in Cu. (p. 302)
(A least squares fit to the linearized e-function)

NLIN3.PAS: Same as above but now with a nonlinearized
e-function. (p. 310)

ERFSIMP.PAS: The Gaussian Error Function by Simpson's rule. (p. 323)

ERFD.PAS: An infinite series expansion for the Gaussian
error function. (p. 326)

ERFD3.PAS: The Gaussian error function and its complement. (p. 330)

ERF4.PAS: An improved Gaussian error function. (p. 334)

TSTGAM.PAS: A program to test the Gamma function. (p. 341)

TSTBES.PAS: A program to test the Bessel function. (p. 344)

BESY.PAS: An evaluation of the Bessel function of the 2nd
kind. (p. 349)