exa2_63.m

来自「有关matlab r2007的实用教程」· M 代码 · 共 5 行

A=[12 -3 3; -16 3 -1; 1 1 1];    %系数矩阵
b=[15 -13 6]';    %右端向量
[L,U]=lu(A)    %对方程组进行LU分解
x=U\(L\b)     %根据LU矩阵求解方程组的解