各种数值积分培训教案.pptVIP

#include stdio.h double func(double x) { double y; y=x; return(y); } double hermite_gauss(int n) { double s; int i,k; static double x[]={0.26356031972,1.41340305911,3.59642577104, 7.08581000586,12.64080084428}; static double w[]={0.52175561058,0.39866681108,0.07594244968, 0.00361175868,0.00002336997}; s=0.0; for(i=0;in;i++) s=s+w[i]*func(x[i]); return(s); } main() { double s; int n; clrscr(); printf(input n=); scanf(%d,n); s=hermite_gauss(5); printf(s=%lf\n,s); } 程序运行结果： input n=5↙ s= 1.000000 5.7 龙贝格（Romberg）积分法 5.7.1 龙贝格积分法的基本思想 前面讲述的各种求积方法是插值求积的思想，而龙贝格积分法的基本思想是，使用一个诸如梯形求积法等代数精度较低的求积公式，相继以步长 和 求得定积分的两个近似结果，然后再做它们适当的线性组合，就可以得到一个代数精度更高的公式。 5.7.2 实现龙贝格积分法的基本步骤 （1）输入区间 的端点 的值,最大迭代次数 以及容许误差ε； （2）计算区间 的长度 ； （3）用梯形积分法计算积分近似值 ； （4）对 计算 对 计算 ， 如果 ，则退出循环。 （5）如果 ，则继续；否则输出无解信息，转（7）； （6）输出 的值； （7）结束。 表5.1龙贝格求积算法元素进行运算的顺序 ? 实现龙贝格积分法的N—S图，如图5.11所示。 例5.9 使用龙贝格求积公式求下列定积分的值。 #include stdio.h #include math.h #define DFS_N 20 /* 等分数 */ #define MAX_N 10 /* 最大循环次数 */ #define eps 0.00001 /* 容许误差 */ double func(double x) { double y; y=4.0/(1+x*x); return(y); } double sum(double aa,double bb,long int n) { double h,s; int i; h=(bb-aa)/n; s=0.0; for(i=1;in;i++) s+=func(aa+i*h); s=s+(func(aa)+func(bb))/2.0; return(h*s); } void romberg(double a,double b) { double s,t[MAX_N+1][2]; int i,flag=0; long int n=DFS_N,m; t[0][1]=sum(a,b,n); n*=2; for(m=1;mMAX_N;m++) { for(i=0;im;i++) t[i][0]=t[i][1]; t[0][1]=sum(a,b,n); n*=2; for(i=1;i=m;i++) t[i][1]=t[i-1][1]+(t[i-1][1]-t[i-1][0])/(pow(2,2*m)-1); if(fabs(t[m][1]-t[m-1][1])eps) { printf(t[%ld][0]=%lf\n,m,t[m][1]); flag=1; break; } } if(flag==0) prin