#求解非線性方程組2x1-x2^2=1,x1^2-x2=2
from scipy.optimize import fsolve #匯入求解方程組的函式
from scipy import integrate
def func(x:list):
    x1=x[0]
    x2=x[1]
    xy=[2*x1-x2**2-1,x1**2-x2-2]
    result=fsolve(xy,x)  #輸入初始值，並求解
    return result
def g(y):
    return (1-y**2)**0.5
pi_2,err=integrate.quad(g,-1,1)  #積分結果和誤差
print(pi_2*2)

相關推薦