数值分析计算方法程序汇总.docVIP

（一）秦九韶法 例：已知p=2x2-3x-4，求x=1时，p=? 程序： #include stdio.h void main() { float a[50],b,x; int n,i,k; scanf(%d%f,n,x); for(i=0;i=n;i++) scanf(%f,a[i]); b=a[0]; k=1; while(k=n) { b=x*b+a[k]; k=k+1; } printf(%f,b); } 结果：-5.000000 (二）复化中矩形公式 例：求=[g(y1)+g(y2)+···+g(ym)]h的值(已知g(y)=yn/2-1e-y/2 ，n=11，h=0.1，ym=20)。 程序： #include stdio.h #include math.h main() { double y,h,s; int n,i; n=11;h=0.1; s=0;y=0; for(i=0;i=200;i++) { y=y+h; s=s+pow(y,n/2.0-1)*exp(-y/2.0); } s=s*h; printf(%f\n,s); } 结果：2266.139467 （三）复化梯形公式 例：求 程序： #include stdio.h void main() { double a,b,s,h,x; int i,n; a=-1.0;b=1.0;n=10; h=(b-a)/n;x=a;s=x*x/2; for(i=1;in;i++) { x=x+h; s=s+x*x; } s=s+b*b/2; s=s*h; printf(s=%f\n,s); } 结果：s=0.680000 （四）复化辛普森公式 例：求 程序： #include stdio.h void main() { double a,b,c,s,h,x,y; int i,n; a=0.0;b=1.0;n=10;s=0.0; h=(b-a)/n;x=a;y=x+h;c=(x+y)/2; for(i=1;i=n;i++) { s=s+x*x*x+4*c*c*c+y*y*y; x=x+h; y=y+h; c=c+h; } s=s*h/6; printf(s=%f\n,s); } 结果：s=0.250000 （五）复化高斯公式 例：求 程序： #include stdio.h #include math.h main() {double a,b,h,s,x1,x2; int i,n; a=0;b=2;n=20;s=0; h=(b-a)/n; for(i=0;in;i++) {x1=a+i*h+h/2*(1/1.732+1); x2=a+i*h+h/2*(1-1/1.732); s=s+x1*x1*x1+x2*x2*x2; } s=h/2*s; printf(s=%f\n,s); } 结果：s=4.000000 (六）二维中矩形公式 例：求 程序： #include stdio.h #include math.h main() {double a,b,c,d,x[50],y[50],hx,hy,s; int i,j,n,k; a=0.0;b=1.0;c=0.0;d=1.0; n=k=10; hx=(b-a)/n;hy=(d-c)/k; x[0]=a+hx/2;y[0]=c+hy/2;s=0.0; for(j=0;jk;j++) { for(i=0;in;i++) { s=s+x[i]*x[i]+y[j]*y[j]; x[i+1]=x[i]+hx; } y[j+1]=y[j]+hy; } s=hx*hy*s; printf(s=%f\n,s); } 结果：s=0.665000 (七）迭代法 例：求x=x2的解。 程序： #include stdio.h #includemath.h main() { double x,xl,y,yl; int i,j; x=0.5;xl=x; y=0.5;yl=y; for(i=0;;i++) {x=x*x; if(fabs(xl-x)0.0001) break; else xl=x; } for(j=0;;j++) {y=sqrt(y); if(fabs(yl-y)0.0001) break; else yl=y; } printf(x=%f,y=%f\n,x,y); } 结果：x=0.000000，y=0.999915 （八）牛顿迭代法 y=f(x),求f(x*)=0