Dear R,

If I have remembered correctly, a square matrix is singular if and only if
its determinant is zero. I am a bit confused by the following code error.
Can someone give me a hint?

> a <- matrix(c(1e20,1e2,1e3,1e3),2)
> det(a)
[1] 1e+23
> solve(a)
Error in solve.default(a) :
  system is computationally singular: reciprocal condition number = 1e-17


Thanks in advance!


Feng

-- 
Feng Li
Department of Statistics
Stockholm University
106 91 Stockholm, Sweden
http://feng.li/

	[[alternative HTML version deleted]]

