COMMON IA(66), X(50), Y(101), C (10), N(A10,21), D(10), V(101)

READ (2,1) NUM,M

1FORMAT(213)

DO 1000 RTM= 1,3

READ (2,2) (X(I),Y(I),I=Í,NÜM)

2FORMAT (213)

CALL AJCUR(NUM?M)

N=M+1

CALLGAUJO(N)

CALL PREGR(NUM,N,M)

1000CONTINUÉ

CALLEXIT

end

SUBROUTINE AJCUR(NUM,M)

COMMON IA(66), X(50), Y(101), C(10), A(10,2i), B(10), V(101)

N=M+1

M1=2*M+1

DO 2 J=l, MI

SUM=0

DO 3 1=1, NUM

3SUM=SUM + X(I)**(J-1)

IF(J-M-1)4,4,5

4l=i

L1=J

GO TO 6

5L=J-M

L1=M+1

6DO 2 I=L, L1

K=J-I+1

A(I,K)=SUM

DO 7 J=l, N

B(J)=O

DO 7 1=1, NUM

7B(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

1002A(I,K)=0

GO TO 1001

DO 1001 J=I,N

1003A(I,K)=1

1001CONTINUÉ

WRITE(3,20)

20FORMAT (1H1, ////, 38X'METODO DE GAUSS-JORDAN PARA SOLUCIONAR SELS.7/)

DO 1 L=1,N

DO 1 L1=1,N

1B(L,L1)=A (L,L1)

DO 2 I=1,N

R=A(I,I)

DO 3 K=I,M

3A(I,K)=A(I,K)/R:

DO 2 J=l, N

RR=A(J,I)

IF(J-I)4,2,4

4DO 5 K=I,Y

5A(J,K)=A(J,K)-A(I,K)**RR

2CONTINUE

DO 7 I=1,N

7DO(I)=A(I,M)

DO 6 I=1,N

c(I,J)=O

DO 6 K=1,N

L2=K+N

6C(I,J)=c(I,J)+A(I,L2)* B(KJ)

WR1TE(3,1004)

1004FORMAT(//,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)

25FORMAT(12X,7F12,4)

8CONTINUE

RETIRN

END

SUBROUTINE PREGR(NUM,N,M)

COMMON IA(66),X(50),Y(101),C(10),A(10,21),B(10),V(101)

WRITE(3,9)(X(I),Y(I),I=1,NUM)

9FORMAT(////,21X’ PAREJA DE DATOS TOMADAS DEL LABORATORIO,//,123X’X’IIX’Y’//,(12X,8F12.5))

WRITE(3,10)(B(I),I=1,N)

10FORMAT(//,21X’COEFICIENTES DE LA ECUACION BUSCADAEXP=0,1,2,…,1…N’//,(12X,6E17.9))

XMIN=0.

ESCHZ=X(NUM)/100.

DO 8 I=1,101

8V(I)=V(I)+B(J+1)*(XVAR,ESCHZ)

RETURN

END