# L-12 MCS 507 Mon 18 Sep 2023 : simpson4pi1.py

"""
Application of the Simpson rule on one process,
as reference to see if multiprocessing will speedup
this computation.
"""

from scipy.integrate import simps
from scipy import sqrt, linspace, pi
x = linspace(0, 1, 10**8)
y = sqrt(1-x**2)
I = 4*simps(y, x)
print(I, abs(I - pi))