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

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

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

		REAL x,h

C     Boucle d'integration

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

C     Boucle de calcul
        DO WHILE (x .LT. b)
          aire = aire + h*(fn(x+h)+fn(x))/2
          x = x+h
        ENDDO
  
		END