[ublas] How to solve systems of linear equations?

Subject: [ublas] How to solve systems of linear equations?
From: Rui Maciel (rui.maciel_at_[hidden])
Date: 2010-03-09 11:23:53

I've started using boost::ublas but unfortunately it has been a bumpy ride.  I 
intend to use this library to solve systems of linear equations but it appears 
that there isn't any documentation on lu.hpp (that header is only mentioned in 
a test case) and google shows that the documentation doesn't even mention 
lu_factorize anywhere[¹].  To make things a bit harder, the ublas code has 
practically zero comments.

So, what am I missing?  How is boost::ublas used to solve systems of linear 
equations? 

Thanks in advance,
Rui Maciel

[¹] 
http://www.google.pt/search?q=site%3Ahttp%3A%2F%2Fwww.boost.org%2Fdoc%2Flibs%2F1_41_0%2Flibs%2F+lu_factorize