gusucode.com > 《精通matlab 7》包括配套光盘 > 配套光盘里的程序附件/Ex-35/Ex3522.m
%例35-22:LU分解法求解下列方程组。 A=[3,1,-4,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; %系数矩阵A b=[13,-9,6,0]'; %矩阵B [L,U]=lu(A); %LU分解 x=U\(L\b) % 求解 [L,U ,P]=lu(A); % 采用LU的第2种格式分解 x=U\(L\P*b) % 求解