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.17 7.83
[2,]0 1.0 3.00 5.00
[3,]0 0.0 1.00 2.00
the row echelon is apparently incorrect.
--
View this
, 2007 5:40 PM
> To: r-help@r-project.org
> Subject: Re: [R] Row-Echelon Form
>
> > I append the function below, along with some other simple
> > linear-algebra functions.
>
> I never thanked you, by the way, John; have you considered
> incurring the overhead of prod
> I append the function below, along with some other simple
> linear-algebra functions.
I never thanked you, by the way, John; have you considered incurring
the overhead of producing a formal R package?
Otherwise, it's worth its SLOC in gold.
__
R-help
3 matches
Mail list logo