计算方法复合梯形公式C语言代码(共2页).doc

合肥大学方名 0904013027#include#include#include#define EPS 1e-6double f(double n);double AutoTrap(double a,double b,double c);double f(double n)return 4/(1+n*n);double AutoTrap(double a,double b,double c)int i,j,k;double s=0,T,t,h;T=(b-a)*(f(a)+f(b)/2;printf(T(1)=%fn,T);for(i=1;i+) s=0; k=pow(2,i); /新的步长for(j=1;j=pow(2,i-1);j+)h=(2*j-1)*(b-a)/k;