习题4-2 求幂级数展开的部分和 (20 分)
已知函数ex可以展开为幂级数1+x+x2/2!+x3/3!+⋯+xk/k!+⋯。现给定一个实数x,要求利用此幂级数部分和求ex的近似值,求和一直继续到最后一项的绝对值小于0.00001。
输入格式:
输入在一行中给出一个实数x∈[0,5]。
输出格式:
在一行中输出满足条件的幂级数部分和,保留小数点后四位。
输入样例:
1.2
输出样例:
3.3201
代码有一个问题,当使用dowhile时,无需加上最后一项,若使用while循环,则需补上,否则会报错
#include<stdio.h>
double fact(int a)
{
double j,sum=1;
if(a==0)
sum=1.0;
for(j=1;j<=a;j++)
{
sum*=j;
}
return sum;
}
double MC(double x,int c)
{
double mc=1;
while(c)
{
c=c-1;
mc=mc*x;
}
return mc;
}
int main()
{
int i=0;
double x,sum=0.0,b;
scanf("%lf",&x);
do
{
b=MC(x,i)/fact(i);
sum+=b;
i++;
}while(b>=0.00001);
printf("%.4lf",sum);
return 0;
}