228 VI = STEP + STEPI * 1.D-10 



229 GO fo Ob 



230 62 IF (V) 142,301 ,70 



231 142 IF(STEP) 44,44,143 



232 C SEARCH FOR MODE 



233 143 SIZE3 = -1 . 



234 SIZE2=0 



235 IJ FLAG=1 



236 V=-V 



237 IF (NXTRA) 55.55,54 



238 55 NXTRA = 20 



239 54 XTRA = NXTRA 



240 HOP = (STEP - V) / XTRA 



241 HOPI=0. 



242 IF(STEPI.NE.O. ) HOP I = ( STEP I-VI ) /XTRA 



243 DO 47 IlI = 1 , NXTRA 



244 CALL SETUP 



245 DET = VEL 



246 DETI = VELI 



247 CALL DETNT(N, VEL, VELI) 



248 DELTA = VEL 



249 n^l TT - \/ci I 



250 SIZE = DELTA*DELTA + DELTI*DELTI 



251 PRINT 56, V, VI, SIZE, VEL, VELI 



252 56 FORMAT (2F12.3, 3D17.5) 



253 IF ( (SIZE2.lt. SIZES) .AND. (SIZE. GT.SIZE2)) GO TO 45 



254 46 SIZE3rSIZE2 



255 SIZE2=SIZE 



256 V = V + HOP 



257 VI=VI+HOPI 



258 GO TO 47 



259 45 V = V - HOP 



260 TEMP = HOP / (SIZE - SIZE2) 



261 DELTI = TEMP * (DET * VELI - DETI * VEL) 



262 TEMP = .500 * (SIZE3 - SIZE) / (SIZE3 + SIZE - SIZE2 - SIZE2) 



263 DELTA = HOP + TEMP 



264 IF(HOPI .EQ.O) GO TO 49 



265 VI=VI-HOPI 



266 DELTAI=HOPI*TEMP 



267 GO TO 49 



268 47 CONTINUE 



269 Kl FLAG=0 



270 NXTRA=0 



271 GO TO 52 



272 53 SIZE 2 = -1 . 



273 SIZE=0 



274 GO TO 46 



275 44 NXTRA = NXTRA - 1 



276 V = 3. * (PHASE V(NN-1) - PHASE V(NN-2)) + PHASE V(NN-3) 



277 VI = 3.* (PHASI V(NN-1) - PHASI V(NN-2)) + PHASI V(NN-3) 



278 STEP = (PHASE V(NN-1) - PHASE V ( NN-2 ) ) * .0001 



279 70 CALL SETUP 



280 CALL DETNT(N, DET, DETI) 



281 80 FORMAT (/, 2D20.11, 4D13.4) 



282 VEL = DET 



283 VELI = DETI 



284 DELTA = STEP 



48 



