数值计算方法实验程序源代码.docVIP

上机实验 1编程求多项式 P(x)=(((0.0625x+0.0425) x+1.912) x+2.1296在x=1.0处的值。 #include stdio.h main() {static float a[]={2.1296,1.912,0.425,0.0625}; float y; int i; float x=1.0; clrscr(); y=a[3]; for (i=2;i=0;i--) y=y*x+a[i]; printf(x=%4.2f,y=%6.4f,x,y); } 2. 二分法：方程f(x)=x3-x-1=0， 利用逐步搜索法确定一个有根区间。 #include stdio.h #include math.h #define f(x) (x*(x*x-1)-1) #define e 0.005 main() { int i=0; float x,a=1,b=1.5,y=f(a); if(y*f(b)=0) { printf(\nThe range is error!); return; } else do { x=(a+b)/2; printf(\nx%d=%6.4f,i,x); i++; if (f(x)==0) break; if(y*f(x)0) b=x; else a=x; }while(fabs(b-a)e); printf(\nx=%4.2f,x); } #include stdio.h #include math.h #define f(x) ((x*x-1)*x-1) main() { int i; float x,a=1,b=1.5,y=f(a); if(y*f(b)=0) { printf(\nThe range is error!); return; } else for(i=0;i=6;i++)/*次数限制结束*/ { x=(a+b)/2; printf(\nx%d=%6.4f,i,x); if (f(x)==0) break; if(y*f(x)0) b=x; else a=x; } printf(\nx=%4.2f,x); } 3. 迭代法： (1) 求方程f(x)=x-10x+2=0的一个根 #include stdio.h #include math.h main() { float x0,x1=1; int i=1; do {x0=x1; x1=log10(x0+2); printf(\nx%d=%6.4f,i,x1); i++; }while(fabs(x1-x0)=0.00005); printf(\nx=%6.4f,x1); printf(\nf(x)=%6.4f,fabs(x1-pow(10,x1)+2)); } (2) 求方程x=e-x在0.5附近的根。 #include stdio.h #include math.h main() { float x0,x1=0.5; int i=1; do {x0=x1; x1=exp(-x0); printf(\nx%d=%7.5f,i,x1); i++; }while(fabs(x1-x0)0.001); printf(\nx=%5.3f,x1); } 4.牛顿迭代法： (1) 求方程x=e-x在0.5附近的根。 #include stdio.h #include math.h main() { float x0,x1=0.5; int i=0; printf(\nx%d=%7.5f,i,x1); i++; do { x0=x1; x1=x0-(x0-exp(-x0))/(1+x0); printf(\nx%d=%7.5f,i,x1); i++; }while(fabs(x1-x0)=0.001); printf(\nx=%5.3f,x1); } (2)求方程f(x)=x3-2x2-4x-7=0在[3, 4]中的根的近似值 , 精度要求为ε=10-2。 #include stdio.h #include math.h main() { float x0,x1=4; int i=0; printf(\nx%d=%5.3f,i,x1); i++; do { x0=x1; x1=x0-