matlab输入输出参数都是矩阵,编制一个matlab函数文件，输入参数是一个m*n的矩阵，输出......

function [A,jb] = rref(A,tol)

%RREF   Reduced row echelon form.

%   R = RREF(A) produces the reduced row echelon form of A.

%

%   [R,jb] = RREF(A) also returns a vector, jb, so that:

%       r = length(jb) is this algorithm's idea of the rank of A,

%       x(jb) are the bound variables in a linear system, Ax = b,

%       A(:,jb) is a basis for the range of A,

%       R(1:r,jb) is the r-by-r identity matrix.

%

%   [R,jb] = RREF(A,TOL) uses the given tolerance in the rank tests.

%

%   Roundoff errors may cause this algorithm to compute a different

%   value for the rank than RANK, ORTH and NULL.

%

%   Class support for input A:

%      float: double, single

%

%   See also RANK, ORTH, NULL, QR, SVD.

%   Copyright 1984-2005 The MathWorks, Inc.

%   $Rev