3.1.m

来自「学习使用各种数值方法计算数值代数问题。其中包括计算矩阵的特征值问题」· M 代码 · 共 12 行

b=[1;1];
for p=1:17
A=[0.1^p 1;2 -3];
x=Gau(a,b);
x1=[4/(2+3*0.1^p) (2-0.1^p)/(2+3*0.1^p)];
normax(p)=max(sum(abs(A')))
e(p)=norm(x1-x)
r(p)=norm(a*x'-b)
g(p)=max(abs(A(1,1)),abs(A(1,2)),abs(A(1,3)),abs(A(1,4)))
ro(p)=max(normax(p)/g(p))
end