exa3.m

来自「solve power flows by gauss seidel method」· M 代码 · 共 5 行

A = [4  -2  -10;  2  10  -12;  -4  -6  16];
B = [-10;  32;  -16];
[L,U] = lu(A)
X = inv(U)*inv(L)*B