untitled2.m

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

format short
A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10];
[L1,U1]=luNopiv(A,0.1)
D=[5 0 0 0;0 0.2 0 0;0 0 2 0;0 0 0 0.5];
L=L1*(D^(1/2))
Lt=L';
a=L*Lt;
b=[14.5 20 29 32.5]';
y=L\b
x=Lt\y
x2=A\b