理学试验直接法求解线性方程组.pptx

试验3解线性方程组的直接法/*DirectMethodforSolvingLinearSystems*/求解§1高斯消元法/*GaussianElimination*/?高斯消元法：思路首先将A化为上三角阵/*upper-triangularmatrix*/，再回代求解/*backwardsubstitution*/。=

Guass列主元消去法步骤基本步骤1.输入A,b,TOL2.对于k=1:N2.1选主元P2.2若|P|=TOL则停止，输出无解信息,否则2.3换行2.4消元3.回代求解4.输出解x

#includestdio.h#includemath.h#includestdlib.h#includeconio.h#defineN4#defineEPS1e-5voidGuass(doubleA[][N+1]){inti,j,k=0,I,t;doublew=0,p,b;for(k=0;kN-1;k++){/*选主元*/ p=A[k][k]; I=k; for(i=k+1;iN;i++) if(fabs(A[i][k])fabs(p)) {p=A[i][k]; I=i; }if(fabs(p)EPS){printf(thenumberiserror!i