①プログラムが短くなりますか?
②x,y,zにするためのアドバイスお願いします。
4x+y <= 9,x+2y >= 4,2x-3y >= -6 ----> Ax>=b
(参考)
sympyのcoeffの使い方を教えて下さい
from sympy import *
var('x y')
def my_Str2Inequality(my_Str):
if my_Str.find("<=") > 0:
tmp = my_Str.split("<=")
ans = sympify(tmp[1] + "-" + tmp[0])
elif my_Str.find(">=") > 0:
tmp = my_Str.split(">=")
ans = sympify(tmp[0] + "-" + tmp[1])
else:
ans=""
return ans
def my_xykeisu(my_Pol2):
f = symbols('f', cls=Function)
f = sympify(my_Pol2)
my_ans0=f.coeff(x**1)
my_ans1=f.coeff(y**1)
p = Poly(f, [x,y])
my_ans2=p.coeff_monomial(x**0)
my_ans=[my_ans0,my_ans1,my_ans2]
return my_ans
my_inequality=['4*x+y <= 9','x+2*y >= 4','2*x-3*y >= -6']
my_Str=my_Str2Inequality(my_inequality[0])
my_list0=my_xykeisu(my_Str)
my_Str=my_Str2Inequality(my_inequality[1])
my_list1=my_xykeisu(my_Str)
my_Str=my_Str2Inequality(my_inequality[2])
my_list2=my_xykeisu(my_Str)
A=Matrix([my_list0,my_list1,my_list2])
print("#",my_inequality)
print("#",A)
# ['4*x+y <= 9', 'x+2*y >= 4', '2*x-3*y >= -6']
# Matrix([[-4, 1, 9], [1, 2, -4], [2, -3, 6]])