zhunganfa.m

来自「基于matlab的一些程序」· M 代码 · 共 8 行

a=[4 -1 0 0 0;-1 4 -1 0 0;0 -1 4 -1 0;0 0 -1 4 -1;0 0 0 -1 4];
b=[6 7 9 1 1]';
[l,u,p]=lu(a);
m=inv(l);
y=m*b;
n=inv(u);
x=m*y;
x,y