Уравнение, которое вы пытаетесь решить, является (или может быть переставлено) полиномом порядка 37:
37
125060⋅r⋅(r + 1)
────────────────── - 4000
37
(r + 1) - 1
SymPy изо всех сил пытается найти аналитическое решение, но в целом аналитические решения гарантированно существуют только до 4-го порядка (и даже тогда могут быть нелепо сложными).
Я думаю, что вы действительно хотите здесь числовое решение, поэтому nsolve
должен сделать эту работу:
>>> nsolve(expr,r,1)
0.00915469685511422
Как дробная месячная процентная ставка, соответствующая годовой ставке (x100x12) в размере ~ 11%.