C     SSAR3                                                             SR3   1 
C     SMALLEST SPACE ANALYSIS-R3, A NONMETRIC METHOD FOR OFF-DIAGONAL   SR3   2 
C     SUBMATRICES AND FOR MULTIDIMENSIONAL UNFOLDING.  THIS             SR3   3 
C     RESEARCH IN NONMETRIC METHODS IS SUPPORTED IN PART BY A NSF GRANT SR3   4 
C     (GS-929) TO L. GUTTMAN AND J. C. LINGOES.  PROGRAMMED IN FORTRAN  SR3   5 
C     II BY LINGOES 12/15/65.                                           SR3   6 
C                                                                       SR3   7 
C     DECK SET-UP FOR G-L(SSAR-III) -                                   SR3   8 
C                                                                       SR3   9 
C        1.  SYSTEM ID CARD/S.                                          SR3  10 
C        2.  BINARY PROGRAM.                                            SR3  11 
C        3.  TITLE CARD (PUNCH A 1 IN COLUMN 1 AND ANY BCD TITLE IN COL-SR3  12 
C            UMNS 2-72, WHICH WILL BE PRINTED OUT FOR EACH PAGE OF OUT- SR3  13 
C            PUT).                                                      SR3  14 
C        4.  PARAMETER CARD, 9 4-COLUMN FIELDS CONTAINING THE FOLLOWING SR3  15 
C            INFORMATION SERIATUM -                                     SR3  16 
C            A)  NS= THE NUMBER OF I-SCALES OR ROWS OF OFF-DIAGONAL SUB-SR3  17 
C                MATRIX, (.LE. 33),                                     SR3  18 
C            B)  NV=THE NUMBER OF STIMULI OR THE NUMBER OF COLUMNS IN   SR3  19 
C                THE OFF-DIAGONAL SUBMATRIX (.LE.33),                   SR3  20 
C            C)  MIND=MINIMUM NUMBER OF DIMENSIONS DESIRED.  IF ZERO OR SR3  21 
C                BLANK PROGRAM WILL DETERMINE,                          SR3  22 
C            D)  MAXD=MAXIMUM NUMBER OF DIMENSIONS DESIRED (.LE. 8 OR   SR3  23 
C                THE ROW OR COLUMN ORDER OF SUBMATRIX, WHICHEVER IS     SR3  24 
C                SMALLER),                                              SR3  25 
C            E)  ISIM=ZERO OR BLANK IF DISTANCES/DISSIMILARITIES AND 1  SR3  26 
C                IF SIMILARITIES/PROXIMITIES,                           SR3  27 
C            F)  MU=0 OR BLANK IF RANKS ARE TO BE MAINTAINED OVER ROWS  SR3  28 
C                AND COLUMNS,  MU=1 IF RANKS MAINTAINED OVER ROWS ONLY  SR3  29 
C                AS IN COOMBS' MULTIDIMENSIONAL UNFOLDING ANALYSIS,     SR3  30 
C            G)  IFD=1 IF DISTANCE MATRIX IS TO BE PRINTED FOR 2 OR MORESR3  31 
C                DIMENSIONS, OTHERWISE SET TO ZERO OR LEAVE BLANK,      SR3  32 
C            H)  IFC=1 IF COORDINATES ARE TO BE PUNCHED FOR 2 OR MORE   SR3  33 
C                DIMENSIONS, OTHERWISE SET TO ZERO OR LEAVE BLANK,      SR3  34 
C            I)  IFCB=1 IF CITY-BLOCK MODEL IS TO BE USED (N. B.  IF    SR3  35 
C                THIS ALTERNATIVE IS USED MIND MUST BE .G. 0), OTHERWISESR3  36 
C                SET TO ZERO OR LEAVE BLANK FOR EUCLIDEAN MODEL.        SR3  37 
C        5.  FORMAT CARD (DESCRIBING IN F-NOTATION HOW DATA IS PUNCHED).SR3  38 
C        6.  DATA - EACH SET OF CARDS CONTAINS ALL THE DATA FOR ONE ROW SR3  39 
C            OF SUBMATRIX.  NS SETS WILL BE INPUT.                      SR3  40 
C        7.  REPEAT ITEMS 3-6 FOR ADDITIONAL RUNS.                      SR3  41 
C                                                                       SR3  42 
C     *** REFERENCES - LINGOES, J. C.  AN IBM 7090 PROGRAM FOR GUTTMAN- SR3  43 
C                                  LINGOES SMALLEST SPACE ANALYSIS-RIII.SR3  44 
C                                  BEHAV. SCI., 1966,11, 323.           SR3  45 
C                                                                       SR3  46 
C     *** NOTE COMMENT CARDS FOR PROGRAM CHANGES ON OTHER SYSTEMS.      SR3  47 
C                                                                       SR3  48 
      DIMENSION DIST(66,66), EVAL(4356), C(66,66), EVEC(4356), X(66,8), SR3  49 
     1XP1(66,9), IND(66,9), STRESS(100), FMT(18), INDEXI(33,33), INDEXJ(SR3  50 
     233,33), INDI(33,33), INDJ(33,33), PD(33,33), DP(33,33), NOTYR(33,3SR3  51 
     33), NOTYC(33,33), NTR(33), NTC(33), PHI(66)                       SR3  52 
      COMMON INDEXI,DIST,EVAL,INDEXJ,X,XP1,PD,INDI,INDJ,DP,NOTYR,NOTYC  SR3  53 
      EQUIVALENCE (DIST,EVEC), (EVAL,C), (XP1,IND)                      SR3  54 
C                                                                       SR3  55 
C     SUBROUTINES NEEDED=EIGEN,PLOT, AND MXOUT                          SR3  56 
C                                                                       SR3  57 
C     TAPE ASSIGNMENTS                                                  SR3  58 
      ITAPE=5                                                           SR3  59 
      JTAPE=6                                                           SR3  60 
C     READ TITLE, PARAMETERS, AND FORMAT                                SR3  61 
1     READ (ITAPE,196)                                                  SR3  62 
      READ (ITAPE,197) NS,NV,MIND,MAXD,ISIM,MU,IFD,IFC,IFCB             SR3  63 
      READ (ITAPE,198) (FMT(J),J=1,18)                                  SR3  64 
C     INPUT DATA MATRIX IN TRANSPOSED FORM                              SR3  65 
      DO 2 J=1,NS                                                       SR3  66 
2     READ (ITAPE,FMT) (PD(I,J),I=1,NV)                                 SR3  67 
      MU=MU+1                                                           SR3  68 
      GO TO (3,4), MU                                                   SR3  69 
3     IT=10                                                             SR3  70 
      GO TO 5                                                           SR3  71 
4     IT=30                                                             SR3  72 
5     DO 6 I=1,NS                                                       SR3  73 
      DO 6 J=1,NV                                                       SR3  74 
6     DP(I,J)=PD(J,I)                                                   SR3  75 
C     PRINT INPUT VALUES                                                SR3  76 
      WRITE (JTAPE,196)                                                 SR3  77 
      WRITE (JTAPE,199)                                                 SR3  78 
      DO 7 J=1,NS                                                       SR3  79 
      WRITE (JTAPE,200) J                                               SR3  80 
7     WRITE (JTAPE,201) (PD(I,J),I=1,NV)                                SR3  81 
      MD=66                                                             SR3  82 
      ND=9                                                              SR3  83 
C     INITIALIZE                                                        SR3  84 
      NSP1=NS+1                                                         SR3  85 
      FNV=NV                                                            SR3  86 
      FNS=NS                                                            SR3  87 
      FNVP1=NV+1                                                        SR3  88 
      NSM1=NS-1                                                         SR3  89 
      NSPNV=NS+NV                                                       SR3  90 
      NSNVM1=NSPNV-1                                                    SR3  91 
      NSNVP1=NSPNV+1                                                    SR3  92 
      FN=NS+NV-1                                                        SR3  93 
      FNR=FN+1.                                                         SR3  94 
C     ORDER ROWS OF INPUT MATRIX                                        SR3  95 
      DO 9 J=1,NS                                                       SR3  96 
C     *** INSERT CALL TO FAST SORT SUBROUTINE HERE THAT TAKES ADVANTAGE SR3  97 
C     OF PARTIAL ORDERING.  SEE G-L(SSA-II).                            SR3  98 
      CALL SORT (ISIM,NV,1,PD(1,J),INDEXI(1,J))                         SR3  99 
      S=0.                                                              SR3 100 
      DO 8 I=1,NV                                                       SR3 101 
      S=S+1.                                                            SR3 102 
      KK=INDEXI(I,J)+NS                                                 SR3 103 
8     DIST(J,KK)=S                                                      SR3 104 
9     CONTINUE                                                          SR3 105 
      LL=NS                                                             SR3 106 
C     ORDER COLUMNS OF INPUT MATRIX                                     SR3 107 
      DO 11 J=1,NV                                                      SR3 108 
C     *** SEE NOTES RE SORT SUBROUTINE.                                 SR3 109 
      CALL SORT (ISIM,NS,1,DP(1,J),INDI(1,J))                           SR3 110 
      S=0.                                                              SR3 111 
      LL=LL+1                                                           SR3 112 
      DO 10 I=1,NS                                                      SR3 113 
      S=S+1.                                                            SR3 114 
      KK=INDI(I,J)                                                      SR3 115 
10    DIST(LL,KK)=S                                                     SR3 116 
11    CONTINUE                                                          SR3 117 
C     PREPARE MATRIX FOR OBTAINING INITIAL CONFIGURATION                SR3 118 
      STR=2.*FN                                                         SR3 119 
      A=(2.*FNV+FNS)/STR                                                SR3 120 
      DO 12 I=1,NSM1                                                    SR3 121 
      DIST(I,I)=1.                                                      SR3 122 
      IP1=I+1                                                           SR3 123 
      DO 12 J=IP1,NS                                                    SR3 124 
      DIST(I,J)=A                                                       SR3 125 
12    DIST(J,I)=A                                                       SR3 126 
      DIST(NS,NS)=1.                                                    SR3 127 
      A=(2.*FNS+FNV)/STR                                                SR3 128 
      DO 13 I=NSP1,NSNVM1                                               SR3 129 
      DIST(I,I)=1.                                                      SR3 130 
      IP1=I+1                                                           SR3 131 
      DO 13 J=IP1,NSPNV                                                 SR3 132 
      DIST(I,J)=A                                                       SR3 133 
13    DIST(J,I)=A                                                       SR3 134 
      DIST(NSPNV,NSPNV)=1.                                              SR3 135 
C     CHECK FOR TIES IN ROWS                                            SR3 136 
      NOTR=0                                                            SR3 137 
      DO 21 I=1,NS                                                      SR3 138 
      NTR(I)=0                                                          SR3 139 
      JJ=1                                                              SR3 140 
      II=-1                                                             SR3 141 
      KK=JJ                                                             SR3 142 
14    LL=1                                                              SR3 143 
      MM=INDEXI(JJ,I)                                                   SR3 144 
      MMM=MM+NS                                                         SR3 145 
      S=DIST(I,MMM)                                                     SR3 146 
15    KK=KK+1                                                           SR3 147 
      NN=INDEXI(KK,I)                                                   SR3 148 
      IF (PD(MM,I)-PD(NN,I)) 17,16,17                                   SR3 149 
16    LL=LL+1                                                           SR3 150 
      NOTR=NOTR+1                                                       SR3 151 
      NNN=NN+NS                                                         SR3 152 
      S=S+DIST(I,NNN)                                                   SR3 153 
      IF (KK-NV) 15,18,18                                               SR3 154 
17    IF (LL-2) 20,18,18                                                SR3 155 
18    II=II+2                                                           SR3 156 
      NTR(I)=NTR(I)+2                                                   SR3 157 
      NOTYR(I,II)=LL                                                    SR3 158 
      NOTYR(I,II+1)=JJ                                                  SR3 159 
      S=S/FLOAT(LL)                                                     SR3 160 
      IJ=JJ                                                             SR3 161 
      DO 19 KL=1,LL                                                     SR3 162 
      MMM=INDEXI(IJ,I)+NS                                               SR3 163 
      DIST(I,MMM)=S                                                     SR3 164 
19    IJ=IJ+1                                                           SR3 165 
20    JJ=KK                                                             SR3 166 
      IF (JJ-NV) 14,21,21                                               SR3 167 
21    CONTINUE                                                          SR3 168 
C     CHECK FOR TIES IN COLUMNS                                         SR3 169 
      NOTC=0                                                            SR3 170 
      MMM=NS                                                            SR3 171 
      DO 29 I=1,NV                                                      SR3 172 
      MMM=MMM+1                                                         SR3 173 
      NTC(I)=0                                                          SR3 174 
      JJ=1                                                              SR3 175 
      II=-1                                                             SR3 176 
      KK=JJ                                                             SR3 177 
22    LL=1                                                              SR3 178 
      MM=INDI(JJ,I)                                                     SR3 179 
      S=DIST(MMM,MM)                                                    SR3 180 
23    KK=KK+1                                                           SR3 181 
      NN=INDI(KK,I)                                                     SR3 182 
      IF (DP(MM,I)-DP(NN,I)) 25,24,25                                   SR3 183 
24    LL=LL+1                                                           SR3 184 
      NOTC=NOTC+1                                                       SR3 185 
      S=S+DIST(MMM,NN)                                                  SR3 186 
      IF (KK-NS) 23,26,26                                               SR3 187 
25    IF (LL-2) 28,26,26                                                SR3 188 
26    II=II+2                                                           SR3 189 
      NTC(I)=NTC(I)+2                                                   SR3 190 
      NOTYC(I,II)=LL                                                    SR3 191 
      NOTYC(I,II+1)=JJ                                                  SR3 192 
      S=S/FLOAT(LL)                                                     SR3 193 
      IJ=JJ                                                             SR3 194 
      DO 27 KL=1,LL                                                     SR3 195 
      MM=INDI(IJ,I)                                                     SR3 196 
      DIST(MMM,MM)=S                                                    SR3 197 
27    IJ=IJ+1                                                           SR3 198 
28    JJ=KK                                                             SR3 199 
      IF (JJ-NS) 22,29,29                                               SR3 200 
29    CONTINUE                                                          SR3 201 
C     SYMMETRIZE                                                        SR3 202 
      DO 30 I=1,NS                                                      SR3 203 
      KK=NS                                                             SR3 204 
      DO 30 J=1,NV                                                      SR3 205 
      KK=KK+1                                                           SR3 206 
      DIST(KK,I)=(DIST(KK,I)+DIST(I,KK))/STR                            SR3 207 
30    DIST(I,KK)=DIST(KK,I)                                             SR3 208 
C     OBTAIN C-MATRIX                                                   SR3 209 
      DO 32 I=1,NSPNV                                                   SR3 210 
      DO 32 J=1,NSPNV                                                   SR3 211 
      IF (I-J) 31,32,31                                                 SR3 212 
31    DIST(I,I)=DIST(I,I)+DIST(I,J)                                     SR3 213 
32    CONTINUE                                                          SR3 214 
C     COMPLEMENT C                                                      SR3 215 
      DO 33 I=1,NSNVM1                                                  SR3 216 
      IP1=I+1                                                           SR3 217 
      DO 33 J=IP1,NSPNV                                                 SR3 218 
      DIST(I,J)=1.-DIST(I,J)                                            SR3 219 
33    DIST(J,I)=DIST(I,J)                                               SR3 220 
C     OBTAIN INITIAL CONFIGURATION                                      SR3 221 
      KK=MAXD+1                                                         SR3 222 
C     CALL ON HOUSEHOLDER SUBROUTINE                                    SR3 223 
      CALL EIGEN (DIST,C,NSPNV,STRESS,KK,MD,X(1,1),X(1,4),X(1,5),X(1,6),SR3 224 
     1X(1,7),X(1,8))                                                    SR3 225 
      DO 34 J=2,KK                                                      SR3 226 
      STRESS(J-1)=STRESS(J)                                             SR3 227 
      DO 34 I=1,NSPNV                                                   SR3 228 
34    X(I,J-1)=C(I,J)                                                   SR3 229 
C     INITIALIZE                                                        SR3 230 
      ITER=100                                                          SR3 231 
      PRESTR=1.                                                         SR3 232 
      A=1.                                                              SR3 233 
      ASSIGN 55 TO N1                                                   SR3 234 
      LFACT=1                                                           SR3 235 
C     DETERMINE WHETHER TO GO UP OR DOWN                                SR3 236 
      IF (MIND) 35,35,42                                                SR3 237 
35    DIM=NSNVM1*2                                                      SR3 238 
      DIM=FNR**2/DIM                                                    SR3 239 
      II=0                                                              SR3 240 
      DO 37 J=1,MAXD                                                    SR3 241 
      IF (STRESS(J)-DIM) 38,36,36                                       SR3 242 
36    II=II+1                                                           SR3 243 
37    CONTINUE                                                          SR3 244 
38    IF (II) 39,39,40                                                  SR3 245 
39    MAXD=1                                                            SR3 246 
      GO TO 41                                                          SR3 247 
40    MAXD=II                                                           SR3 248 
41    MIND=MAXD                                                         SR3 249 
      LFACT=0                                                           SR3 250 
      ITER=25                                                           SR3 251 
42    M=MIND-1                                                          SR3 252 
C     BEGIN LEAST-SQUARES-RANK-IMAGES ITERATIONS                        SR3 253 
43    II=0                                                              SR3 254 
      JJ=0                                                              SR3 255 
      M=M+1                                                             SR3 256 
      NTIER=NOTR                                                        SR3 257 
      NTIEC=NOTC                                                        SR3 258 
      NCY=1                                                             SR3 259 
44    II=II+1                                                           SR3 260 
45    ITCT=0                                                            SR3 261 
C     COMPUTE DISTANCES                                                 SR3 262 
46    IF (IFCB) 47,50,47                                                SR3 263 
47    DO 49 I=1,NS                                                      SR3 264 
      KK=0                                                              SR3 265 
      DO 49 J=NSP1,NSPNV                                                SR3 266 
      KK=KK+1                                                           SR3 267 
      DIST(I,KK)=0.                                                     SR3 268 
      DO 48 K=1,M                                                       SR3 269 
48    DIST(I,KK)=DIST(I,KK)+ABS(X(I,K)-X(J,K))                          SR3 270 
49    DIST(J,I)=DIST(I,KK)                                              SR3 271 
      GO TO 53                                                          SR3 272 
50    DO 52 I=1,NS                                                      SR3 273 
      KK=0                                                              SR3 274 
      DO 52 J=NSP1,NSPNV                                                SR3 275 
      KK=KK+1                                                           SR3 276 
      DIST(I,KK)=0.                                                     SR3 277 
      DO 51 K=1,M                                                       SR3 278 
51    DIST(I,KK)=DIST(I,KK)+(X(I,K)-X(J,K))**2                          SR3 279 
C     COMPUTE DISTANCE                                                  SR3 280 
      DIST(I,KK)=SQRT(DIST(I,KK))                                       SR3 281 
52    DIST(J,I)=DIST(I,KK)                                              SR3 282 
53    IF (ITCT) 120,54,120                                              SR3 283 
54    GO TO N1, (55,58)                                                 SR3 284 
55    DO 56 J=1,NV                                                      SR3 285 
C     SORT DISTANCES WITHIN ROWS AND WITHIN COLUMNS                     SR3 286 
56    CALL SORT (0,NS,1,DIST(1,J),INDJ(1,J))                            SR3 287 
      DO 57 J=1,NS                                                      SR3 288 
C     *** CALL ON SORT SUBROUTINE.                                      SR3 289 
57    CALL SORT (0,NV,1,DIST(NSP1,J),INDEXJ(1,J))                       SR3 290 
      ASSIGN 58 TO N1                                                   SR3 291 
      GO TO 61                                                          SR3 292 
58    DO 59 J=1,NV                                                      SR3 293 
C     CALL ON SORT                                                      SR3 294 
59    CALL SORT (0,NS,0,DIST(1,J),INDJ(1,J))                            SR3 295 
      DO 60 J=1,NS                                                      SR3 296 
C     *** CALL ON SORT.                                                 SR3 297 
60    CALL SORT (0,NV,0,DIST(NSP1,J),INDEXJ(1,J))                       SR3 298 
61    IF (NTIER) 72,72,62                                               SR3 299 
C     OPTIMALLY PERMUTE TIES ON FIRST AND LAST ITERATION FOR ROWS       SR3 300 
62    NTIER=0                                                           SR3 301 
      DO 71 I=1,NS                                                      SR3 302 
      IF (NTR(I)) 63,71,63                                              SR3 303 
63    KK=NTR(I)                                                         SR3 304 
      DO 70 J=1,KK,2                                                    SR3 305 
      LL=NOTYR(I,J)                                                     SR3 306 
      MM=NOTYR(I,J+1)                                                   SR3 307 
      NN=MM                                                             SR3 308 
      DO 64 K=1,LL                                                      SR3 309 
      IJ=INDEXI(NN,I)                                                   SR3 310 
      EVAL(K)=DIST(I,IJ)                                                SR3 311 
      IND(K,1)=IJ                                                       SR3 312 
64    NN=NN+1                                                           SR3 313 
      IJ=LL-1                                                           SR3 314 
65    JI=0                                                              SR3 315 
      DO 67 L=1,IJ                                                      SR3 316 
      LP1=L+1                                                           SR3 317 
      DO 67 N=LP1,LL                                                    SR3 318 
      IF (EVAL(N)-EVAL(L)) 66,67,67                                     SR3 319 
66    S=EVAL(L)                                                         SR3 320 
      EVAL(L)=EVAL(N)                                                   SR3 321 
      EVAL(N)=S                                                         SR3 322 
      JI=1                                                              SR3 323 
      NN=IND(L,1)                                                       SR3 324 
      IND(L,1)=IND(N,1)                                                 SR3 325 
      IND(N,1)=NN                                                       SR3 326 
67    CONTINUE                                                          SR3 327 
      IF (JI) 68,68,65                                                  SR3 328 
68    NN=MM                                                             SR3 329 
      DO 69 KL=1,LL                                                     SR3 330 
      INDEXI(NN,I)=IND(KL,1)                                            SR3 331 
69    NN=NN+1                                                           SR3 332 
70    CONTINUE                                                          SR3 333 
71    CONTINUE                                                          SR3 334 
72    IF (NTIEC) 83,83,73                                               SR3 335 
C     OPTIMALLY PERMUTE TIES FOR FIRST AND LAST ITERATION FOR COLUMNS   SR3 336 
73    NTIEC=0                                                           SR3 337 
      MMM=NS                                                            SR3 338 
      DO 82 I=1,NV                                                      SR3 339 
      MMM=MMM+1                                                         SR3 340 
      IF (NTC(I)) 74,82,74                                              SR3 341 
74    KK=NTC(I)                                                         SR3 342 
      DO 81 J=1,KK,2                                                    SR3 343 
      LL=NOTYC(I,J)                                                     SR3 344 
      MM=NOTYC(I,J+1)                                                   SR3 345 
      NN=MM                                                             SR3 346 
      DO 75 K=1,LL                                                      SR3 347 
      IJ=INDI(NN,I)                                                     SR3 348 
      EVAL(K)=DIST(MMM,IJ)                                              SR3 349 
      IND(K,1)=IJ                                                       SR3 350 
75    NN=NN+1                                                           SR3 351 
      IJ=LL-1                                                           SR3 352 
76    JI=0                                                              SR3 353 
      DO 78 L=1,IJ                                                      SR3 354 
      LP1=L+1                                                           SR3 355 
      DO 78 N=LP1,LL                                                    SR3 356 
      IF (EVAL(N)-EVAL(L)) 77,78,78                                     SR3 357 
77    S=EVAL(L)                                                         SR3 358 
      EVAL(L)=EVAL(N)                                                   SR3 359 
      EVAL(N)=S                                                         SR3 360 
      JI=1                                                              SR3 361 
      NN=IND(L,1)                                                       SR3 362 
      IND(L,1)=IND(N,1)                                                 SR3 363 
      IND(N,1)=NN                                                       SR3 364 
78    CONTINUE                                                          SR3 365 
      IF (JI) 79,79,76                                                  SR3 366 
79    NN=MM                                                             SR3 367 
      DO 80 KL=1,LL                                                     SR3 368 
      INDI(NN,I)=IND(KL,1)                                              SR3 369 
80    NN=NN+1                                                           SR3 370 
81    CONTINUE                                                          SR3 371 
82    CONTINUE                                                          SR3 372 
C     PERFORM RANK-IMAGE PERMUTATIONS FOR ROWS                          SR3 373 
83    DO 84 J=1,NS                                                      SR3 374 
      DO 84 I=1,NV                                                      SR3 375 
      KK=INDEXI(I,J)                                                    SR3 376 
      LL=INDEXJ(I,J)                                                    SR3 377 
84    PD(KK,J)=DIST(J,LL)                                               SR3 378 
C     RANK-IMAGE PERMUTATIONS FOR COLUMNS                               SR3 379 
      DO 85 J=1,NV                                                      SR3 380 
      DO 85 I=1,NS                                                      SR3 381 
      KK=INDI(I,J)                                                      SR3 382 
      LL=INDJ(I,J)                                                      SR3 383 
85    DP(KK,J)=DIST(LL,J)                                               SR3 384 
      GO TO (86,88), MU                                                 SR3 385 
C     SYMMETRIZE                                                        SR3 386 
86    DO 87 I=1,NS                                                      SR3 387 
      DO 87 J=1,NV                                                      SR3 388 
      DP(I,J)=(DP(I,J)+PD(J,I))/2.                                      SR3 389 
87    PD(J,I)=DP(I,J)                                                   SR3 390 
88    IF (JJ) 148,89,148                                                SR3 391 
C     TIE ROW TIES                                                      SR3 392 
89    IF (NOTR) 96,96,90                                                SR3 393 
90    DO 95 I=1,NS                                                      SR3 394 
      IF (NTR(I)) 91,95,91                                              SR3 395 
91    KK=NTR(I)                                                         SR3 396 
      DO 94 J=1,KK,2                                                    SR3 397 
      LL=NOTYR(I,J)                                                     SR3 398 
      MM=NOTYR(I,J+1)                                                   SR3 399 
      NN=MM                                                             SR3 400 
      S=0.                                                              SR3 401 
      DO 92 K=1,LL                                                      SR3 402 
      IJ=INDEXI(NN,I)                                                   SR3 403 
      S=S+PD(IJ,I)                                                      SR3 404 
92    NN=NN+1                                                           SR3 405 
      NN=MM                                                             SR3 406 
      S=S/FLOAT(LL)                                                     SR3 407 
      DO 93 L=1,LL                                                      SR3 408 
      IJ=INDEXI(NN,I)                                                   SR3 409 
      PD(IJ,I)=S                                                        SR3 410 
93    NN=NN+1                                                           SR3 411 
94    CONTINUE                                                          SR3 412 
95    CONTINUE                                                          SR3 413 
96    IF (NOTC) 103,103,97                                              SR3 414 
C     TIE COLUMN TIES                                                   SR3 415 
97    DO 102 I=1,NV                                                     SR3 416 
      IF (NTC(I)) 98,102,98                                             SR3 417 
98    KK=NTC(I)                                                         SR3 418 
      DO 101 J=1,KK,2                                                   SR3 419 
      LL=NOTYC(I,J)                                                     SR3 420 
      MM=NOTYC(I,J+1)                                                   SR3 421 
      NN=MM                                                             SR3 422 
      S=0.                                                              SR3 423 
      DO 99 K=1,LL                                                      SR3 424 
      IJ=INDI(NN,I)                                                     SR3 425 
      S=S+DP(IJ,I)                                                      SR3 426 
99    NN=NN+1                                                           SR3 427 
      NN=MM                                                             SR3 428 
      S=S/FLOAT(LL)                                                     SR3 429 
      DO 100 L=1,LL                                                     SR3 430 
      IJ=INDI(NN,I)                                                     SR3 431 
      DP(IJ,I)=S                                                        SR3 432 
100   NN=NN+1                                                           SR3 433 
101   CONTINUE                                                          SR3 434 
102   CONTINUE                                                          SR3 435 
C     TEST FOR ZERO DISTANCE                                            SR3 436 
103   DO 104 J=1,NS                                                     SR3 437 
      KK=INDEXJ(1,J)                                                    SR3 438 
      IF (DIST(J,KK)) 104,105,104                                       SR3 439 
104   CONTINUE                                                          SR3 440 
      GO TO 108                                                         SR3 441 
105   DO 106 I=2,NV                                                     SR3 442 
      KK=INDEXJ(I,J)                                                    SR3 443 
      IF (DIST(J,KK)) 107,106,107                                       SR3 444 
106   CONTINUE                                                          SR3 445 
107   LL=INDEXJ(NV,J)                                                   SR3 446 
      GL=DIST(J,LL)/DIST(J,KK)                                          SR3 447 
      GO TO 109                                                         SR3 448 
108   GL=0.                                                             SR3 449 
109   GO TO (117,110), MU                                               SR3 450 
110   DO 116 I=1,NS                                                     SR3 451 
      MM=NS                                                             SR3 452 
      DO 116 J=1,NV                                                     SR3 453 
      MM=MM+1                                                           SR3 454 
      IF (GL) 111,115,111                                               SR3 455 
111   IF (DP(I,J)+DIST(I,J)) 113,112,113                                SR3 456 
112   C(MM,I)=1.                                                        SR3 457 
      GO TO 116                                                         SR3 458 
113   IF (DIST(I,J)) 115,114,115                                        SR3 459 
114   C(MM,I)=GL                                                        SR3 460 
      GO TO 116                                                         SR3 461 
115   C(MM,I)=DP(I,J)/DIST(I,J)                                         SR3 462 
116   CONTINUE                                                          SR3 463 
      NCY=0                                                             SR3 464 
      GO TO 127                                                         SR3 465 
117   IF (NOTR+NOTC) 120,120,118                                        SR3 466 
C     SYMMETRIZE                                                        SR3 467 
118   DO 119 I=1,NS                                                     SR3 468 
      DO 119 J=1,NV                                                     SR3 469 
119   PD(J,I)=(DP(I,J)+PD(J,I))/2.                                      SR3 470 
120   DO 126 I=1,NS                                                     SR3 471 
      MM=NS                                                             SR3 472 
      DO 126 J=1,NV                                                     SR3 473 
      MM=MM+1                                                           SR3 474 
      IF (GL) 121,125,121                                               SR3 475 
121   IF (PD(J,I)+DIST(I,J)) 123,122,123                                SR3 476 
122   C(MM,I)=1.                                                        SR3 477 
      GO TO 126                                                         SR3 478 
123   IF (DIST(I,J)) 125,124,125                                        SR3 479 
124   C(MM,I)=GL                                                        SR3 480 
      GO TO 126                                                         SR3 481 
125   C(MM,I)=PD(J,I)/DIST(I,J)                                         SR3 482 
126   CONTINUE                                                          SR3 483 
C     COMPUTE CORRECTION MATRIX                                         SR3 484 
127   DO 128 I=NSP1,NSPNV                                               SR3 485 
      C(I,I)=FNV                                                        SR3 486 
      DO 128 J=1,NS                                                     SR3 487 
128   C(I,I)=C(I,I)+C(I,J)                                              SR3 488 
      DO 129 J=1,NS                                                     SR3 489 
      C(J,J)=FNS                                                        SR3 490 
      DO 129 I=NSP1,NSPNV                                               SR3 491 
129   C(J,J)=C(J,J)+C(I,J)                                              SR3 492 
      DO 130 I=NSP1,NSPNV                                               SR3 493 
      DO 130 J=1,NS                                                     SR3 494 
130   C(I,J)=1.-C(I,J)                                                  SR3 495 
C     APPLY CORRECTIONS                                                 SR3 496 
      DO 132 K=1,M                                                      SR3 497 
      DO 132 I=1,NS                                                     SR3 498 
      XP1(I,K)=X(I,K)*C(I,I)                                            SR3 499 
      DO 131 J=NSP1,NSPNV                                               SR3 500 
131   XP1(I,K)=XP1(I,K)+X(J,K)*C(J,I)                                   SR3 501 
132   XP1(I,K)=XP1(I,K)/FNR                                             SR3 502 
      DO 134 K=1,M                                                      SR3 503 
      DO 134 I=NSP1,NSPNV                                               SR3 504 
      XP1(I,K)=X(I,K)*C(I,I)                                            SR3 505 
      DO 133 J=1,NS                                                     SR3 506 
133   XP1(I,K)=XP1(I,K)+X(J,K)*C(I,J)                                   SR3 507 
134   XP1(I,K)=XP1(I,K)/FNR                                             SR3 508 
      DO 135 K=1,M                                                      SR3 509 
      DO 135 I=1,NSPNV                                                  SR3 510 
135   X(I,K)=XP1(I,K)                                                   SR3 511 
      IF (NCY) 137,136,137                                              SR3 512 
136   NCY=1                                                             SR3 513 
      GO TO 120                                                         SR3 514 
137   ITCT=ITCT+1                                                       SR3 515 
      IF (IT-ITCT) 138,138,46                                           SR3 516 
C     COMPUTE NORMALIZED PHI                                            SR3 517 
138   S=0.                                                              SR3 518 
      STR=0.                                                            SR3 519 
      DO 139 I=1,NS                                                     SR3 520 
      DO 139 J=1,NV                                                     SR3 521 
      S=S+DIST(I,J)**2                                                  SR3 522 
139   STR=STR+(DIST(I,J)-PD(J,I))**2                                    SR3 523 
      STR=STR/(2.*S)                                                    SR3 524 
      STRESS(II)=STR                                                    SR3 525 
C     TEST FOR TERMINATION                                              SR3 526 
      IF (II-6) 44,140,140                                              SR3 527 
140   IF (STRESS(II-5)-STRESS(II)-1E-06) 142,142,141                    SR3 528 
141   IF (II-ITER) 44,142,142                                           SR3 529 
142   IF (LFACT) 145,145,143                                            SR3 530 
143   IF (PRESTR-STR) 1,1,144                                           SR3 531 
144   PRESTR=STR                                                        SR3 532 
C     PRINT RESULTS                                                     SR3 533 
145   WRITE (JTAPE,196)                                                 SR3 534 
      WRITE (JTAPE,202) M,(KK,KK=1,M)                                   SR3 535 
      WRITE (JTAPE,203)                                                 SR3 536 
      JJ=1                                                              SR3 537 
      NTIER=NOTR                                                        SR3 538 
      GO TO (147,146), MU                                               SR3 539 
146   NTIEC=0                                                           SR3 540 
      GO TO 45                                                          SR3 541 
C     COMPUTE INDIVIDUALIZED PHI'S                                      SR3 542 
147   NTIEC=NOTC                                                        SR3 543 
      GO TO 45                                                          SR3 544 
148   S=0.                                                              SR3 545 
      DO 149 I=1,NS                                                     SR3 546 
      PHI(I)=0.                                                         SR3 547 
      DO 149 J=1,NV                                                     SR3 548 
      S=S+DIST(I,J)**2                                                  SR3 549 
149   PHI(I)=PHI(I)+(DIST(I,J)-PD(J,I))**2                              SR3 550 
      STR=0.                                                            SR3 551 
      GO TO (153,150), MU                                               SR3 552 
150   S=2.*S                                                            SR3 553 
      DO 151 J=1,NS                                                     SR3 554 
      PHI(J)=PHI(J)/S                                                   SR3 555 
151   STR=STR+PHI(J)                                                    SR3 556 
      DO 152 J=NSP1,NSPNV                                               SR3 557 
152   PHI(J)=0.                                                         SR3 558 
      GO TO 156                                                         SR3 559 
153   KK=NS                                                             SR3 560 
      DO 154 J=1,NV                                                     SR3 561 
      KK=KK+1                                                           SR3 562 
      PHI(KK)=0.                                                        SR3 563 
      DO 154 I=1,NS                                                     SR3 564 
154   PHI(KK)=PHI(KK)+(DIST(I,J)-DP(I,J))**2                            SR3 565 
      S=4.*S                                                            SR3 566 
      DO 155 I=1,NSPNV                                                  SR3 567 
      PHI(I)=PHI(I)/S                                                   SR3 568 
155   STR=STR+PHI(I)                                                    SR3 569 
156   IF (M-2) 162,157,157                                              SR3 570 
157   IF (IFCB) 162,158,162                                             SR3 571 
C     ROTATE TO PRINCIPAL AXIS                                          SR3 572 
158   DO 160 I=1,NSPNV                                                  SR3 573 
      DO 160 J=I,NSPNV                                                  SR3 574 
      DIST(I,J)=0.                                                      SR3 575 
      DO 159 K=1,M                                                      SR3 576 
159   DIST(I,J)=DIST(I,J)+X(I,K)*X(J,K)                                 SR3 577 
160   DIST(J,I)=DIST(I,J)                                               SR3 578 
C     CALL ON HOUSEHOLDER SUBROUTINE                                    SR3 579 
      CALL EIGEN (DIST,C,NSPNV,STRESS,M,MD,XP1(1,1),XP1(1,4),XP1(1,5),XPSR3 580 
     11(1,6),XP1(1,7),XP1(1,8))                                         SR3 581 
      DO 161 J=1,M                                                      SR3 582 
      DO 161 I=1,NSPNV                                                  SR3 583 
161   X(I,J)=C(I,J)                                                     SR3 584 
C     NORMALIZE COORDINATES TO LIE IN RANGE -1 TO +1                    SR3 585 
162   A=0.                                                              SR3 586 
      DO 164 K=1,M                                                      SR3 587 
      C(1,K)=X(1,K)                                                     SR3 588 
      C(2,K)=X(1,K)                                                     SR3 589 
      DO 163 I=2,NSPNV                                                  SR3 590 
      C(1,K)=AMIN1(C(1,K),X(I,K))                                       SR3 591 
163   C(2,K)=AMAX1(C(2,K),X(I,K))                                       SR3 592 
164   A=AMAX1(A,(C(2,K)-C(1,K)))                                        SR3 593 
      A=2./A                                                            SR3 594 
      DO 165 K=1,M                                                      SR3 595 
      DO 165 I=1,NSPNV                                                  SR3 596 
165   XP1(I,K)=A*(X(I,K)-C(1,K))-1.                                     SR3 597 
C     COMPUTE CENTRALITY INDEX                                          SR3 598 
      KK=M+1                                                            SR3 599 
      DO 169 I=1,NSPNV                                                  SR3 600 
      XP1(I,KK)=0.                                                      SR3 601 
      DO 166 K=1,M                                                      SR3 602 
166   XP1(I,KK)=XP1(I,KK)+X(I,K)**2                                     SR3 603 
      IF (IFCB) 167,168,167                                             SR3 604 
167   XP1(I,KK)=A*XP1(I,KK)                                             SR3 605 
      GO TO 169                                                         SR3 606 
168   XP1(I,KK)=A*SQRT(XP1(I,KK))                                       SR3 607 
169   CONTINUE                                                          SR3 608 
      DO 170 I=1,NS                                                     SR3 609 
170   WRITE (JTAPE,204) I,PHI(I),(XP1(I,K),K=1,KK)                      SR3 610 
      WRITE (JTAPE,206)                                                 SR3 611 
      DO 171 I=NSP1,NSPNV                                               SR3 612 
171   WRITE (JTAPE,204) I,PHI(I),(XP1(I,K),K=1,KK)                      SR3 613 
      CAL=SQRT(1.-(1.-STR)**2)                                          SR3 614 
      WRITE (JTAPE,205) STR,II,CAL                                      SR3 615 
      IF (LFACT) 172,184,172                                            SR3 616 
172   IF (M-2) 194,173,173                                              SR3 617 
C     CALL ON PLOT SUBROUTINE                                           SR3 618 
173   CALL PLOT (XP1,DIST(1,1),DIST(1,2),DIST(1,3),DIST(1,4),DIST(1,5),NSR3 619 
     1SPNV,M,MD,ND)                                                     SR3 620 
      IF (IFC) 174,176,174                                              SR3 621 
174   DO 175 I=1,NSPNV                                                  SR3 622 
175   PUNCH 207, (XP1(I,J),J=1,M)                                       SR3 623 
176   IF (IFD) 177,194,177                                              SR3 624 
177   IF (IFCB) 181,178,181                                             SR3 625 
178   DO 180 I=1,NSNVM1                                                 SR3 626 
      DIST(I,I)=0.                                                      SR3 627 
      IP1=I+1                                                           SR3 628 
      DO 180 J=IP1,NSPNV                                                SR3 629 
      DIST(J,I)=0.                                                      SR3 630 
      DO 179 K=1,M                                                      SR3 631 
179   DIST(J,I)=DIST(J,I)+(XP1(I,K)-XP1(J,K))**2                        SR3 632 
180   DIST(J,I)=SQRT(DIST(J,I))                                         SR3 633 
      GO TO 183                                                         SR3 634 
181   DO 182 I=1,NSNVM1                                                 SR3 635 
      DIST(I,I)=0.                                                      SR3 636 
      IP1=I+1                                                           SR3 637 
      DO 182 J=IP1,NSPNV                                                SR3 638 
      DIST(J,I)=0.                                                      SR3 639 
      DO 182 K=1,M                                                      SR3 640 
182   DIST(J,I)=DIST(J,I)+ABS(XP1(I,K)-XP1(J,K))                        SR3 641 
183   DIST(NSPNV,NSPNV)=0.                                              SR3 642 
      CALL MXOUT (DIST,NSPNV,1,MD)                                      SR3 643 
      GO TO 194                                                         SR3 644 
184   IF (IFCB) 185,186,185                                             SR3 645 
185   WRITE (JTAPE,208)                                                 SR3 646 
      GO TO 1                                                           SR3 647 
C     DETERMINE HOW MANY DIMENSIONS TO DROP                             SR3 648 
186   LL=1                                                              SR3 649 
      KK=M                                                              SR3 650 
      DO 188 J=2,KK                                                     SR3 651 
      S=STRESS(J)/STRESS(1)                                             SR3 652 
      IF (S-.3) 189,189,187                                             SR3 653 
187   LL=LL+1                                                           SR3 654 
188   CONTINUE                                                          SR3 655 
189   M=LL-1                                                            SR3 656 
      MAXD=LL                                                           SR3 657 
      IF (ITER-100) 190,191,191                                         SR3 658 
190   ITER=ITER+25                                                      SR3 659 
191   IF (KK-LL) 192,192,43                                             SR3 660 
192   LFACT=1                                                           SR3 661 
      IF (M) 1,1,193                                                    SR3 662 
193   M=M+1                                                             SR3 663 
      GO TO 173                                                         SR3 664 
194   IF (.000001-STR) 195,1,1                                          SR3 665 
195   IF (M-MAXD) 43,1,1                                                SR3 666 
C     *** FORMAT STATEMENTS ***                                         SR3 667 
C                                                                       SR3 668 
196   FORMAT (72H                                                       SR3 669 
     1                 )                                                SR3 670 
197   FORMAT (9I4)                                                      SR3 671 
198   FORMAT (18A4)                                                     SR3 672 
199   FORMAT (19H0INPUT COEFFICIENTS)                                   SR3 673 
200   FORMAT (4H0ROW,I3)                                                SR3 674 
201   FORMAT (10F13.4)                                                  SR3 675 
202   FORMAT (34H0G-L(SSAR-III) COORDINATES FOR M =,I3,1H./1H0,9X10H DIMSR3 676 
     1ENSION,10I10)                                                     SR3 677 
203   FORMAT (1H ,130(1H-)/6H0  ROW,8X6HN. PHI)                         SR3 678 
204   FORMAT (I5,5XF10.5,2P11F10.3)                                     SR3 679 
205   FORMAT (17H0NORMALIZED PHI =,F7.5,4H FOR,I4,11H ITERATIONS/28H COESR3 680 
     1FFICIENT OF ALIENATION =,E12.6)                                   SR3 681 
206   FORMAT (8H0 COLUMN)                                               SR3 682 
207   FORMAT (9F8.5)                                                    SR3 683 
208   FORMAT (27H0ERROR - IFCB=1 AND MIND=0.)                           SR3 684 
      END                                                               SR3 685-
C     SORT/ORDER                                                        SRT   1 
C     THE FOLLOWING SUBROUTINE CAN BE CHANGED TO SORT OR ORDER BY RE-   SRT   2 
C     MOVING THE APPROPRIATE CARD NAMING THE SUBROUTINE.  THIS FORTRAN  SRT   3 
C     IV SUBROUTINE (A N LOG(BASE 2) N) SORT WILL ORDER A SET OF FLOAT- SRT   4 
C     ING POINT NUMBERS FROM LOW TO HIGH (ISIM=0) OR FROM HIGH TO LOW   SRT   5 
C     (SIM=1) WITHOUT DISTURBING THE NUMBERS TO BE SORTED.  THE DIMEN-  SRT   6 
C     SION OF LIST SHOULD LIE BETWEEN 2 LOG (BASE 2) N AND N.   PRO-    SRT   7 
C     GRAMMED BY ARDEN/LINGOES, U. OF M. COMPUTING CENTER (10/18/67).   SRT   8 
C                                                                       SRT   9 
      SUBROUTINE SORT (ISIM,NEL,ISW,DIST,INDEX)                         SRT  10 
      DIMENSION DIST(1), INDEX(1), LIST(200)                            SRT  11 
C                                                                       SRT  12 
1     N=NEL                                                             SRT  13 
      IF (ISW) 6,6,2                                                    SRT  14 
2     DO 3 J=1,N                                                        SRT  15 
3     INDEX(J)=J                                                        SRT  16 
      IF (ISIM) 5,4,5                                                   SRT  17 
4     ASSIGN 10 TO N1                                                   SRT  18 
      GO TO 6                                                           SRT  19 
5     ASSIGN 16 TO N1                                                   SRT  20 
6     I1=1                                                              SRT  21 
      IN=N                                                              SRT  22 
      LN=0                                                              SRT  23 
7     I=I1-1                                                            SRT  24 
      K=I                                                               SRT  25 
      DN=DIST(INDEX(IN))                                                SRT  26 
8     K=K+1                                                             SRT  27 
      DF=DIST(INDEX(K))                                                 SRT  28 
      IF (DF.NE.DN) GO TO 9                                             SRT  29 
      IF (K-IN) 8,18,18                                                 SRT  30 
9     J=IN+1                                                            SRT  31 
      AV=.5*ABS(DF-DN)+AMIN1(DF,DN)                                     SRT  32 
      IF (AV.EQ.AMAX1(DF,DN)) AV=AMIN1(DF,DN)                           SRT  33 
      GO TO N1, (10,16)                                                 SRT  34 
10    I=I+1                                                             SRT  35 
      IF (DIST(INDEX(I)).LE.AV) GO TO 10                                SRT  36 
11    J=J-1                                                             SRT  37 
      IF (DIST(INDEX(J)).GT.AV) GO TO 11                                SRT  38 
12    IF (I.GT.J) GO TO 13                                              SRT  39 
      ITEMP=INDEX(I)                                                    SRT  40 
      INDEX(I)=INDEX(J)                                                 SRT  41 
      INDEX(J)=ITEMP                                                    SRT  42 
      GO TO N1, (10,16)                                                 SRT  43 
13    IF (I.EQ.IN) GO TO 14                                             SRT  44 
      IF (J.EQ.I1) GO TO 15                                             SRT  45 
      LIST(LN+1)=J+1                                                    SRT  46 
      LIST(LN+2)=IN                                                     SRT  47 
      LN=LN+2                                                           SRT  48 
      IF (LN.EQ.200) GO TO 21                                           SRT  49 
14    IN=I-1                                                            SRT  50 
      GO TO 7                                                           SRT  51 
15    I1=J+1                                                            SRT  52 
      GO TO 7                                                           SRT  53 
16    I=I+1                                                             SRT  54 
      IF (AV.LT.DIST(INDEX(I))) GO TO 16                                SRT  55 
17    J=J-1                                                             SRT  56 
      IF (AV.LT.DIST(INDEX(J))) GO TO 12                                SRT  57 
      GO TO 17                                                          SRT  58 
18    IF (LN) 19,20,19                                                  SRT  59 
19    IN=LIST(LN)                                                       SRT  60 
      I1=LIST(LN-1)                                                     SRT  61 
      LN=LN-2                                                           SRT  62 
      GO TO 7                                                           SRT  63 
20    RETURN                                                            SRT  64 
21    CALL ERROR                                                        SRT  65 
      END                                                               SRT  66-
1UNIDIMENSIONAL CASE WITH COMPLETE DATA.                                        
  11   5   1   2       1   1                                                    
(36F2.0)                                                                        
 1 2 3 4 5                                                                      
 2 1 3 4 5                                                                      
 3 1 2 4 5                                                                      
 4 1 2 3 5                                                                      
 4 2 1 3 5                                                                      
 4 3 1 2 5                                                                      
 4 3 2 1 5                                                                      
 5 3 2 1 4                                                                      
 5 4 2 1 3                                                                      
 5 4 3 1 2                                                                      
 5 4 3 2 1                                                                      
1FOUR POINTS IN TWO DIMENSIONS - COOMB'S (P. 157).                              
  18   4   2   3       1                                                        
(4F2.0)                                                                         
 1 2 3 4                                                                        
 4 3 2 1                                                                        
 1 2 4 3                                                                        
 4 3 1 2                                                                        
 1 3 4 2                                                                        
 4 2 1 3                                                                        
 2 1 3 4                                                                        
 3 4 2 1                                                                        
 3 1 2 4                                                                        
 2 4 3 1                                                                        
 3 2 1 4                                                                        
 2 3 4 1                                                                        
 1 3 2 4                                                                        
 1 4 2 3                                                                        
 1 4 3 2                                                                        
 2 4 1 3                                                                        
 2 3 1 4                                                                        
 3 4 1 2                                                                        
1CHILDREN'S PREFERENCES FOR 4 RADIO STATIONS-TBL.7.,P.179,COOMB'S.              
  18   4       3       1                                                        
(36F2.0)                                                                        
 1 2 3 4                                                                        
 1 3 2 4                                                                        
 2 3 1 4                                                                        
 2 1 3 4                                                                        
 3 2 1 4                                                                        
 3 1 2 4                                                                        
 1 2 4 3                                                                        
 1 3 4 2                                                                        
 2 3 4 1                                                                        
 2 1 4 3                                                                        
 3 2 4 1                                                                        
 1 4 2 3                                                                        
 1 4 3 2                                                                        
 2 4 1 3                                                                        
 2 4 3 1                                                                        
 3 4 2 1                                                                        
 4 2 1 3                                                                        
 4 3 2 1                                                                        
