exa3.m

来自「市面上电力系统的范例 直得收藏学习的ㄧ本书」· 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