human
I am working on the file "Oraux_CentraleSupelec_PSI__Juin_2019.ipynb".
The first task is:
Pour calculer l'intégrale, on utilise scipy.integrate.quad :
Can you write Python program for it?
 ... Give me a the Python program you requested?
assistant
#Here is the Python program you requested.

def In(x, n):
    def fn(x):
        return f(x) ** n
    return integrate.quad(fn, 0, 1)[0]

def Sn(x):
    return np.sum([In(Xs, n) * x**n for n in range(0, n+1)], axis=0)