习题4-2 求幂级数展开的部分和 (20 分)
已知函数ex可以展开为幂级数1+x+x2/2!+x3/3!+⋯+xk/k!+⋯。现给定一个实数x,要求利用此幂级数部分和求ex的近似值,求和一直继续到最后一项的绝对值小于0.00001。
输入格式:
输入在一行中给出一个实数x∈[0,5]。
输出格式:
在一行中输出满足条件的幂级数部分和,保留小数点后四位。
输入样例:
1.2
输出样例:
3.3201
分析:
可以展开为幂级数1+x+x2/2!+x3/3!+⋯+xk/k!+⋯,即1+ /1! + /2! +......,找其中的规律,然后主函数用来累加,子函数用来计算子项值。
代码:
#include<stdio.h>
#include<math.h>
double fun(double N,double m);//子函数
int main()
{
double N,sum=1;
scanf("%lf",&N);
int m=1;
while(fun(N,m) >= 0.00001)
{
sum = sum + fun(N,m);//计算累加和
m++;
}
sum = sum + fun(N,m);//确保最后一项加进去
printf("%.4lf\n",sum);
return 0;
}
double fun(double N,double m)
{
double fact=1;
double answer=0;
for(int i=1;i<=m;i++)
{
fact = fact * i;//阶乘
}
answer = pow(N,m)/fact;//子项值
return answer;
}