试值法求利率（非线性方程求解）

问题描述

如果在240个月内每月付款300美元，使用试值法在利率区间[a,b] 内，求能够满足在这240个月之后使本金和利息的总值达到50万美元（允许误差范围为0.0001美元）的利率值，精确到小数点后第d位。

输入形式

在屏幕上输入3个数，依次为利率区间左端点值a、右端点值b和精确到小数点后d位。各数间都以一个空格分隔。测试用例的输入满足：b>a>0, 1<=d<=8, d为正整数。

输出形式

输出两行数据，第一行为迭代次数，第二行为求得的利率，保留d位小数。

样例输入

0.15 0.16 8

样例输出

5

0.15753931

样例说明

输入：左端点值为0.15，右端点值为0.16，求得的利率精确到小数点后8位。输出：表示经5次迭代，求得满足条件的利率值为0.15753931。

代码

#试值法

def f(I):
    P = 300
    N = 240
    target = 500000
    return P / (I / 12) * ((1 + I / 12)**N - 1)-target

def regula(a,b,d,epsilon):
    n = 0
    delta = 0.5 * 10 ** (-d)
    for each in range(1000000):
        n+=1
        c = b - f(b)*(b-a) / ( f(b)-f(a) )
        if f(c) == 0:
            break
        elif f(a)*f(c) < 0:
            b = c
            c = b - f(b)*(b-a) / ( f(b)-f(a) )
        else:
            a = c
            c = b - f(b)*(b-a) / ( f(b)-f(a) )
        if b-a<delta:
            break
        if abs(f(c))<epsilon:
            break
    return n,c

def main():
    a, b, d = input().split(' ')
    a = float(a)
    b = float(b)
    d = int(d)
    epsilon = 1e-4
    n, ans = regula(a,b,d,epsilon)
    print(n)
    print(round(ans,d))

if __name__=='__main__':
    main()