Trailing-Edge
-
PDP-10 Archives
-
decuslib20-01
-
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