function[X]=rowElimination(A,b)% row principal element elimination solution% input:系数矩阵A double n*n,常数项b double n*1% output:解向量X double n*1Ab=[Ab];%double n*(n+1)[n,n]=size(A);p=eye(n);%向量p代表置换矩阵,左乘即代表交换行后原来每行的位置L=eye(n);fori=1:n-1[a,...