optimize求解非线性方程组求解

SciPy的optimize模块提供了许多数值优化算法

 

import numpy as np
from scipy import optimize

z = np.array([60,20])
a = np.array([[4, 2], [0.5,1]])
b = np.array([16, 3.5])
# 表示两个解均大于0
x1_bound = x2_bound = (0, None)

res = optimize.linprog(z, A_ub=a, b_ub=-b, bounds=(x1_bound, x2_bound))

print(res)