From: Sven Hammarling
Date: Mon, 11 Dec 2006 15:40:32 +0000
It sounds as though the LAPACK routine DGESV is what you want. It does
not compute inv(A) explicitly, but it does solve the equations A X = B.
See: www.netlib.org/lapack/index.html for information on downloading
and a link to the Users' Guide.