У меня есть стандартная проблема в экономике, которую я хотел бы решить, используя sympl's nonlinsolve.
Я хочу символически найти стационарное состояние модели и с этой целью решить нелинейную систему уравнений.
Я был рад видеть, что у sympy есть nonlinsolve, который должен решать такие проблемы, но пока я не в состоянии это сделать.
Ниже вы найдете ту часть проблемы, которая не решается сама по себе:
from sympy.solvers.solveset import nonlinsolve
from sympy import var
var('c k mc n r_k eta w alpha', Rational = True, positive=True)
eq1 = -r_k + alpha * mc * k**(-1 + alpha) * n**(1 - alpha)
eq2 = -w + mc * (1 - alpha) * k**alpha * n**(-alpha)
eq3 = c**-1 * w - n**eta
eq_system = [eq1,eq2,eq3]
nonlinsolve(eq_system , [n,k,w])
alpha и eta являются параметрами, тогда как другие символы являются переменными.Я знаю, что это рациональные положительные числа, но по причине, которую я пока не понимаю, код выдает:
ValueError: alpha: Экспонент должен быть положительным целым числом
Ни то, ни другоея понимаю, почему альфа должна быть целым числом, и это не поможет, если я объявлю альфу как таковую.
Любой намек на то, как я могу получить осмысленное решение, используя sympy?
Большое спасибо,
Thore