复化梯形公式

复化梯形公式

复化梯形公式重在理解梯形公式的概念，计算f(X)=

#include<stdio.h>          //梯形公式
#include<math.h>
//#define f(x) (sin(x)/x)
double f(double x)
{return x>0? (sin(x)/x):1;
}
int paw(int x,int y)
{int i,a=1;if(y==0)return 1;for(i=0;i<y;i++)a=a*2;return a;
}
int main()
{int i,j,k,n;double a,b,h,sum,x,tn[100],s[100];printf("a=");scanf("%lf",&a);printf("b=");scanf("%lf",&b);printf("k=");scanf("%d",&k);for(j=0;j<=k;j++){   n=paw(2,j);h=(b-a)/n;x=a;sum=0;for(i=1;i<n;i++){x=x+h;sum=sum+f(x);}tn[j]=(h/2)*(f(a)+f(b)+2*sum);}for(i=0;i<=k;i++)printf("%d  tn(%d)=%fn",i,paw(2,i),tn[i]);return 0;
}

这里有几点需要注意的地方，第一个就是公式里的x若为0，则需要返回1，还有就是math.h中的pow(,)里面的参数须为double，因此只能自己写一个pow函数
梯形图