Programa Cramer

!PROGRAMA PARA CALCULAR LAS RAICES DE UN SISTEMA
!DE ECUACIONES LINEALES 3X3 POR EL METODO DE CRAMER

PROGRAM CRAMERDIMENSION A(3,3),AUX(3,3),B(3),X(3),DetX(3)
N=3
PRINT*,"TECLEA LOS VALORES DE LOS COEFICIENTES DEL SISTEMA "
PRINT*,"PORRENGLONES DE IZQUIERDA A DERECHA"
READ*,((A(I,J),J=1,N,1),I=1,N,1)
PRINT*,"TECLEA LOS VALORES DE LOS TERMINOS INDEPENDIENTES"READ*,(B(I),I=1,N,1)
PRINT*,"LA MATRIZ DEL COEFICIENTE ES"
DO I=1,N,1
WRITE (6,10)(A(I,J),J=1,N,1)
10 FORMAT (15XF7.3)
END DOCALL DETE (A,DETA)
DO K=1,N
DO I=1,N
DO J=1,N
AUX(I,J)=A(I,J)
END DO
END DO
DO I=1,N
AUX(I,K)=B(I)
END DO
CALLDETE(AUX,DetB)
DetX(K)=DetB
END DO
DO I=1,N
X(I)=DETX(I)/DetA
END do
PRINT*," LA SOLUCION EL SISTEMA ES"
DO I=1,N,1WRITE(6,20)I,DETX(I)
20 FORMAT (1(/),20X,"Det",I1,")=",F9.5)
END DO
PRINT*,"Det(A)=",DetA
PRINT*,"LA SOLUCION DEL SISTEMA ES"
DOI=1,N,1
WRITE(6,30)I,X(I)
30 FORMAT (1(/),20X,"X(" ,I1,")=",F9.5)
END DO
END PROGRAM!**************************************************************************
SUBROUTINE DETE(A,DET)
DIMENSION A(3,3)P=((A(1,1)*A(2,2)*A(3,3))+(A(1,2)*A(2,3)*A(3,1))+(A(1,3)*A(2,1)*A(3,2)))
S=((A(1,3)*A(2,2)*A(3,1))+(A(1,1)*A(2,3)*A(3,2))+(A(1,2)*A(2,1)*A(3,3)))
DET=P-S

END SUBROUTINE