c语言计算级数展开式 e^x=1+x+x^2/2!+……+x^n/n!

方法一:

#include<stdio.h>
#include<conio.h>
#include<math.h>                  
float fact(int a)              //定义一个求阶乘的方法 
{
    float s = 1;
    for (int i = 1; i <= a; i++)
    {
        s = s * i;             //求一个数的阶乘,并返回
    }
    return s;
}

int main() 
{
    float x, s = 1, t = 1;
    int n, i;
    printf("请输入x和n的值:");
    scanf("%f %d", &x, &n);    //输入展开式中x和n的值

    for ( i = 1; i <= n; i++)  //求n次循环的值
    {
        t = pow(x,i) / fact(i);//调用函数
        s += t;                //求和
    }
    printf("%f\n", s);

    getch();                   //使运行框停留
    return 0;
}

方法二:

/
//实验4.3 4.5 4.6
//4.3 求出10 ~10000之内能同时被2, 3, 7整除的数并输出
//4.5 根据级数公式前20项求e,缺1(教材错误)
//4.6 根据无穷级数公式求e,e即为自然对数的底
/
#include "stdio.h"  //for printf
#include "conio.h"  //for getch
#include "math.h"

double Fact(long num);  //求整数 num 的阶乘

#define MAX(M, N)  ((M) > (N) ? (M) : (N))  //求最大值的最简洁表达

void main()
{
    //4.3 
    int n = 0;
    for(int i = 10; i <= 10000; i ++)
    {
        if( ((i % 2) == 0) && ((i % 3) == 0) && ((i % 7) == 0) )
        {
            printf("%-7d", i);

            n ++;

            if((n % 10) == 0)
            {
                printf("\n");
            }
        }
    }

    //4.5 
    //前20项目 1 + 1/2! + 1/3! + 1/4! + ... + 1/n!
    int n1 = 20;
    double sum = 0;
    for(int j = 1; j <= n1; j ++)
    {
        double njie = Fact(j);
        sum += 1.0 / njie;
    }
    printf("\n循环%d次, %.10lf", j, sum);


    //4.6
    //e^x = 1 + x + x2 / 2! + x3 / 3! + ... + x^n/n! + ...
    //前后项精度满足 1e-10
    double sum1 = 0;
    //double x = -3.;
    double x = 1.;
    double pre = 0;
    for(j = 0; ; j ++)
    {
        double njie = Fact(j);
        double cur = pow(x, j) / njie;
        sum1 += cur;

        if(j > MAX(3, fabs(x)) && fabs(pre - cur) <= 1e-10)
        //if(j > 4 && fabs(pre - cur) <= 1e-10)  //x = 5 error
        //if(j > 3 && fabs(pre - cur) <= 1e-10)  //x = 4 error
        //if(fabs(pre - cur) <= 1e-10)  //error
            break;

        pre = cur;
    }
    printf("\n循环%d次, %.10lf", j, sum1);



    getch();  //暂停
}

//num >= 13 num! 用long存不下
double Fact(long num)
{
    double s = 1.0;
    for(int i = 2; i <= num; i ++)
    {
        s *= i;
    }
    return s;
}

最后,希望我的分享对你有所帮助,别忘为我点赞哦!愿我们一起努力,奔赴梦想!敬请关注,持续更新!

  • 2
    点赞
  • 18
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

小小_扫地僧

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值