shiyan51.m

来自「用Matlab软件以及雅克比迭代和高斯－赛德尔迭代解方程组Ax=b」· M 代码 · 共 12 行

A=[10 3 1;2 -10 3;1 3 10],
%A=[0 3 1;2 -10 3;1 3 10],
b=[14 -5 14]',
x=A\b,
[L1,U1]=lu(A);
L1,U1,
A1=L1*U1,
[L2,U2,P]=lu(A);
L2,U2,P,
A2=L2*U2,
A3=inv(P)*A2