c    basic Gaussian elimination with pivoting  
c
      parameter (n = 2)
      dimension a(n,n)
      integer p(n)
      data (a(1,j),j=1,n) /10.0e-8,1.0/
      data (a(2,j),j=1,n) /1.0,1.0/
      data (p(i),i=1,n) /2,1/
c
      print *
      print *,' Gaussian elimination with pivoting'
      print *
c
      do k=1,n-1
         do i=k+1,n      
            z = a(p(i),k)/a(p(k),k)       
            a(p(i),k) = 0.0
            do j=k+1,n    
               a(p(i),j) = a(p(i),j) - z*a(p(k),j)      
      end do
      end do
      end do
c
      do 6 i=1,n
         do 5 j=1,n
            print 7,i,j,a(i,j)
 5       continue
 6    continue
c
 7    format (3x,'a(',i2,',',i2,') =',e13.6)
      stop
      end