% exm041_1.m
randn('state',0);
A=gallery('randsvd',100,2e13,2);	%產(chǎn)生條件數(shù)為2e13的100階隨機矩陣
x=ones(100,1);				%指定真解
b=A*x;							%為使Ax=b方程一致，用A和 x 生成 b 向量。
disp('矩陣 A 的條件數(shù) ')
disp(cond(A))

disp('按任意鍵繼續(xù)！  用逆陣法求解所需時間、解的絕對誤差、方程相對誤差。 ')
pause

tic					             %啟動計時器Stopwatch Timer
xi=inv(A)*b;		       % xi 是用"求逆"法解恰定方程所得的解。
ti=toc					       %關(guān)閉計時器，并顯示解方程所用的時間。
eri=norm(x-xi)			  %解向量 xi 與真解向量 x 的范-2誤差。
rei=norm(A*xi-b)/norm(b)	%方程的范-2相對殘差

disp(' 按任意鍵繼續(xù)！  用左除法求解所需時間、解的絕對誤差、方程相對誤差。 ')
pause

tic
xd=A\b;			              %"左除"求方程解。
td=toc
erd=norm(x-xd)
red=norm(A*xd-b)/norm(b)