原创复合辛普生积分法

 2008-4-13 22:53  2082 4 4 分类: 工业电子

#include "stdio.h"

double composite_simpson(double(*f)(double),double a,double b,int n);
double myfun(double x);

void main()
{
double(*fun)(double);
double a,b,S4;
int n;

a=0;
b=1;
n=4;
fun=myfun;
S4=composite_simpson(fun,a,b,n);
printf("\n积分值为:%f\n",S4);

}

/********************************************************************
* 用复合辛普生法计算函数f(x)从a到b的积分值
* 输入: f--函数f(x)的指针
*       a--积分下限
*       b--积分上限
*       n--分段数
* 输出: 返回值为f(x)从a点到b点的积分值
*******************************************************************/
double composite_simpson(double(*f)(double),double a,double b,int n)
{
double s,h,x;
int i;
printf("x\t\tf(x)\t\ts\n");
s=(*f)(a)-(*f)(b);
h=(b-a)/(2*n);
x=a;
for(i=1;i<2*n;i+=2)
{
x+=h;
s+=4*(*f)(x);
printf("%f\t%f\t%f\n",x,(*f)(x),s*h/3);
x+=h;
s+=2*(*f)(x);
printf("%f\t%f\t%f\n",x,(*f)(x),s*h/3);

}
   return(s*h/3);
}

double myfun(double x)
{
double y;

y=4.0/(1+x*x);

return(y);

}

写原创有奖励！2024面包板原创奖励正在进行中