program test

      integer :: i,INFO
      real , dimension(3,3) :: A
      real , dimension(3) :: B
      integer , dimension(3) :: IPIV

      A(1,:)=(/ 1,2,3 /)
      A(2,:)=(/ 4,5,6 /)
      A(3,:)=(/ 3,8,9 /)

      B = (/1,0,-1/)

      do i = 1,3
         print *, A(i,1),A(i,2),A(i,3) 
      enddo

      print *," "
      do i = 1,3
         print *, B(i)
      enddo


      call SGESV(3,1,A,3,IPIV,B,3,INFO)


      Print *, "Solved the system.  INFO = ", INFO
      Print *, "INFO = 0 is success, INFO > 0 means A is singular", INFO
      do i = 1,3
         print *, A(i,1),A(i,2),A(i,3)
      enddo

      print *," "
      do i = 1,3
         print *, B(i)
      enddo


      end