exa2_64.m

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

A=[16 4 8; 4 5 -4; 8 -4 22];
b=[28 5 26]';
R=chol(A)    %进行Cholesky分解
x=R\(R'\b)        %根据R矩阵求解方程组的解