以下程序段定义一个函数,以文件名liezy.m存盘。 %列主元的高斯消去法 function x=liezy(a,b) % a为方程组的系数矩阵,b为右端项 zg=[a b]; n=length(b); x=zeros(n,1);%zg为方程组的增广矩阵 for k=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行交换 for i=k+1:n %这个循环完成消元的过程 l= 【1】 ; zg(i,k:n+1)= 【2】 ; end end b=zg(1:n,n+1);a=zg(1:n,1:n); x(n)=b(n)/a(n,n); for k=n-1:-1:1%这个循环为回代过程 x(k)=(b(k)-a(k,k+1:n)*x(k+1:n))/ 【3】 ; end end 在命令窗口调用函数: >> a=[2 2 3;4 7 7;-2 4 5]; >> b=[3;1;-7]; >> x=liezy( 【4】 ) x = 2.0000 -2.0000 1.0000