Как упоминал @dantechguy, ответов бесконечное множество, но это не значит, что мы не можем заставить python сказать нам это наверняка. Лучший маршрут для решения систем уравнений - sympy
. Проверьте это здесь: Sympy
Следующая команда решит вашу систему уравнений и расскажет вам для каждой переменной и сообщит вам границы каждой.
from sympy.solvers import solve
from sympy import S
x1,x2,x3 = S('x1 x2 x3'.split())
Eq = [1*x1 + 3*x2 + 5*x3-1200, x1>0, x2>0,x3>0]
sol = solve(Eq, x1),solve(Eq, x2),solve(Eq, x3)
display(sol)
Это выводит:
((0 < x1) & (0 < x2) & (0 < x3) & (x1 < oo) & (x2 < oo) & (x3 < oo) & Eq(x1, -3*x2 - 5*x3 + 1200),
(0 < x1) & (0 < x2) & (0 < x3) & (x1 < oo) & (x2 < oo) & (x3 < oo) & Eq(x2, -x1/3 - 5*x3/3 + 400),
(0 < x1) & (0 < x2) & (0 < x3) & (x1 < oo) & (x2 < oo) & (x3 < oo) & Eq(x3, -x1/5 - 3*x2/5 + 240))
Если вы работаете в jupyter, используйте следующее, чтобы все выглядело красиво набранным с помощью LATEX:
display(solve(Eq, x1))
display(solve(Eq, x2))
display(solve(Eq, x3))
введите описание изображения здесь