MathGroup Archive 2004

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Matrix Inversion


If your matrix were badly conditioned, another system would probably 
also have trouble. Often the matrix merely needs rescaling. I use the 

scaledInverse[m_?MatrixQ] :=
   With[{scale = MapIndexed[If[#2[[1]] == #2[[2]], 1/Sqrt[#1], 0] & , m, 
     scale . Inverse[scale . m . scale] . scale

Tom Burton

  • Prev by Date: Drawing diagrams easily
  • Next by Date: Re: Mandelbrot Set & Mathematica
  • Previous by thread: Matrix Inversion
  • Next by thread: simultaneous fit of multiple data sets