复合梯形公式c语言程序.doc

《复合梯形公式c语言程序.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、#include"stdio.h"#includeintmain(){inta,b,n;doubleS1,S2,r=1.0;doublef(doublex);doublejifen(doublea,doubleb,intn);printf("输入所求区间：");scanf("%ld%ld",&a,&b);S1=((b-a)/2.0)*(f(a)+f(b));n=1;while(r>pow(10,-5))//判断积分值是否达到误差限制{n*=2;S2=jifen(a,b,n);r

2、=S2-S1;S1=S2;}printf("n=%ld:",n);printf("S1=%lf:",S1);return0;}doublef(doublex){doubley;y=4.0/(1+x*x);returny;}doublejifen(doublea,doubleb,intn){doublex[1000],sum=0,S;for(intk=0;k<=n;k++){x[k]=a+k*(b-a)/(n*1.0);}for(k=1;k<=n-1;k++){sum+=f(x[k]);}S

3、=(b-a)/(2*n*1.0)*(f(a)+2*sum+f(b));returnS;}