Esope source files

   1 : C GQUAD     SOURCE    CHAT      05/01/13    00:19:19     5004           
   2 :       FUNCTION GQUAD(F,A,B,N)                                                   
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 :       EXTERNAL F                                                                
   6 :       DIMENSION X(1),W(1)                                                       
   7 : C                                                                               
   8 : C     ******************************************************************        
   9 : C                                                                               
  10 : C     GQUAD IS SET EQUAL TO THE INTEGRAL OF THE FUNCTION F OVER THE IN-         
  11 : C     TERVAL (A,B), AS APPROXIMATED BY THE GAUSS-LEGENDRE QUADRATURE            
  12 : C     FORMULA OF ORDER N.                                                       
  13 : C                                                                               
  13 : C     CALLS ... SUBROUTINE D106XX.                                              
  14 : C                                                                               
  15 : C     ******************************************************************        
  16 : C                                                                               
  16 :       CALL D106XX(1,F,A,B,N,X,W)                                                
  17 :       GQUAD=X(1)                                                                
  18 :       RETURN                                                                    
  19 :       END                                                                       
  20 :