c语言编程定积分sinx,C语言用辛普森公式求sinx在0到π上的定积分的源程序

满意答案

chenan2242089

推荐于 2018.05.01

用梯形法估算，再用辛普森法。

fsimpf 积分函数

a,b 积分下上限，eps 精度。

#include

#include

double fsimpf(double x)

{

return sin(x) ;

}

double fsimp(double a,double b,double eps)

{

int n,k;

double h,t1,t2,s1,s2,ep,p,x;

n=1; h=b-a;

t1=h*(fsimpf(a)+fsimpf(b))/2.0;

s1=t1;

ep=eps+1.0;

while (ep>=eps)

{

p=0.0;

for (k=0;k<=n-1;k++)

{

x=a+(k+0.5)*h;

p=p+fsimpf(x);

}

t2=(t1+h*p)/2.0;

s2=(4.0*t2-t1)/3.0;

ep=fabs(s2-s1);

t1=t2; s1=s2; n=n+n; h=h/2.0;

}

return(s2);

}

void main()

{

double a,b,eps,t;

a=0.0; b=3.141592653589793238; eps=0.0000001;

a definite integral by Simpson Method.

t=fsimp(a,b,eps);

printf("%e\n",t);

printf("\n Press any key to quit...");

getch();

}

00分享举报