计算方法第7章.docVIP

第7章 算法的程序实现 实验一 秦九韶和对分法 1.利用秦九韶算法计算多项式的值。 例如计算3x2+2x+1，当x=-1时值为2。 (1) C程序填空 #include stdio.h void main() { float a[20],y,x; int i,n; printf(请输入多项式的次数 n:\n); 【1】 ; printf(请按降幂顺序输入%d个系数\n,n+1); for(i=n;i=0;i--) scanf(%f,a[i]); printf(请输入x\n); scanf(%f,x); 【2】 ; for(i=n-1;i=0; 【3】 ) y= 【4】 ; printf(多项式值为%f\n,y); } 2. 用对分法求出方程在区间[3,4]内的根，精度要求为10-5，方程的根为3.631981。 （1）C程序填空 #include stdio.h float f(float x) {return x*x*x-2*x*x-4*x-7;} main() {float a,b,eps=1e-5,c; scanf(%f%f,a,b); while( 【1】 ) { 【2】 ; if(f(c)==0) break; else if 【3】 b=c; else 【4】 ; } printf(root=%f\n,c); } 3. 求方程的实根的上、下界，扫描法实现根的隔离，并用对分法求出所有的实根，精度要求为10-5。此方程4个实根的近似值分别为： root=-2.246983 root=-0.554953 root=0.801944 root=2.000003 提示：实根的上下界可以输入，用扫描法实现根的隔离，然后再对隔根区间使用对分法。 （1）C程序填空 #include stdio.h float f(float x) {return x*x*x*x-5*x*x+x+2;} float df(float a,float b) {float c; /*对分法求根部分*/ 【1】 return c; } main() {float a,b,x,h,eps,c; int i=0; a=-6; b=6;h=0.1;eps=1e-5; x=a; while(x=b-h/2) { if( 【2】 ) {printf([%f,%f]\n,x,x+h); i++; c= 【3】 ; printf(x%d=%f\n,i,c); } x= 【4】 ; } } 实验二 牛顿法和弦割法 1. 用牛顿法求a的立方根，精度要求为0.0000。 #include stdio.h 【1】 main () { 【2】 ; int a; printf(input a\n); scanf(%d,a); if (a==0) {printf(a=0\n);exit(0);} x1=a; do { 【3】 ; x1= 【4】 ; } while(fabs(x0-x1)=1e-5); printf (root=%f\n,x1); } 2. 编写牛顿法方程程序。 () 在2附近的根。 (2) 在1附近的根 (3) 在1附近的根。 计算结果分别为(1) 2.374424 (2) 1.497300 (3) 0.5671 3. 用弦割法解方程，精度要求为10-4，取初值x0=0，x1=1。 #include stdio.h #include math.h main () { float x0=0,x1=1,x2; int i=1; float f2(float x ),f1(float x0,float x1); do { x2=f1(x0,x1);