PDP-10 Archive: decus/20-0025/detex.for from decus_20tap1

Trailing-Edge - PDP-10 Archives - decus_20tap1_198111 - decus/20-0025/detex.for

There is 1 other file named detex.for in the archive. Click here to see a list.

C	THE PROGRAM DETED.F4 IS A SAMPLE PROGRAM, WHEN USED IN
C	CONJUNCTION WITH DETEX(LOADED AND SAVED AS DETED.SAV)
C	WILL PRODUCE A RUNNING EXAMPLE OF THE USE OF DETEX.
      FUNCTION DETE(A,NARG)
      DIMENSION A(25,25)
1     N=NARG
      SIGN=1.
      NMIN1=N-1
      IF (NMIN1)  5001, 401, 11
11    DO 391 J1=1,NMIN1
101   TEMP=0.
      DO 121 J2=J1,N
	IF(ABS(A(J2,J1))-TEMP)121,111,111
  111	TEMP=ABS(A(J2,J1))
      IBIG=J2
121   CONTINUE
      IF (TEMP)  5001, 131, 201
131   DETE=0.
      GO TO 5001
201   IF (J1-IBIG)  211, 301, 211
211   DO 221 J2=J1,N
      TEMP=A(J1,J2)
      A(J1,J2)=A(IBIG,J2)
221   A(IBIG,J2)=TEMP
      SIGN=-SIGN
301   N1=J1+1
      DO 321 J2=N1,N
      TEMP=A(J2,J1)/A(J1,J1)
      DO 321 J3=N1,N
321   A(J2,J3)=A(J2,J3)-A(J1,J3)*TEMP
391   CONTINUE
401   DETE=1.
      DO 421 J1=1,N
421   DETE=DETE*A(J1,J1)
      DETE=DETE*SIGN
5001  RETURN
      END