数值计算方法实验程序源代码.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)