已知函数ex可以展开为幂级数1+x+x2/2!+x3/3!+⋯+xk/k!+⋯。现给定一个实数x,要求利用此幂级数部分和求ex的近似值,求和一直继续到最后一项的绝对值小于0.00001。
输入格式:
输入在一行中给出一个实数x∈[0,5]。
输出格式:
在一行中输出满足条件的幂级数部分和,保留小数点后四位。
输入样例:
1.2
输出样例:
3.3201
#include<stdio.h>
#include<math.h>
double fact(int n);
int main()
{
int n;
double x,item,sum;
scanf("%lf",&x);
item = x;
sum = 1;
n = 1;
while(fabs(item)>=0.00001){
item = pow(x,n)/fact(n);
sum = sum+item;
n++;
}
printf("%.4f",sum);
return 0;
}
double fact(int n)
{
int i;
double result;
result = 1;
for(i=1;i<=n;i++){
result = result*i;
}
return result;
}
错误做法及原因:
while(fabs(item)>=0.00001){
sum = sum+item;
n++;
item = pow(x,n)/fact(n);
}
这种写法当x取得最大值时,会和结果有0.00001的误差,分析原因为:少加一项,最后一项为0.0000093,fabs(item)判断出小于精度值,这一项还是被加进去才能算对,也就是说要包含第一个小于精度值的数,这个数也就是作为整个式子的最后一项。错误做法是先判断再求和,判断不符合就不加这个数,所以不对。