问题描述
如果在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()