数值分析实验报告第七章复化梯形公式.docxVIP

数值分析 实验报告 第七章复化辛普森公式 一，题目：用复化梯形公式计算 4 dx ，试用误差余项表达式进行分析，欲保证误差不 1x2 超过指定的 ，须至少将区间 [-4,4] 多少等分？将此结果与实际计算结果进行比较。 二，源程序 1． 用逐次分半技术 # includestdio.h # define e 0.1 /* 定义函数 */ double f(double x) { double y; y=1.0/(1.0+x*x); return(y); } void main() { int n,j,k; double T,m,t,H,E; for(n=2;;n++) { m=0; t=0; for(j=1;j=n-1;j++) m=m+f(-4.0+8.0/n*j); T=(8.0/(double)(2*n))*(f(-4)+f(4)+2.0*m);/* 复化梯形 n 次公式 */ for(k=1;k=n;k++) t=t+f(-4.0+(2.0*k-1.0)*4.0/3.0); H=8.0*t/(double)n; /* 分半后的节点的矩形的面积和 */ E=((T+H)/2-T)/3.0;/* 误差 */ if(E=eE=-e) break; else continue; } printf(%d\n,n-1); printf(%1.10f\n,E); printf(%1.1