sekibun - 経済統計練習帳

import numpy as np
from scipy import integrate
def computePi(x):
    return 4/(1+x**2)
z=integrate.quad(computePi, 0, 1)
print(z)