4. sin(x) + cos(x2) + 1 함수를 적분하고 구간 0~5의 적분 값 찾기
>>> 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) # 결과의 첫 번째 값이 적분 값