program main
      use gauss_seidel
      implicit none
      integer i, n
      real, allocatable :: A(:,:), b(:), x(:)
	  real :: prec
	  
	  read (*,*) n
	  read (*,*) prec
	  allocate(A(n,n))
	  allocate(b(n))
	  allocate(x(n))
	  read (*,*) A
	  read (*,*) b
	  read (*,*) x

      do i=1,3000000
        call solve_system(A, b, x)
        if (norm2(matmul(A, x) - b) < prec) then
          exit
        endif
      enddo
      
	  write(*,*) x
      write(*,*) matmul(A, x)
	  end