C语言求定积分

这里以函数f=sinx为例子，首先思路：求定积分用的是元素法，把面积分为一个个的小矩形，分的越密集，则结果越正确。那现在就用这种方法，利用C语言求定积分。
下面附上代码：

#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <string.h>
#include <math.h>
using namespace std;
int main()
{
    int i;
    double a,b;
    double wide=0.001;//wide表示划分的单位宽度，wide越小结果越精确，n是积分中矩形的宽
    printf("please input a,b:\n");
    scanf("%lf%lf",&a,&b);
    int n=(b-a)/wide;//a，b为积分的上限和下限
    double s=0;//s为所有小矩形面积的和，即积分的值
    for(i=0; i<n; i++)
    {
        s+=sin(a+i*wide)*wide;//a+i*wide是横坐标的值
    }
    printf("%.2lf\n",s);
    return 0;
}

下面做一下升级，方便函数变化时，灵活运用。

#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <string.h>
#include <math.h>
using namespace std;
double fun(double x)
{
    return sin(x);
}
double jifen(double a,double b,double (*p)(double))//(*p)(double)为指向函数的指针，相当于一个函数的名字作为形参了
{
    int i;
    double wide=0.001;//wide表示划分的单位宽度，wide越小结果越精确，n是积分中矩形的宽
    dint n=(b-a)/wide;//a，b为积分的上限和下限
    double s=0;//s为所有小矩形面积的和，即积分的值
    for(i=0;i<n;i++)
    {
        s+=(*p)(a+i*wide)*wide;//a+i*wide是横坐标的值
    }
    return s;
}
int main()
{

    double a,b;
    printf("please input a,b:\n");
    scanf("%lf%lf",&a,&b);
    printf("%.2lf\n",jifen(a,b,fun));
    return 0;
}

注意：double (*p)(double)表示p是指向函数的指针，被指的函数是一个实型函数，有一个实型的形参，p是指向函数的指针变量。同以往的指针不同的是，p指向的是一个函数，代表一个函数的地址，p就代表这个函数的名字。(p)(a+widei)相当于sin(a+widei)。这样jifen函数就是一个通用求定积分的函数了。里面的函数定义随意变，而所需改变的代码很少。这样大大提高了效率