题目来源:大工慕课 链接
作者:Caleb Sung
题目要求
求1+1/2+1/4+1/7+1/11+1/16+1/22+…的值,直到最后一项的值小于10^-5。
其中a_n=1/(n+1/a_n-1), a_0=1, n=1,2,3…
注意事项
总和的最后一项是小于10^-5的,因此在代码中我在循坏体外多写了一次。
为方便理解,我把每一项的值都列了出来,项值和结果精确到小数点后10位。
解答示范
#include<stdio.h>
#include<math.h>
int main(){
double sum=1, tmp=1;
int i;
printf("第1项:1\n");
for(i=1; tmp>pow(10, -5); i++){
tmp = 1/(i+1/tmp);
sum += tmp;
printf("第%d项:%.10lf\n", i+1, tmp);
}
i++;
tmp = 1/(i+1/tmp);
sum += tmp;
printf("\n--------------------------------\n");
printf("Result: %.10lf\n", sum);
return 0;
}
运行结果
前面的几百项太长了,这里省略