[R] 'solve' algorithm
W. Beldman
wgbeldman at student.han.nl
Wed Jan 14 10:13:17 CET 2004
I'm looking for the algorithm of the solve function (base package). As I
understand the help text neither Linpack nor Lapack's ZESV is used (in the
simple example below) but what algorithm is? Something like Gauss Jordan
Elimination should be applicable to solve a linear system
In particular I want to be able to translate it into i.e. Delphi, but in the
example below I'm not convinced that I understand what's happening. vD seems
rather small...
TIA,
W. Beldman
vD <- matrix(0.305625, 1, 1)
> vD
[,1]
[1,] 0.305625
> v <- matrix(c(0.348125, .296875, .305625), 1, 3)
> v
[,1] [,2] [,3]
[1,] 0.348125 0.296875 0.305625
> A <- solve(vD, v)
> A
[,1] [,2] [,3]
[1,] 1.139059 0.9713701 1
More information about the R-help
mailing list