牛顿方法不起作用-未来价值增长年金

Question

我对牛顿法的代码是不同的。

在试图以“相当接近的猜测”开始之后，我试图使用牛顿法在增长的年金到期（立即）公式的终值中求解未知的利率（高精度）。但是，我的代码是发散的（即，如果正确运行，下面的答案应该为10％。公式为fv =（1 + i）* pp *（（1 + i）** n-（1 + g）** n）/（ig）。

fv = 735.30
pp = 100
i = .095 #starting guess
g = .05
n = 5

def newton_raphson_method(fv,pp,i,g, n):
    newton_raphson_i = i
    for num in range(1,15):
        newton_raphson_i = i - (((1+i)*pp*((1+i) ** n - (1+g) ** n ) - fv * i + fv *g) / ((n+1) * pp *(1+i)**n - fv))
        i = newton_raphson_i
newton_raphson_method(fv,pp,i,g, n)