Matrix: finding x in A.x=b



On Friday 23 November 2007 10:36, Mansur Marvanov wrote:

> I'm sorry, I forget mention, that I need to solve A.x=b by Gauss method.

> invert(A).b; <<< it's work, but is it close to Gauss?

As for Maxima, invert(A) uses the adjoint method. Gaussian methods involve 
elimination scheme (partial---in classic method with substitution---or 
full---in Gauss-Jordan algorithm), which is not the case of using matrix 
inverse explicitly.

> LU-method - it's another method with what I should work:)

LU decomposition is more convenient and useful than Gaussian methods, so it's 
not a big fault for Maxima if there are no implementation of classic 
elimination algorithms.

It's better to write your own code if you need to use Gaussian methods in 
educational works.

-- 
Alexey Beshenov <al at beshenov.ru>
http://beshenov.ru/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 189 bytes
Desc: not available
Url : http://www.math.utexas.edu/pipermail/maxima/attachments/20071124/4b9e24b9/attachment.pgp