Re: [Help-gsl] Ax=b

[Brian Gough]

address@hidden writes:
 > Other question: how can I solve linear system A x = b with indinite
 > solutions?
 > 
 > For example: we suppose to have:
 > x+2y+3z=3 2x+3y+8z=4
 3x+2y+17z=1
 > 
 > Thi system have infinite solutions:
 > x=-1-7t y=2+2t z=t t in R
 > 
 > Can i use gsl_linalg_LU_solve or gsl_linalg_complex_LU_solve to
 > solve this kind of linear system?

The LU routines only work with non-singular matrices.
Look at other matrix decompositions such as the SVD.

regards
-- 
Brian Gough