10  COM T[64],Q[64],M2,T1,N4
20  COM N2,N3,G[10],E[64],X[64],Y[64]
30  REM ***** AEVS VERSION 1, AUGUST 5, 1971     R.J.M. TAYLOR
40  REM ***** TRANSLATED FROM 'FORTRAN PROGRAMMING FOR THE BEHAVIORAL
50  REM ***** SCIENCES ' BY D.J. VELDMAN
60  C=0
70  DIM A[64],V[64],Z[64]
80  FILES W,A
90  REM COMPUTE TRACT
100  T=0
110  FOR I=1 TO N2
120  READ #1,I*2-1
130  MAT  READ #1;A[N2]
140  T=T+A[I]
150  NEXT I
160  FOR K=1 TO N3
170  FOR I=1 TO N2
180  X[I]=Y[I]=1
190  NEXT I
200  E[K]=1
210  E1=1
220  FOR M=1 TO 25
230  FOR I=1 TO N2
240  READ #2,I*2-1
250  MAT  READ #2;V[N2]
260  V[K]=X[I]/E[K]
270  Z[I]=Y[I]/E1
280  READ #2,I*2-1
290  MAT  PRINT #2;V
300  NEXT I
310  FOR I=1 TO N2
320  X[I]=0
330  READ #1,I*2-1
340  MAT  READ #1;A[N2]
350  FOR I1=1 TO N2
360  READ #2,I1*2-1
370  MAT  READ #2;V[N2]
380  X[I]=X[I]+A[I1]*V[K]
390  NEXT I1
400  Y[I]=0
410  FOR I1=1 TO N2
420  READ #1,I1*2-1
430  MAT  READ #1;A[N2]
440  Y[I]=Y[I]+A[I]*Z[I1]
450  NEXT I1
460  NEXT I
470  E2=0
480  FOR I1=1 TO N2
490  READ #2,I1*2-1
500  MAT  READ #2;V[N2]
510  E2=E2+X[I1]*V[K]
520  NEXT I1
530  E[K]=SQR(ABS(E2))
540  NEXT M
550  E1=0
560  FOR I1=1 TO N2
570  E1=E1+X[I1]*Z[I1]
580  NEXT I1
590  E1=SQR(ABS(E1))
600  IF E2<C*C THEN 820
610  REM DEFLATE R
620  B=0
630  FOR I1=1 TO N2
640  READ #2,I1*2-1
650  MAT  READ #2;V[N2]
660  B=B+V[K]*Z[I1]
670  NEXT I1
680  D=E[K]/B
690  FOR I=1 TO N2
700  FOR J=1 TO N2
710  READ #1,I*2-1
720  MAT  READ #1;A[N2]
730  READ #2,I*2-1
740  MAT  READ #2;V[N2]
750  A[J]=A[J]-V[K]*Z[J]*D
760  READ #1,I*2-1
770  MAT  PRINT #1;A
780  NEXT J
790  NEXT I
800  NEXT K
810  GOTO 830
820  N3=K-1
830  REM COMPUTE % OF TRACE
840  FOR I=1 TO N3
850  X[I]=E[I]/T*100
860  NEXT I
870  E3=0
880  FOR K1=1 TO N3
890  E3=E3+X[K1]
900  NEXT K1
910  PRINT '10'13'10'13
920  PRINT '10'13'10'13'10'13'10'13'13
930  PRINT "TRACE =";.0001*INT(T*10000+.5)
940  PRINT '10'13
950  PRINT .001*INT(E3*1000+.5);"PCT. OF TRACE EXTRACED BY";N3;"ROOTS."
960  PRINT '10'13
970  CHAIN "$DISC2"
980  END 
