ea564.m

来自「一本关于工程线性代数及其matlab编程方面的好书。」· M 代码 · 共 11 行

% 《工程线性代数（MATLAB版）》第五章例题5.6.4程序ea564
%
A=input('输入系数矩阵A'),
[U0,ip]=rref(A)  
[m,n]=size(A);
is=1:n; is(ip)=[],			% 根据ip求is的语句
lp=length(ip),ls=length(is)
X=zeros(n,n-lp);
X(ip,:)=-U0(1:lp,is);
X(is,:)=eye(ls)