我們利用Simpson's Rule (綠色部分的程式) 即可求得定區間的積分 :

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define F(x) (double) (2*pow(x,4)+3*pow(x,2)-15)
#define Del_X (double) 0.01

int main(void)
{
        double Area=0,x,a,b;
        char buf[12],*c;

        printf("\nStart : "); c=gets(buf); a=atof(c);
        printf("\nEnd : "); c=gets(buf); b=atof(c);

        /* Simpson's Rule */

       for(x=a;x<b;x+=Del_X)
       {
              Area+=F(x)+4*F(x+(Del_X/2))+F(x+Del_X);
       }
       Area=Area*Del_X/6;

       printf("\nArea ( %f , %f ) : %f",a,b,Area);

       return 0;
}

對於不同的函數只要更改粉紅色部分定義函數 F(x) 的程式碼即可