285 DELTI = STEPI 



286 Ih (DhLIA .Nt.O.) GO TO 49 



287 IF (DELTI .EO.O.) DELTA = .01 



288 49 SIZE2 = 100. 



289 RX = DET**2 + DETI**2 



290 IF (K6 .LT. 3) PRINT 80, V, VI, DET , DETI, SIZE, CNTR 



291 J = 



292 48 J = J + 1 



293 IF (J .GT. 15) GO TO 51 



294 V = V + DELTA 



295 VI = VI + DELTI 



296 IF (VI) 83,84,85 



297 83 DELTI = DELTI - VI 



298 84 VI = 1 .D-18 



299 85 CALL SETUP 



300 NNN = N + N - 1 



301 DO 82 lA = 1 ,NNN 



302 DO 82 IB = 1,4 



303 BKIA.IB) = 0(IA,IB) 



304 82 B( lA, IB) = A( lA, IB) 



305 CALL DETNT(N, DET, DETI) 



306 IF ^K6 .NE. 1 ) GO TO 72 



307 71 PRINT 81, V, VI, DET, DETI, SIZE, CNTR 



308 81 FORMAT (2D20.11, 4D13.4) 



309 72 IF (NXTRA .LT. 0) GO TO 51 



310 TEMNR = DET ♦ DELTA - DETI * DELTI 



311 TEMNI = DETI ♦ DELTA + DET * DELTI 



312 TEMDR = VEL - DET 



313 TEMDI = VELI - DETI 



314 TEMDEN = TEMDR*TEMDR + TEMDI*TEMDI 



315 IF (TEMDEN . EQ . 0.) GO TO 51 



316 TEMRNU = TEMNR*TEMDR + TEMNI*TEMDI 



317 TEMINU = TEMNI*TEMDR - TEMNR+TEMDI 



318 DELTA = TEMRNU/TEMDEN 



319 DELTI = TEMINU/TEMDEN 



320 c * * * THE NEXT CONSTANT DEPENDS ON WORD LENGTH AND SIZE OF PHASE VELOCITY 



321 IF (ABS(DELTA) . LT . 1.D-14) GO TO 51 



322 SIZE = DELTA*DELTA + DELTI*DELTI 



323 IF ( (SIZE.GT.SIZE2) .AND.( J.GT.3) ) GO TO 51 



324 92 SIZE2 = SIZE * 2. 

 335 VEL - DET 



326 VELI = DETI 



327 GO TO 48 



328 C FIND DEPTH FUNCTIONS 



329 51 IF (INSUR .EQ. 0) GO TO 61 



330 TRE = (DET**2 + DETI**2) / RX 



331 IF (TRE .LT. 1E-10) GO TO 61 



332 PRINT 998, NN , TRE 



333 998 FORMAT ( 5H MODE ,I4,23H FAILED TO CONVERGE — , E9.2) 



334 GO TO 999 



335 61 IF (MPCH . EQ . 0) GO TO 63 



336 IF (NXTRA . LT . 0) GO TO 63 



337 TEM1 = V * 1 .04 



338 C0L(1 ) = TEM1 



339 TEMP = C0L(1 ) 



340 C0L(2) = (TEM1 - TEMP) * 1.D10 



341 TEM1 = VI * 1 .04 



49 



