untitled3.m

来自「这是用来解线性方程的」· M 代码 · 共 12 行

A=[2 1 0 0;1 4 1 0;0 1 4 1;0 0 1 2];
[L1,U1]=luNopiv(A,0.1)
F=[-3 6 14 -2]';
x=A\F
y=U1'\F
x1=L1'\y