B(J)=O
DO 7 1=1, NUM
7 B(J)=B(J)+Y(I)*(X(I)**(J-1))
RETURN
END
SUBROUTINE GAUJO(N)
COMMON IA(66),X(50), Y( 101), F(10), A(10,21), D(10)
DIMENSIÓN B(10,10), c(10,10)
M=2*N+1
DO 1001 J=1,N
K=J+N
IF (I,J)1002,1003,1002
1002 A(I,K)=0
GO TO 1001
DO 1001 J=I,N
1003 A(I,K)=1
1001 CONTINUÉ
WRITE(3,20)
20 FORMAT (1H1, ////, 38X'METODO DE GAUSS-JORDAN PARA
SOLUCIONAR SELS.7/)
DO 1 L=1,N
DO 1 L1=1,N
1 B(L,L1)=A (L,L1)
DO 2 I=1,N
R=A(I,I)
DO 3 K=I,M
3 A(I,K)=A(I,K)/R:
DO 2 J=l, N
RR=A(J,I)
IF(J-I)4,2,4
4 DO 5 K=I,Y
5 A(J,K)=A(J,K)-A(I,K)**RR
2 CONTINUE
DO 7 I=1,N
7 DO(I)=A(I,M)
DO 6 I=1,N
c(I,J)=O
DO 6 K=1,N
L2=K+N
6 C(I,J)=c(I,J)+A(I,L2)* B(KJ)
WR1TE(3,1004)
1004 FORMAT(//,21X’ PRODUCTO DE LA INVERSA POR EL COEFICIENTE
DE SELS DE IDENTIDAD ‘/)
DO 8 I=1,N
WRITE(3,25)(C(I,J),J=1,N)
25 FORMAT(12X,7F12,4)
8 CONTINUE
RETIRN
END