试值法求利率python

题目：
【问题描述】如果在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。

代码如下：

import math
def f(x):
    sum = 12 * 300 * ((1 + x / 12) ** 240 - 1) / x - 500000
    return sum
def main():
    a,b,n=map(float,input().split())
    k = 0
    n=int(n)
    a1 = f(a)
    b1 = f(b)
    for i in range(100):
        c = b - b1 * (b - a) / (b1 - a1)
        c1 = f(c)
        if c1 == 0 or abs(b - a)< 0.5 * 10 ** (-n) or abs(c1)< 0.0001:
            print(k)
            print(round(c,n))
            break
        elif c1 * a1 < 0:
            b = c
            b1 = c1
        else:
            a = c
            a1 = c1
        k=k+1

main()

运行结果：