Scott, it seems to have bug in your code: > A=matrix(c(1,-2,3,9,-1,3,0,-4,2,-5,5,17),ncol=4,byrow=T) > ref(A) [,1] [,2] [,3] [,4] [1,] 1 -2.5 2.166667 7.833333 [2,] 0 1.0 3.000000 5.000000 [3,] 0 0.0 1.000000 2.000000
the row echelon is apparently incorrect. -- View this message in context: http://r.789695.n4.nabble.com/row-echelon-form-tp833608p2529183.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.