单纯形法: #导入包 from scipy import optimize import numpy as np #确定c,A,b,Aeq,beq c = np.array([115,90]) A = np.array([[10,20],[4,16],[15,10]]) b = np.array([200,128,220]) #Aeq = np.array([[1,-1,1]]) #beq = np.array([2]) #求解 res = optimize.linprog(-c,A,b) print(res) 输出结果: 大M法: #导入包 from scipy import opt