PROGRAM HA7_LAPACK_LGS_solve
      implicit none
      
C     Loest das LGS
C         Ax=b

C     Konstanten
      integer N,M
      parameter (N=3,M=3)

C     Variablen fuer das Hauptprogramm
      double precision A(N,M)
      double precision b(N)
      integer i,j

C     Variablen die von DGESV intern benötigt werden
      integer ipiv(N)
      integer info

      data A / 1.0D0, 2.0D0, 3.0D0, 4.0D0, 5.0D0, 6.0D0, 
     +         7.0D0, 8.0D0, 10.0D0 /
      data b / 1.0D0, 2.0D0, 3.0D0 /


      write (*,*) 'A= '      
      write (*,99) ((A(i,j), j=1,M), i=1,N)

      write (*,*)

      write (*,*) 'b= '
      do 15, i = 1,N 
            write (*, 98) b(i)
 15   continue 

      write (*,'(//)')

C     der Funktionsaufruf zum loesen des LGS
      call DGESV(N, 1, A, N, ipiv, b, N, info)

      if (info .NE. 0) then
         write (*,*) 'Fehler in DGESV, INFO=', info
      else
         write (*,*) 'Loesung x  ist in b gespeichert'
         write (*,*) 'x='
         write (*, 98) (b(i), i=1,N)

         write (*,'(//)')
         write (*,*) 'A enthaelt jetzt die LU Faktoren der permutierten'
         write (*,*) ' Matrix'
         write (*,*) 'LU= '      
         write (*,99) ((A(i,j), j=1,M), i=1,N)
      endif

 98   format (' ',E8.2E2)
 99   format (3(' ',E8.2E2,' '))

      END