%exm041_2.m

%（1）生成矩陣 及 ，并用三種方法求解
A=gallery(5)				
A(:,1)=[];
y=[1.7 7.5 6.3 0.83 -0.082]';
x=inv(A'*A)*A'*y,
xx=pinv(A)*y,
xxx=A\y

%（2）計算三個解的范數
nx=norm(x),
nxx=norm(xx),
nxxx=norm(xxx)

%（3）比較三種解法的方程誤差
e=norm(y-A*x),
ee=norm(y-A*xx),
eee=norm(y-A*xxx)