>>> import numpy as np
>>> from scipy.integrate import quad
>>> x = np.sort(np.random.randn(150) * 4 + 4).clip(0, 5)  # 0~5 사이의 랜덤한 값 생성
>>> func = lambda x: np.sin(x) * np.cos(x**2) + 1         # 적분할 함수
>>> y = func(x)

>>> fsolution = quad(func, 0, 5)              # 함수 quad로 적분하기
>>> print(fsolution)
(5.100345067540932, 1.2589916699365199e-08)   # 결과의 첫 번째 값이 적분 값