我在Matlab中有这个功能:
function [L,A]=LU_factor(A,n)
L=eye(n);
for k=1:n
if (A(k,k) == 0) Error('Pivoting is needed!'); end
L(k+1:n,k)=A(k+1:n,k)/A(k,k);
A(k+1:n,:) = A(k+1:n,:) - L(k+1:n,k)*A(k,:)
end
当我使用命令行创建矩阵A = rand(4,4)
并计算[L U] = LU_factor(A,4)
时,一切都很好,我得到了理想的结果。但是当我尝试发布它时,我得到了这个错误:
使用LU_factor时出错(第2行)没有足够的输入参数。
如何在文件中发布我的作品,使其同时包含函数和示例。 (对不起,我对Matlab来说是全新的)。
答案 0 :(得分:0)
第一种方法是拥有一个主要功能,然后发布:
function main
A = rand(4);
n = 2;
[L,A]=LU_factor(A,n)
end
function [L,A]=LU_factor(A,n)
L=eye(n);
for k=1:n
if (A(k,k) == 0)
error('Pivoting is needed!');
end
L(k+1:n,k)=A(k+1:n,k)/A(k,k);
A(k+1:n,:) = A(k+1:n,:) - L(k+1:n,k)*A(k,:)
end
end
this answer已经提到了另一种方式。