10  COM X[103,22],M[19,19],U[19,19],Q[20],V[20],P[78]
30  COM M$[60],N$[72]
40  COM N,K,N8,K8,N9,K9,Q9,Q7,Q5,Q4,Q3,Q2,Q1
45  COM I3,I4,U9,X$[20]
50  REM:10MAY74
60  REM:>SCAT
100  DEF FNA()=INT(18.5+5*)
150  DIM C[5,35],D[450],E[450]
400  Z1=Q[8]
410  Z2=Q[9]
420  F=I3
430  L=I4
500  N1=0
510  FOR I=F TO L
520  IF X[I,K9+2]=0 THEN 560
530  N1=N1+1
540  E[N1]=(X[I,Z1]-X[N9+1,Z1])/X[N9+2,Z1]
550  D[N1]=(X[I,Z2]-X[N9+1,Z2])/X[N9+2,Z2]
560  NEXT I
7000  PRINT 
7001  IMAGE #,"  HOR. OBS. IN ROW ",3D," AN OUTLIER; PLOTTED AS "
7002  IMAGE #,"  VERT. OBS. IN ROW ",3D," AN OUTLIER; PLOTTED AS "
7005  MAT C=ZER
7010  FOR I=1 TO N1
7015  D[I]=FNA(D[I])
7020  E[I]=FNA(E[I])
7025  IF D[I]>.9 THEN 7040
7030  D[I]=1
7032  GOSUB 8000
7035  PRINT  USING 7001;J
7036  PRINT "-3.4"
7040  IF D[I]<35.1 THEN 7055
7045  D[I]=35
7048  GOSUB 8000
7050  PRINT  USING 7001;J
7051  PRINT "3.4"
7055  IF E[I]>.9 THEN 7070
7060  E[I]=1
7062  GOSUB 8000
7065  PRINT  USING 7002;J
7066  PRINT "-3.4"
7070  IF E[I]<35.1 THEN 7083
7075  E[I]=35
7078  GOSUB 8000
7080  PRINT  USING 7002;J
7081  PRINT "3.4"
7083  NEXT I
7084  PRINT 
7085  FOR I1=7 TO 1 STEP -1
7086  FOR I=1 TO N1
7087  IF E[I]<I1*5+.1 AND E[I]>I1*5-5 THEN 7090
7088  GOTO 7092
7090  C[I1*5-E[I]+1,D[I]]=C[I1*5-E[I]+1,D[I]]+1
7092  NEXT I
7094  FOR I2=1 TO 5
7095  I=I1*5-I2+1
7100  IF ABS(INT((I-18)/5)-(I-18)/5)>.01 THEN 7120
7105  PRINT  USING "#,DD";(I-18)/5
7115  GOTO 7125
7120  PRINT "  -";
7125  FOR J=1 TO 35
7130  IF C[I2,J]<.01 THEN 7160
7135  PRINT TAB(1+2*J)"";
7140  IF C[I2,J]>9.1 THEN 7155
7145  PRINT  USING "#,D";C[I2,J]
7150  GOTO 7160
7155  PRINT TAB(1+2*J)"*";
7160  NEXT J
7165  PRINT 
7170  NEXT I2
7171  MAT C=ZER
7172  NEXT I1
7175  PRINT "   - - '";
7180  FOR I=1 TO 6
7182  PRINT " - - - - '";
7184  NEXT I
7185  PRINT " - -"
7186  PRINT "      -3";
7187  FOR I=-2 TO 3
7188  PRINT "        ";
7189  PRINT  USING "#,2D";I
7190  NEXT I
7195  PRINT LIN(2)
7200  PRINT "           ","MEAN","STD. DEV."
7210  PRINT "VERT.  VAR.",X[N9+1,Z1],X[N9+2,Z1]
7220  PRINT "HORIZ. VAR.",X[N9+1,Z2],X[N9+2,Z2]
7230  PRINT "SAMPLE SIZE =";Q[11]
7240  IF F=Q4 AND L=N THEN 9998
7250  PRINT 
7260  PRINT "BUT FOR THE"Q[12]"OBSERVATIONS PLOTTED ABOVE:"
7270  PRINT 
7280  PRINT "           ","MEAN","STD. DEV."
7290  PRINT "VERT.  VAR.",Q[13],Q[15]
7300  PRINT "HORIZ. VAR.",Q[14],Q[16]
7310  PRINT "SAMPLE SIZE =";Q[12]
7320  GOTO 9998
8000  N3=0
8010  FOR J=F TO L
8020  IF X[J,K9+2]=0 THEN 8050
8030  N3=N3+1
8040  IF N3=I THEN 8060
8050  NEXT J
8060  RETURN 
9998  CHAIN "$IDA",150
9999  END 
