function LLS()
Ae=5;																			%<2>
Av=2;																			%<3>
[A,b]=LLSDATA(Ae,Av);
x=A\b;
%以下都是為獲得較清晰的顯示而編寫
S='恰定';
if Ae>Av
   S='超定';
elseif Ae<Av
   S='欠定';
end
disp(['外部獨立程序求' S '方程 Ax=b 解，其中'])
cs=blanks(Ae)';			%產生一個由方程數決定的空白字符列
ns=fix(Ae/2);				%從零方向最接近空白列長之半的數
As=cs;As(ns)='A';			
bs=cs;bs(ns)='b';
es=cs;es(ns)='=';
disp([As cs es cs num2str(A) cs cs cs cs cs bs cs es cs num2str(b)])
disp('方程的解 ')
nxs=fix(Av/2);
cxs=blanks(Av)';
xs=cxs;xs(nxs)='x';
exs=cxs;exs(nxs)='=';
disp([xs cxs exs cxs num2str(x)])