subroutine inicia(a,b,n)
        double precision a(n,n),b(n,n)
        do i=1,n
                do j=1,n
                        a(i,j)=1.0
                        b(i,j)=0
                enddo
        enddo
        return
        end

        subroutine confere(a,b,n)
        double precision a(n,n),b(n,n)
        do i=1,n
                do j=1,n
                        if ( b(i,j) .ne. a(j,i) ) print *,"Erro: ",i,j 
                enddo
        enddo
        return
        end