实现用python的optimize库的fsolve对非线性方程组进行求解。
#coding=utf-8
#Scrip的使用,提供了真正的矩阵运算
#求解线性方程组2x1-x2^2=1,x1^2-x2=2
from scipy.optimize import fsolve #导入求解方程组的函数
def f(x): #定义要求解的方程组
x1 = x[0]
x2 = x[1]
return [2*x1 - x2**2 - 1, x1**2 -x2 -2]
result = fsolve(f,[1,1]) #出入初值[1,1]并求解
print(result) #输出结果