以下程序段定义一个函数,以文件名liezy.m存盘。%列主元的高斯消去法functionx=liezy(a,b)%a为方程组的系数矩阵,b为右端项zg=[ab];n=length(b);x=zeros(n,1);%zg为方程组的增广矩阵fork=1:n-1[Y,j]=max(abs(zg(k:n,k)));%Y是最大值,j是最大值的序号r=j+k-1;%r是最大值所在的行t=zg(k,:);zg(k,:)=zg(r,:);zg(r,:)=t;%将第k行和第r行交换fori=k+1:n%这个循环完成消元的过程l=【1】;zg(i,k:n+1)=【2】;endendb=zg(1:n,n+1);a=zg(1:n,1:n);x(n)=b(n)/a(n,n);fork=n-1:-1:1%这个循环为回代过程x(k)=(b(k)-a(k,k+1:n)*x(k+1:n))/【3】;endend在命令窗口调用函数:>>a=[223;477;-245];>>b=[3;1;-7];>>x=liezy(【4】)x=2.0000-2.00001.0000