C+ C NAME: C MidExp C PURPOSE: C CATEGORY: C CALLING SEQUENCE: subroutine MidExp( 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) = Funk(-alog(X))/X B = exp(-AA) A = 0. include 'midpnt_inc.h' return end