342 C0L(3) = TEM1 



343 TEMP = C0H3) 



344 C0L(4) = (TEM1 - TEMP) * 1.D10 



345 C0L(5) = -NN 



346 PUNCH 64, (COL(I), I = 1,5) 



347 64 FORMAT (5110) 



348 63 AF( 1 ,NN) = A( 1 ,3) 



349 AG(1 ,NN) = 0(1,3) 



350 BF( 1 ,NN) =-A( 1 ,4) 



351 BG( 1 ,NN) =-Q( 1 ,4) 



352 PHASE V(NN) = V 



353 PHASI V(NN) = VI 



354 IF (KB .EO. 1 ) GO TO 73 



355 74 PRINT 81, V, VI, DET, DET I , SIZE, CNTR 



356 73 LL = N - 1 



357 IF (LL-1 ) 95,96,97 



358 96 1=0 



359 GO TO 98 



360 97 DO 110 J = 2, LL 



361 1 = J + J - 2 



362 TEMM? :- A(I,2)+AF( J-1 ,NN) - 0(I,2)*AG( J-1 ,NN) 



363 TEMNI = 0(I,2)+AF( J-1 ,NN) + A(I,2)*AG( J-1 ,NN) 



364 TEMDR = A(I,3)*A( 1+1 ,4) - 0(I,3)*0( 1+1 .4) - 



365 1 A(I,4)*A( 1+1 ,3) +0(I,4)*0( 1+1 ,3) 



366 TEMDI = 0(I,3)*A( 1+1 ,4) + A(I,3)*0( 1+1 ,4) - 



367 1 0(I,'?)*A( 1 + 1 ,3) -A(I,4)*(3( 1 + 1 ,3) 



368 TEMDEN = TEMDR+TEMDR + TEMDI*TEMDI 



369 TEMRNU = TEMNRtTEMDR + TEMNI*TEMDI 



370 TEMINU = TEMNI*TEIV1DR - TEMNR*TEMDI 



371 TEMP = TEMRNU / TEMDEN 



372 TEMPI = TEMINU / TEMDEN 



373 BF(J,NN) = -(TEMP*A( 1+1 ,4) - TEMPI*0( 1+1 ,4)) 



374 BG(J,NN) = -(TEMPI*A( 1+1 ,4) + TEMP*0( 1+1 ,4)) 



375 AG(J,NN) = TEMPI*A( 1+1 ,3) + TEMP+0( 1+1 .3) 



376 110 AF(J,NN) = TEMP*A( 1+1 ,3) - TEMPI*0( 1+1 ,3) 



377 98 TEMNR = -(A(I+2,2) * AF(LL,NN) -0(1+2,2) * AG(LL,NN)) 



378 TEMNI = -(0( 1+2 , 2 ) * AF ( LL . NN ) + A( 1+2 , 2 ) ♦ AG( LL , NN ) ) 



379 TEMDEN = A( 1+2 ,3)*A( 1+2 ,3) + 0( 1+2 ,3)*0( 1+2 ,3) 



380 TEMRNU = TEMNR*A( 1+2 ,3) + TEMNI*0( 1+2 ,3) 



381 TEMINU ^ TEMNI*A( 1+2 ,3) - TEMNR*0( 1+2 ,3) 



382 BF(N,NN) = TEMRNU / TEMDEN 



383 BG(N,NN) = TEMINU / TEMDEN 



384 95 AF(N,NN) = 0. 



385 AG(N,NN) = 0. 



386 C FIND NORMALIZING FACTOR 



387 D(NN) = 2.12429296D0 * RH0(1)**3 / 0(1) 



388 DI(NN) = 0. 



389 DO 111 I = 2,N 



390 TEMRSP = AF( 1-1 ,NN)*B( 2*1-2 ,2) - AG( 1-1 ,NN)*BI( 2*1-2 ,2) 



391 1 BF( 1-1 ,NN)*B( 2*1-2 ,1) - BG( 1-1 ,NN)*BI( 2*1-2 ,1) 



392 TEMISP = AG( 1-1 ,NN)*B( 2*1-2 ,2) + AF( 1-1 ,NN)*BI( 2*1-2 ,2) 



393 1 BG( 1-1 ,NN)*B( 2*1-2 ,1) + BF( 1-1 ,NN)*BI( 2*1-2 ,1) 



394 AX1 = TEMRSP*TEMRSP - TEMISP*TEMISP. 



395 AX1I = TEMRSP * TEMISP 



396 AX1I = AX1 I + AX1 I 



397 TEMDR = (G(I-1)**2 + GI(I-1)**2) 



398 TEMDI = G(I)**2 + GI(I)**2 



50 



