algoritmo

program trans
implicit none
real, dimension(:,:), allocatable:: A
real, dimension(:), allocatable:: B
real, dimension(:), allocatable:: X
integer:: erro, i, j, k, r
real::m, s
!========================================================
print*,"Entre com os coeficientes da matriz A escrevendo linha por linha"
    do i=1,4,1
    read*, (A(i,j), j=1,4,1)
    end do
print*,'==================================================='
print*,'A matriz digitada foi:'
do i=1,4,1
print*,(A(i,j), j=1,4,1)
end do
print*, 'Se a matriz digitada está certa, digite 0, caso contrário, digite 1!'
read*, r
if(r==1) then
else if (r==0) then
print*,"Entre com a matriz dos termos independentes"
    read*,(B(i), i=1,4,1)
print*,'==================================================='
print*,'Os valores digitados foram:'  
          do i=1,4,1    
    print*,B(i)
          end do
!O processo de impressão da matriz foi utilizado para confirmar se os dados !digitados estão corretos. Caso não esteja, há a opção de digitá-los novamente
!==============================================================
!Processo de triangulação ou eliminação
    do k=1,3,1
        do i=k+1,4
            do j=k+1,4
            m=(A(i,k))/(A(K,K))
            A(i,k)=0
            A(i,j)=A(i,j)-m*(A(k,j))
            B(i)=B(i)-m*(B(i))
            end do        
        end do
    end do
!Processo de resolução do sistema
X(4)=B(4)/(A(4,4))
    do k=k-1,1,-1
    s=0
        do j=k+1,4,1
        s=s+(A(k,j))*(X(j))
        X(k)=((B(k))-s)/(A(k,k))
        end do
    end do
print*, '======================================================'
print*, 'Os valores das temperatura são:'
do i=1,4,1
print*, 'temperatura',i,'=', X(i)
end do
else
print*, 'O valor digitado não corresponde ao requerido'
end if
deallocate(A,B,X)
stop
end program trans