مطالب جالب

function [R,k]=rref(A,tol)
% reduced row echelon form
[m,n]=size(A); p=max(m,n);
if nargin<2, tol=p*eps*norm(A,inf); end
  % pad matrix with zeros to make it square
if m~=n, A(p,p)=0; end
[L,U,P]=lu(A);
k=find(abs(diag(U))>tol);
R=U(k,k)\U(k,:);
  % remove added columns if any, add correct number of zero rows
if n
به نقل از:
https://mailman.cae.wisc.edu/pipermail/help-octave/2000-May/005323.html