数值积分方法的C语言实现和验证（数值计算方法）

一  复合梯形公式

二  复合辛普森公式

三  C语言程序实现

#include <stdio.h>
#include <math.h>
typedef float (*funType)(float x);
float fun(float x)       //测试算例函数
{
    return x*x;
}
float SnFun(funType fp,float a,float b,int n)  //辛普森求积公式子函数
{
    float S1=0,S2=0,Sn=0,h;
    int k=0,m;
    h=(b-a)/n;
    while(k<n)
    {
        S1=S1+4*fp(a+(k+0.5)*h);
        k++;
        printf("k=%d S1=%f\n",k,S1);
    }
    k=1;
    while(k<n)
    {
        S2=S2+2*fp(a+k*h);
        k