选主元交换行列消元回代交换解的分量选主元交换行列三角分解计算Pb求解下三角方程组求解上三角方程组计算Qb休息一下休息一下
全主元消去法开始分别输入A矩阵和B矩阵的元素aij(i,j=1,2,...,n)和bi(i-1,2,...,n)对i=1,2,...,n;设置一位数组od(i)存放未知数的次序对k=1,2,...,n,选主元素|aikjk|=max|aij|(k<=i,j<=n)aikjk==0?是否ik==k?否换行,akj=aikj(j=k,k+1,...,n);bk=bikjk==k?是否换列,aik=aijk(i=k+1,k+2,...,n);od(k)=od(jk)计算乘子mik=aik/akk,i=k+1,k+2,...,n是消元计算aij=aij-mik*akj,i,j=k+1,k+2,...nbi=bi-mik*bi,i=k+1,k+2,...n回代求解b...