C+ C NAME: C MidSql C PURPOSE: C CATEGORY: C CALLING SEQUENCE: subroutine MidSql( Funk, AA,BB,S,N ) C INPUTS: C OUTPUTS: C INCLUDE: C include 'midpnt_inc.h' C PROCEDURE: C See Numerical Recipes C MODIFICATION HISTORY: C JUN-1993, Paul Hick (UCSD) C- real Funk real AA real BB real S integer N save IT Func(X) = 2.*X*Funk(AA+X*X) B = sqrt(BB-AA) A = 0. include 'midpnt_inc.h' return end