%N matinv %F $MIR/src/subs/lspoly.for
%D Single precision matrix inversion
%P mchw
%: math
%B
      subroutine matinv(a,b,n,d,determ)   

      implicit none
      integer n,d
      real a(d,n),b(n),determ

  This is a single precision version of matrix inversion with accom-
   panying solution of linear equations.     
   a(n,n) is the matrix to be inverted. the inverse is returned into 
   the same array. b(n) is the constant vector, i.e. a(n,n)*x(n)=b(n)
   the solution x(n) is returned into array b(n). the determinant  of
   a as determ. Arrays a(d,d) and b(n) must be dimensioned in the    
   main program.      **** note: dimension a(d,d) in main. ****     
