# -*- coding: utf-8 -*-
"""
Created on Thu Feb 13 01:38:28 2014

@author: dconduche
"""

def methode_rect(f, a, b, n):
    I = 0
    pas = (b-a) / float(n)
    for k in range(n):
        I = I + f(a + k*pas)
    return pas*I

def methode_rect_alternatif(f, a, b, n):
    pas = (b-a) / float(n)
    return pas * sum([f(a + k*pas) for k in range(n)])

# Retourne : un float, 
# représentant l'intégrale de f sur [a,b] par la méthode des rectangles en n pas

# Vous pouvez aussi le coder avec des tableaux numpy (en important numpy au préalable), ou sans mémoriser le pas : l'essentiel, c'est la formule.