C	  Integration par la methode des trapezes
C	  Dominique Lefebvre  janvier 2007
C  
C     a = borne inferieure d'integration
C     b = borne superieure d'integration
C     n = nombre de pas
C     aire = surface retournee

		SUBROUTINE IntTrapezes (fn,a,b,n,aire)

		REAL a,b,fn,aire
		INTEGER n
		EXTERNAL fn

		REAL h,app
        INTEGER i

C     Boucle d'integration

C	  Initialisation des variables
    	aire = 0
        h = (b-a)/n

C     Calcul de l'aire approximative du trapeze f(a) f(b)
 		app = (h/2)*(fn(a)+fn(b))

C     Boucle de calcul
        DO i=1, n-1
          aire = aire + fn(a+i*h)
        ENDDO

C     Calcul final de l'aire
        aire = app + aire*h
   
		END