1 





SObKOUllNt UtTNT(N,DET,DETI) 



2 





IMPLICIT DOUBLE PRECISION (A-H,0-Z) 



3 





DOUBLE PRECISION DET, DET I 



4 





COMMON /EXPO/ EXSUM, CNTR , RATI0(25) 



5 





COMMON/DETMNT/ A(25,4), 0(25,4) 



6 





DLOSS = 1 . 



7 





CNTR = 0. 



8 





DET = A(1 ,3) 



9 





DETI = 0(1 ,3) 



10 





LIM = N + N - 3 



11 





DO 100 1=1 ,LIM,2 



12 





J = I 



13 





K = 4 



14 





L = I 



15 





M = 3 



16 





II = 1 



17 





GO TO 500 



18 



10 



d = I + 1 



19 





K = 2 



20 





L = d 



21 





M = 1 



22 





GO TO 600 



23 



30 



d = I + 2 



24 





L = d 



25 





II = 2 



26 





GO TO 600 



27 



40 



L = I + 1 



28 





M = 2 



29 





GO TO 500 



30 



50 



K = 3 



31 





M = 3 



32 





II = 3 



33 





GO TO 600 



34 



60 



K = 4 



35 





IF (I .EO. LIM) GO TO 70 



36 





M = 4 



37 





11=4 



38 





GO TO 600 



39 



70 



K = 3 



40 





II = 2 



41 





GO TO 700 



42 



500 



C = A(L,M)*A( L.M) + 0(L,M)*0(L,M) 



43 



80 



B = (A(d,K)*A(L,M) + 0(d.K)*0(L,M) ) / C 



44 





BI = (0(d,K)*A(L,M) - A(d,K)*0(L,M))/C 



45 





GO TO ( 10,50) , II 



46 



600 



TD = A(d,K) - (A(L,M)*B - 0(L,M)*BI) 



47 





TDI = Q(d,K) - (A(L,M)*BI + 0(L,M)*B) 



48 





TEM = TD**2 + TDI**2 



49 





TEMP = A(d,K)**2 + 0(d,K)*»2 



50 





TEMP = TEM / TEMP 



51 





IF (II .EO. 2) GO TO 92 



52 





IF (II .EO- 4) GO TO 92 



53 





0(d,K) = 0(d,K) * 10.D-18 



54 





A(d,K) = A(d,K) * 10.D-18 



55 





IF (TEMP .GT. 10.D-35) GO TO 92 



56 





CNTR = CNTR + 1 . 



61 



