rref
Gauss-Jordan elimination.
Syntax
- R = rref(A)
- R = rref(A, tol)
- [R, p] = rref(A)
- [R, p] = rref(A, tol)
Input argument
- A - input matrix (double or single)
- tol - tolerance: scalar or max(rows, cols) _ eps(class(A)) _ norm(A, inf) (default)
Output argument
- R - a matrix: reduced row echelon form of A.
- p - a vector: nonzero pivot columns.
Description
R = rref(A) returns the reduced row echelon form of A.
[R, p] = rref(A) returns also the nonzero pivots p.
Bibliography
https://en.wikipedia.org/wiki/Gaussian_elimination
Example
A = [magic(4), eye(4)]
[R, p] = rref(A)
See also
rank.
History
Version | Description |
---|---|
1.0.0 | initial version |
Author
Allan CORNET