[L U]=lu(A) %A为方阵
再问： 这个我知道，我是想问你这个LU分解的程序，可以帮帮我编写这程序吗？
再答： 建议参考该函数的m文件
再问： 找到里面的，可好像不是，你把m文件发给我下，帮帮忙，谢谢
再答： 下面是Matlab提供的源代码,楼主好好思考吧，看完之后记得采纳啊！^.^ function [L,U,P] = lu(A,opt) %Embedded MATLAB Library Function % Copyright 2002-2008 The MathWorks, Inc. %#eml eml_assert(nargin > 0, 'Not enough input arguments.'); eml_assert(isa(A,'float'), ... ['Function ''lu'' is not defined for values of class ''' class(A) '''.']); eml_lib_assert(ndims(A) == 2, 'EmbeddedMATLAB:lu:inputMustBe2D', ... 'Matrix must be 2-D.'); eml_assert(nargin == 1 || ... (ischar(opt) && (strcmp(opt,'vector') || strcmp(opt,'matrix'))), ... 'Second argument must be ''vector'' or ''matrix'''); ONE = ones(eml_index_class); m = cast(size(A,1),eml_index_class); n = cast(size(A,2),eml_index_class); if nargout
