数值分析报告作业(C语言编程实现).pdfVIP

实用标准 #include stdio.h #include math.h double f(double x) { double ans; ans=exp(x); return ans; } void main() { double a=1,b=3,error=0.0001,t[20][20],h,c; int i,j,k,m,n; h=b-a; t[0][0]=h*(f(a)+f(b))/ 2; k=1; while(1) { t[0][k]=0; m=1; for(j=0;jk-1;j++) m=m*2; for(i=1;i=m;i++) t[0][k]=t[0][k]+h*f(a+(i-0.5)*h); t[0][k]=(t[0][k]+t[0][k-1])/ 2; for(j=1;j=k;j++)