C+ C NAME: C MidSqu C PURPOSE: C CATEGORY: C CALLING SEQUENCE: subroutine MidSqu( 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(BB-X*X) B = sqrt(BB-AA) A = 0. include 'midpnt_inc.h' return end