Pythonで複雑な方程式を解きたい！

Question

前提・実現したいこと

Pythonで方程式を解きたい。
具体的には、

P = C * ((1-(1+r)**(-n)))/r + fv*(1+r)**(-n)

という等式をrについて解きたい。

発生している問題・エラーメッセージ

NotImplementedError                       Traceback (most recent call last)
<ipython-input-47-a183add313da> in <module>
      3 P, C, r, n, fv = sympy.symbols('P C r n fv')
      4 eq = sympy.Eq(P, C * ((1-1/(1+r)**n))/r + fv/(1+r)**n)
----> 5 sympy.solve(eq, r)

~\Anaconda3\lib\site-packages\sympy\solvers\solvers.py in solve(f, *symbols, **flags)
   1169     ###########################################################################
   1170     if bare_f:
-> 1171         solution = _solve(f[0], *symbols, **flags)
   1172     else:
   1173         solution = _solve_system(f, symbols, **flags)

~\Anaconda3\lib\site-packages\sympy\solvers\solvers.py in _solve(f, *symbols, **flags)
   1740 
   1741     if result is False:
-> 1742         raise NotImplementedError('\n'.join([msg, not_impl_msg % f]))
   1743 
   1744     if flags.get('simplify', True):

NotImplementedError: multiple generators [r, (r + 1)**n]
No algorithms are implemented to solve equation -C*(1 - (r + 1)**(-n))/r + P - fv*(r + 1)**(-n)

該当のソースコード

# Solve the equation for r
import sympy
from sympy import symbols

P, C, r, n, fv = sympy.symbols('P C r n fv')
eq = sympy.Eq(P, C * ((1-1/(1+r)**n))/r + fv/(1+r)**n)
sympy.solve(eq, r)

試したこと

sympyで方程式が解けることは理解したのですが、エラー吐いたためなかなか前に進まないです。自力で30分くらい格闘したので降参します。。。
pythonで方程式を解ける方がいればアドバイスをいただきたく投稿しました。

補足情報（FW/ツールのバージョンなど）

python==3.7
sympy==1.4

Answer 1

sympy にその方程式を解くプログラムが実装されていないのだと思います。

[Return a checked solution for f in terms of one or more of the

symbols. A list should be returned except for the case when a linear
undetermined-coefficients equation is encountered (in which case
a dictionary is returned).
If no method is implemented to solve the equation, a NotImplementedError
will be raised. In the case that conversion of an expression to a Poly
gives None a ValueError will be raised.](https://github.com/sympy/sympy/blob/bdac8f1c7b7d1e29338676442138aaf39a2ab27b/sympy/solvers/solvers.py#L1381-L1383)

Wolfram Alpha でも解けないですね。変数 n が難しすぎるように見えるので、自力で解くか、n に具体的な値を代入して解くか、解析的に評価するか、数値的に評価するかした方が良いと思います。