数值分析中的各种公式C代码..docVIP

二分法 2.2 算法步骤 步骤1： 准备 计算f(x)在有根区间[a,b]端点处的值f(a)，f(b). 步骤2： 二分 计算f(x)在区间中点 (a+b)/2处的值 f((a+b)/2). 步骤3： 判断 若f((a+b)/2)=0，则(a+b)/2即是根，计算过程结束，否则检验； 若f((a+b)/2)f(a)0，则以(a+b)/2代替b，否则以(a+b)/2代替a. 反复执行步骤2和步骤3，直到区间[a,b]的长度小于允许误差e，此时中点(a+b)/2即为所 求近似根。 2.3 程序流程图 3 实验结果分析 #includestdio.h void main() { float a,b,e,x,c; int k=0,n=1; scanf(%f,%f,%f,a,b,e); while(n==1) { x=(a+b)/2;c=(x*x*x-x-1)*(a*a*a-a-1); if(c0) { b=x; if(b-a=e) { printf(%f,%d\n,x,k); break;} else k++;} else { if(c==0) { printf(%f,%d\n,x,k); break;} else { a=x; if(b-a=e) { printf(%f,%d\n,x,k); break;} else k++; }}}} 高斯塞德尔迭代法求方程组解 高斯主元素消去法求方程组解 2.2 算法步骤 高斯塞德尔迭代法： 步骤1：给定初值 ，精度，最大容许迭代次数M,令k=1。 步骤2：对 i=1,2,…,n 依此计算 步骤3：求出 e=，若 e，则输出 (i=1,2,..,n )，停止计算。否则执行步 步骤4：若kM, k+1k,转步骤2继续迭代。若kM,表明迭代失败，停止计算。 高斯主元素消去法： 步骤1：det1，k=1,2,…,n-1,从第2步做到第7步； 步骤2： 按列主元素 步骤3：如果 =0，则 det0，计算停止。 步骤4：若 ，转步骤5，否则换行： 步骤5：计算乘数 , (i=k_1,…,n) 步骤6：消元计算 步骤7：detdet; 步骤8：若 计算停止。否则 步骤9：回代求解 (i=n-1,…,1) 2.3 程序流程图 3 实验结果分析 高斯塞德尔迭代法： 高斯主元素消去法： 高斯塞德尔迭代法： #includeiostream.h #includeiomanip.h//控制输入的格式，如下边setw(6),1.234567控制了小数点后六位 #includestdio.h #includemath.h #define M 10 int n; void Gauss_seidel(double x[],double x0[],double a[][M],double b[],double eps,int Nmax) { int i,j,s=0; double max; while(sNmax) { for(i=0;in;i++) { x[i]=b[i]; for(j=0;jn;j++) if(j!=i) x[i]=x[i]-a[i][j]*x[j]; x[i]=x[i]/a[i][i]; } max=fabs(x[0]-x0[0]); for(i=1;in;i++) if(fabs(x[i]-x0[i])max) max=fabs(x[i]-x0[i]); if(maxeps) break; for(i=0;in;i++) x0[i]=x[i]; s++; } if(s=Nmax) cout迭代发散！\nendl; else { cout原方程组的解为：\nendl; for(i=0;in;i++) printf(X%d=%f\n,i+1,x[i]); } } void main() { double a[10][10]={{5,