安裝
sudo pip install scipy
應用
from scipy.optimize import fsolve
# 求解非線性方程組和數值積分
# 求解2x1-x2^2=1,x1^2-x2=2
def f(x):
x1 = x[0]
x2 = x[1]
return [2*x1-x2**2-1,x1**2-x2-2]
result = fsolve(f,[1,1])
print(result)
from scipy import integrate
# 數值積分
def g(x):
return (1-x**2)**0.5
pi_2,err=integrate.quad(g,-1,1)
print(pi_2*2)