PROGRAM TOKTIN
C   Ž–…Š€ Š€—…‘’‚€ €„ˆŽ‹…Š’ŽŽ‰ €€€’“›
C
C   NM - ŠŽ‹-‚Ž ‘’Ž‹–Ž‚, Œ - ŠŽ‹-‚Ž ‘’ŽŠ
C
      CHARACTER *15 FIL,FREZ
      DIMENSION •1 •(50,20),•‚(50,20),XRD(20,10),
     1 DEL(20,20),DR(10,10),DPN(20,20),DZ(20,20) ,ZAP(20,20),
     2 CX(20,20),AB(20,20),TG(20,20),TIM(20,20),BK(20,20),
     3 D(20),S(20),KL(20),XCP(20),DIS(20),AKL(20),KOL(20)
      DATA NK/3/
      WRITE(*,*)'‚ˆŒ€ˆ…!  ‚ Ž–…‘‘… €Ž’› ‘Ž‡„€…’‘Ÿ ”€‰‹'
      WRITE(*,*)'C …‡“‹œ’€’€Œˆ, € Š€ ˆ”ŽŒ€–ˆŸ … ‚›„€…’‘Ÿ!'
      WRITE(*,*)' '
      WRITE(*,*)'‚‚…„ˆ’… ˆŒŸ ”€‰‹€ …‡“‹œ’€’Ž‚: '
      READ(*,9) FREZ
      OPEN(UNIT=3,FILE=FREZ,STATUS='NEW')
      WRITE(3,*)'------------------------------------------------'
      WRITE(3,*)'   1-Š ”€‰‹ - Œ€’ˆ–€ ‡€—…ˆ‰ €€Œ…’Ž‚       '
      WRITE(3,*)'           ‚ €’“€‹œŽŒ ‚ˆ„…                   '
      WRITE(3,*)'------------------------------------------------'
      WRITE(3,*)'                                                (
      WRITE(6,2147)
 2147 FORMAT(/10X,'‚‚…„ˆ’… €‡Œ…Ž‘’œ 1-ƒŽ Œ€‘‘ˆ‚€ '/)
      WRITE(6,2148)
      WRITE(6,*)'  '
 2148 FORMAT(/10X,'Œ - ‘’ŽŠ,      NM - ‘’Ž‹–Ž‚    '/)
      READ(*,*) M,NM
      WRITE(3,*) M,NM
      WRITE(6,91)
 91   FORMAT(' ‚‚…„ˆ’… ˆŒŸ 1-ƒŽ ”€‰‹€ „€›•: ')
      READ(*,9) FIL
 9    FORMAT(A15)
      WRITE(3,9) FIL
      OPEN(UNIT=1,FILE=FIL)
      WRITE(3,3000)
 3000 FORMAT(/15X,'ˆ‘•Ž„›… „€›…'/16•,'1-€Ÿ ’€‹ˆ–€')
      DO 3001 I=1,Œ
      READ(1,*) (XIX(I,J),J=1,NM)
 3001 WRITE(3,3002) I,(XIX(I,J),J=1,NM)
 3002 FORMAT(I4,'.  ',10F9.2)
      CLOSE(1)
      WRITE(3,*)'---------------------------------------------'
      WRITE(3,*)'    2-ŽŠ ”€‰‹ - Œ€’ˆ–€ ŽŒˆ€‹œ›•,
      WRITE(3,*)' „Ž“‘’ˆŒ›• ˆ €‡Ž‚›• ‡€—…ˆ‰ €€Œ…’Ž‚
      WRITE(3.*)'---------------------------------------------'
      WRITE(3,*)'  '
      WRITE(*,333)
 333  FORMAT(' ‚‚…„ˆ’… ˆ„Ÿ 2-ƒŽ ”€‰‹€ „€›•: ')
      READ(*,222) FIL
 222  FORMAT(A15)
      WRITE(3,222) FIL
      OPEN(UNIT=1,FILE=FIL)
      WRITE(3,6112)
 6112 FORMAT(/15X,'2-€Ÿ ’€‹ˆ–€')
      DO 6111 I=1,NM
      READ(1,*) (XRD(I,J),J=1,NK)
 6111 WRITE<3,39) I,(XRD(I,J),J=1,NK)
 39   FORMAT(I4,'.  ',10F9.2)
      CLOSE(1)
      WRITE(*,*) '‚‚…„ˆ’… ŽŒ…€ ƒ“ €€Œ…’Ž‚:
      READ(*,*) (KOL(K),K=1,NM)
      WRITE(3,243)
 243  FORMAT(/'            ŽŒ…€ ƒ“ €€Œ…’Ž‚             ')
      WRITE<3,*)'          ',(KOL(K),K=1,NM)
C
C     …‡€‡Œ…›… ‡€—…ˆŸ €€Œ…’Ž‚
C
      DO 905 I=1,NM
      T=XRD(I,3)
      IF (KOL(I).EQ.2) GOTO 906
      DO 907 J=1,M
 907  XIX(J,I)=XIX(J,I)/T
      GOTO 905
 906  DO 908 J=1,M
 908  XIX(J,1)=T/XIX(J,1)
 905  CONTINUE
      WRITE(3,600)
 600  FORMAT(//10X,'BE3PA3METHNE ‡€—…ˆ… €€Œ…’Ž‚')
      DO 601 J=1,M
 601  WRITE(3,602) J,(XIX(J,I),I=1,NM)
 602  FORMAT(1X,I3,'.',2X,10F11.5)
C
C     ‘…„ˆ… ‡€—…ˆŸ €€Œ…’Ž‚
C
      AN=M
      DO 107 J=1,NM
      XCP(J)=0.
      DO 106 I=1,M
 106  XCP(J)=XCP(J)+XIX(I,J)
 107  XCP(J)=XCP(J)/AN
C
C     ‡€—…ˆŸ „ˆ‘…‘ˆ‰ €€Œ…’Ž‚
C
      AN=AN-1
      DO 108 J=1,NM
      D(J)=0.
      DO 109 I=1,M
      XIX(I,J)-XCP(J)
 109  D(J)=D(J)+T*T
      D(J)=D(J)/AN
 108  S(J)=SQRT(D(J))
C
C     ‘’€„€’ˆ‡Ž‚€›… ‡€—…ˆŸ €€Œ…’Ž‚
C
      DO 110 J=1,NM
      DO 110 I=1,Œ
      XB(I,J)=(XIX(I,J)-XCP(J))/S(J)
      WRITE(3,610)
 610  FORMAT(/8X,'*** ‘…„ˆ… ‡€—…ˆŸ €€Œ…’Ž‚ ***')
      WRITE(3,611) (XCP(J),J=1,NM)
 611  FORMAT(10(1X,F11.5))
      WRITE(3,612)
 612  FORMAT(/10X,'‡€—…ˆŸ „ˆ‘…‘ˆ‰ €€Œ…’Ž‚')
      WRITE(3,611) (D(J),J=1,NM)
      WRITE(3,613)
 613  FORMAT(/10X,'‘…„…Š‚€„€’ˆ—…‘Šˆ… Ž’Š‹Ž…ˆŸ €€Œ…’Ž‚')
      WRITE(3,611) (S(J),J=1,NM)
      WRITE(3,614)
 614  FORMAT(/10X,'‘’€„€’ˆ‡Ž‚€›… ‡€—…ˆŸ €€Œ…’Ž‚')
      DO 615 I=1,M
 615  WRITE(3,616) I,(Xe(I,J),J=1,NM)
 616  FORMAT(1X,I3,'.',2X,10F11.5)
C
C     ŠŽ””ˆ–ˆ…’› ŠŽ…‹Ÿ–ˆˆ Œ…†„“ €€Œ…’€Œˆ
C
      AN=AN+1.
      J11=NM-1
      DO 112 L=1,J11
      J=L+1
      DO 112 K=J,NM
      T=0.
      DO 111 I=1,M
 111  T=T+XB(I,L)*XB(I,K)
      T=T/AN
      DEL(L,K)=T
      DEL(K,L)=T
 112  CONTINUE
      R=(2.*NM-SQRT(8.*NM+1.)+1.)/2.
      IZ=INT(R+0.5)
      WRITE(3,625) IZ
 625  FORMAT(/10X,'ŠŽ‹ˆ—…‘’‚Ž ŽŽ™…›• ŽŠ€‡€’…‹…‰  IZ= ',I3)
C
C     Ž“‹…ˆ… „ˆ€ƒŽ€‹œ›• ‹…Œ…’Ž‚ ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›
C
      DO 133 LZ=1,IZ
      WRITE(3,854) LZ
 854  FORMAT(///'      Ž˜‹ˆ €‘—…’› „‹Ÿ ',I3,' ŽŠ€‡€’…‹Ÿ'/)
      DO 139 I=1,NM
 139  DEL(I,I)=0.
      DO 151 I=1,NM
      DO 151 J=1,NM
 151  DR(I,J)=DEL(I,J)
      IF (LZ.EQ.1) WRITE(3,*)'          ŠŽ…‹Ÿ–ˆŽ€Ÿ Œ€’ˆ–€'
      IF (LZ.NE.1) WRITE(3,741) LZ
 741  FORMAT('          Ž‘’€’Ž—€Ÿ ŠŽ…‹Ÿ–ˆŽ€Ÿ Œ€’ˆ–€'/
     * '           „‹Ÿ ',I3,' ŽŠ€‡€’…‹Ÿ')
 176  DO 628 I=1,NM
 628  WRITE(3,629) I,(DR(I,J),J=1,NM)
 629  FORMAT(1X,I3,'.',2X,10F11.5)
C
C     ‘“ŒŒ› ‹…Œ…’Ž‚ Ž ‘’Ž‹–€Œ ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›
C
 120  DO 116 J=1,NM
      XCP(J)=0.
      KL(J)=1
      DO 115 I=1,NM
 115  XCP(J)=XCP(J)+DEL(I,J)
      IF (XCP(J).LT.0.) KL(J)=-1
 116  WRITE(3,99995) KL(J)
99995 FORMAT(1X,'*** Š ***',120)
      DO 78877 J=1,NM
78877 IF (KL(J).EQ.-1) GOTO 758
      GOTO 121
C     …‘‹ˆ Š=-1 (…‘’œ Œˆ“‘›) ˆ„…Œ € Ž€™…ˆ…
 758  WRITE(3,*)' '
      WRITE(3,*)'     ‘“ŒŒ› Ž ‘’Ž‹–€Œ ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›'
      WRITE(3,811) (XCP(J),J=1,NM)
 811  FORMAT(5X,10(1X,F11.5))
C
C    ŒˆˆŒ€‹œ€Ÿ Ž’ˆ–€’…‹œ€Ÿ ‘“ŒŒ€
C       ˆ‡ ‘’Ž‹–Ž‚ ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›
C
      T = 0.
      DO 117 J=1,NM
      IF (KL(J).EQ.1) GOTO 117
      IF (T-XCP(J)) 117,117,118
 118  T=XCP(J)
      L=J
 117  CONTINUE
      WRITE(3,4) XCP(L)
 4    FORMAT(/10X,'ŒˆˆŒ€‹œ€Ÿ Ž’ˆ–€’…‹œ€Ÿ ‘“ŒŒ€ ˆ‡ /10X,
     * '‘’Ž‹–Ž‚ ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›  ',F11.5)
      DPN(LZ,L)=-1
C
C     Ž€™…ˆ… ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›
C
      DO 119 I=1,NM
 119  DEL(I,L)=-DEL(I,L)
      DO 200 I=1,NM
 200  DEL(L,I)=-DEL(L,I)
      WRITE(3,188)
 188  FORMAT(/20X,'Ž€™…ˆ… ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›')
      DO 51 1=1,NM
 51   WRITE(3,1551) (DEL(I,J),J=1,NM)
 1551 FORMAT(11X,10F11.5)
      GOTO 120
C
C     Œ€Š‘ˆŒ€‹œ›… ‹…Œ…’› ‚ ƒ‹€‚“ž
C         „ˆ€ƒŽ€‹œ ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›
C
 121  DO 127 J=1,NM
      T=ABS(DEL(1,J))
      DO 126 I=1,NM
      IF (T-ABS(DEL(I,J))) 125,126,126
 125  T=ABS(DEL(I,J))
 126  CONTINUE
 127  DEL(J,J)=T
      WRITE(3,7)
 7    FORMAT(/6X,'Œ€Š‘ˆŒ€‹œ›… ‹…Œ…’› ‚ ‘’Ž‹–…'/
     * 10•,'ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–› - ‹…Œ. Ž „ˆ”ƒŽ€‹ˆ')
      WRITE(3,3434) (DEL(J,J),J=1,NM)
 3434 FORMAT(10X,10F11.5)
      WRITE(3,761)
 761  FORMAT(/5X,'Œ€’ˆ–€ ‘ Ž„‘’€‚‹…›Œˆ Œ€Š‘ˆŒ€‹œ›Œˆ ',
     * '‹…Œ…’€Œˆ')
      DO 605 I=1,NM
 605  WRITE(3,604) (DEL(I,J),J=1,NM)
 604  FORMAT(10X,10F11.5)
C
C     ‚›„…‹…ˆ… ŽŽ™…›• ŽŠ€‡€’…‹…‰
C
C     ‘“ŒŒ€ ‚ ‘’Ž‹–€• ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›
C                 ‘ “—…’ŽŒ „ˆ€ƒŽ€‹œ›• ‹…Œ…’Ž‚
C
      DO 754 J=1,NM
      XCP(J)=0.
      DO 754 I=1,NM
 754  XCP(J)=XCP(J)+DEL(I,J)
      WRITE(3,756)
 756  FORMAT(/' ‘“ŒŒ€ ‚ ‘’Ž‹–€• ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–›
     * ‘ “—…’ŽŒ „ˆ€ƒŽ€‹œ›• ‹…Œ…’Ž‚')
      WRITE(3,755) (XCP(J),J=1,NM)
 755  FORMAT(10X,10F11.5)
C
C     •ˆ’›‰ ŠŽ…œ, ‘€Œ ŠŽ…œ ‘ŒŽ’ˆ ˆ†…
      T=0.
      DO 1092 I=1,NM
 1092 T=T+XCP(I)
      WRITE(3,*)' '
      WRITE(3,*)'   ‚›„…‹…›… ŽŽ™…›… ŽŠ€‡€’…‹ˆ'
      DO 128 J=1,NM
      DZ(LZ,J)=XCP(J)/SQRT(T)
 128  WRITE(3,*)'   ',DZ(LZ,J)
C
C     Ž‘’€’Ž—€Ÿ ŠŽ…‹Ÿ–ˆŽ€Ÿ Œ€’ˆ–€
C
      DO 523 I=1,NM
      DO 523 J=1,NM
 523  DEL(I,J)=ABS(DEL(I,J))-ABS(DZ(LZ,I)*DZ(LZ,J))
      DO 623 I=1,NM
      DO 623 J=1,NM
 623  IF (DR(I,J).LT.0) DEL(I,J)=-DEL(I,J)
      WRITE(3,3003)
 3003 FORMAT(/10X,'ŠŽ…‹Ÿ–ˆŽ€Ÿ Œ€’ˆ–€ ŽŽ™…›• ŽŠ€‡€’…‹…‰')
 3004 WRITE(3,3005) LZ,(DZ(LZ,J),J=1,NM)
 3005 FORMAT('  LZ=',I3,'.  ',10F11.5)
C
C     ‚Ž‘‘’€Ž‚‹…ˆ… ‡€ŠŽ‚ ‚ ŠŽ…‹Ÿ–ˆŽŽ‰ Œ€’ˆ–…
C        ŽŽ™…›• ŽŠ€‡€’…‹…‰
C
      DO 142 J=1,NM
      IF (DPN(LZ,J)) 141,142,142
 141  K=J
      DO 143 L=K,IZ
 143  DZ(LZ,L)=-DZ(LZ,L)
 142  CONTINUE
      WRITE(3,640)
 640  FORMAT(/5X,' ŠŽ…‹Ÿ–ˆŽ€Ÿ Œ€’ˆ–€ ŽŽ™…›• ŽŠ€‡€’…‹…‰'/
     * 10•,' ‘ ‚Ž‘‘’€Ž‚‹…›Œˆ ‡€Š€Œˆ')
 641  WRITE(3,642) LZ,(DZ(LZ,I),I=1,NM)
 642  FORMAT('  LZ=',I3,'.  ',10F11.5)
C
C     ŠŽ…‹Ÿ–ˆŽ›… ‚Š‹€„›
C
      DO 144 J=1,NM
      XCP(J)=0.0
 144  XCP(J)=XCP(J)+DZ(LZ,J)*DZ(LZ,J)
      WRITE(3,4000)
 4000 FORMAT(/20X,' ŠŽ…‹Ÿ–ˆŽ›… ‚Š‹€„›')
      WRITE(3,4001) (XCP(J),J=1,NM)
 4001 FORMAT(10(1X,F11.5))
C
C     …˜€…Œ Œ€’ˆ–“, ˆ™…Œ ŠŽ””ˆ–ˆ…’› ŠŽ…‹Ÿ–ˆˆ
C
C     Ž„ƒŽ’Ž‚Š€ ‘ˆ‘’…Œ› “€‚…ˆ‰
C
      DO 145 I=1,NM
 145  DR(I,I)=1.0
      K=NM-1
      NM1=NM+1
      DO 321 I=1,NM
 321  D(I)=DZ(LZ,I)
C   DR,D,NM - TO —’Ž €„Ž, A AB,AKL - €‹‹€‘’ „‹Ÿ ‹“—˜…‰ …˜€…ŒŽ‘’ˆ
      CALL MATR(DR,D,NM,AB,AKL,DET)
      IF (DET.EQ.0.) WRITE(*,*)' Ž…„…‹ˆ’…‹œ ‘ˆ‘’…Œ› €‚… “‹ž'
      IF (DET.EQ.0.) GOTO 910
      DO 1534 I=1,NM
 1534 BK(LZ,I)=D(I)
      WRITE(3,345)
 345  FORMAT(/10X, 'ŠŽ””ˆ–ˆ…’› …ƒ…‘‘ˆˆ €€Œ…’Ž‚
     * ‚ ‘’€„€’ˆ‡Ž‚€ŽŒ ‚ˆ„…  ')
      WRITE(3,647) LZ,(BK(LZ,J),J=1,NM)
 647  FORMAT('  LZ=',I3,'.',1X,10F12.4)
C
C     ŠŽ””ˆ–ˆ…’› …ƒ…‘‘ˆˆ €€Œ…’Ž‚ ‚ …‡€‡Œ…ŽŒ ‚ˆ„…
C
      DO 155 I=1,NM
 155  BK(LZ,I)=BK(LZ,I)/S(I)
      WRITE(3,4030)
 4030 FORMAT(/10X,'ŠŽ””ˆ–ˆ…’› …ƒ…‘‘ˆˆ €€Œ…’Ž‚
     * ‚ …‡€‡Œ…ŽŒ ‚ˆ„…')
      WRITE(3,647) LZ,(BK(LZ,J),J=1,NM)
C
C     ŠŽ””ˆ–ˆ…’› …ƒ…‘‘ˆˆ ‚ €’“€‹œŽŒ ‚ˆ„…
C     (2,4) TIM - ŠŽ””ˆ–ˆ…’ ‚ €’“€‹œŽŒ ‚ˆ„…... €‡Ž‚›‰...
C
      DO 160 I=1,NM
      IF (KOL(I).EQ.1) TIM(LZ,I)=BK(LZ,I)/XRD(I,3)
 160  IF (KOL(I).EQ.2) TIM(LZ,I)=BK(LZ,I)*XRD(I,3)
      WRITE(3,7413)
 7413 FORMAT(/20X,'ŠŽ””ˆ–ˆ…’› …ƒ…‘‘ˆˆ ‚ €’“€‹œŽŒ ‚ˆ„… ')
      WRITE(3,647) LZ,(TIM(LZ,I),I=1,NM)
C
C     Œ€’ˆ–€ ŽŽ™…›• ŽŠ€‡€’…‹…‰
C
      DO 1624 I=1,NM
      DO 1624 J=1,NM
      TG(I,J)=0.
 1624 CONTINUE
      DO 1625 K=1,3
      DO 1625 I=1,NM
      IF(KOL(I).EQ.1) TG(LZ,K)=TG(LZ,K)+TIM(LZ,I)*XRD(I,K)
 1625 IF(KOL(I).EQ.2) TG(LZ,K)=TG(LZ,K)+TIM(LZ,I)/XRD(I,K)
      WRITE(3,1010) LZ
 1010 FORMAT(//10X,'    ŽŽ™…›… ŽŠ€‡€’…‹ˆ ‚ €’“€‹œŽŒ ‚ˆ„…'/
     1 '           ŽŒˆ€‹œ›… ... „Ž“‘’ˆŒ›… ... €‡Ž‚›… ... '/
     2 /'           Œ€’…Œ€’ˆ—…‘Š€Ÿ ŒŽ„…‹œ (LZ) ŽŒ… ',I2)
      WRITE(3,1622) (TG(LZ,Š),K=1,3)
 1622 FORMAT(10X,3F11.5)
      L=0
      DO 1 I=1,3
      IF(TG(LZ,I).LT.0.) L=L+1
 1    CONTINUE
      IF(L.EQ.O) GOTO 3
      P=TG(LZ,1)
      IF(P.GT.TG(LZ,2)) P=TG(LZ,2)
      IF(P.GT.TG(LZ,3)) P=TG(LZ,3)
      P=-2.*P
      DO 2 I=1,NK
 2    TG(LZ,I)=TG(LZ,I)*P
 652  FORMAT(1X,'LZ=',I3,'.  ',3F11.5)
C
C     Ž’Ž‘ˆ’…‹œ›… ‡€—…ˆŸ ŽŽ™…›• ŽŠ€‡€’…‹…‰
C
 3    WRITE(3,650)
 650  FORMAT(/8X,'Œ€’ˆ–€ ŽŽ™…›• ŽŠ€‡€’…‹…‰    (’Ž†… ‘€ŒŽ…)'/
     * 13•,'‘Ž ‘„‚ˆƒŽŒ €—€‹€ ŠŽŽ„ˆ€’')
      WRITE(3,652) LZ,(TG(LZ,K),K=1,3)
      IF(TG(LZ,1)-TG(1,2)) 162,164,163
 162  CX(LZ,1)=TG(LZ,3)/TG(LZ,1)
      CX(LZ,2)=TG(LZ,3)/TG(LZ,2)
      GOTO 165
 163  CX(LZ,1)=TG(LZ,1)/TG(LZ,3)
      CX(LZ,2)=TG(LZ,2)/TG(LZ,3)
 165  WRITE(3,655)
 655  FORMAT(/10X,'Ž’Ž‘ˆ’…‹œ›… ‡€—…ˆŸ ŽŽ™…›• ŽŠ€‡€’…‹…‰')
      WRITE(3,656) LZ,(CX(LZ,I),I=1,2)
 656  FORMAT('  LZ=',I3,'.    ',2F11.5)
C
C     Ž–…Šˆ ŽŽ™…›• ŽŠ€‡€’…‹…‰
C
      DO 170 I=1,2
      R=ALOG10(CX(LZ,I)+1.)/ALOG10(2.)
      R=1./EXP(1./EXP(R))
 170  CX(LZ,I)=R
      WRITE(3,660)
 660  FORMAT(/20X,'Ž–…Šˆ ŽŽ™…›• ŽŠ€‡€’…‹…‰')
      WRITE(3,656) LZ,(CX(LZ,I),I=1,2)
 133  CONTINUE
C
C     ‚…‘Ž‚›… ŠŽ””ˆ–ˆ…’› ŽŽ™…›• ŽŠ€‡€’…‹…‰
C
      ’=0.
      DO 171 LZ=1,IZ
 171  T=T+1./ABS.(CX(LZ,1)-CX(LZ,2))
      DO 172 LZ=1,IZ
 172  XCP(LZ)=1./ABS(CX(LZ,1)-CX(LZ,2))/T
      ’=0.
      DO 1721 LZ=1,IZ
 1721 T=T+ABS(1./ALOG10(CX(LZ,1)/CX(LZ,2)))
      DO 1722 LZ=1,IZ
 1722 AKL(LZ)=1./ABS(ALOG10(CX(LZ,1)/CX(LZ,2)))/T
      WRITE(3,665)
 665  FORMAT(///10X,'‚…‘Ž‚›… ŠŽ””ˆ–ˆ…’› ŽŽ™…›• ŽŠ€‡€’…‹…‰'
     * /15•,'„‹Ÿ Œ“‹œ’ˆ‹ˆŠ€’ˆ‚Ž‰ ŒŽ„…‹ˆ')
      DO 6652 LZ=1,IZ
 6652 WRITE(3,6651) LZ,AKL(LZ)
 6651 FORMAT(15X,'X= ',I3,'.  ',F11.5)
      WRITE(3,6653)
 6653 FORMAT(/15•,'„‹Ÿ €„„ˆ’ˆ‚Ž‰ ŒŽ„…‹ˆ')
      DO 6654 LZ=1,IZ
 6654 WRITE(3,66501) LZ,XCP(LZ)
66501 FORMAT(15X,'A= ',I3,'. .',F11.5)
C
C     €‘—…’ „ˆ•Ž’ŽŒˆ—…‘ŠŽ‰ ”“Š–ˆˆ
      CALL DIM(M,KOL,NM,XDR)
C
C     ŠŽŒ‹…Š‘›‰ Šˆ’…ˆ‰ Š€—…‘’‚€
C     D(1) - „‹Ÿ Œ“‹’ˆ‹ˆŠ€’ˆ‚Ž‰ ŒŽ„…‹ˆ  D(2) - „‹Ÿ €„„ˆ’ˆ‚Ž‰ ŒŽ„…‹ˆ
C
      D(1)=1.
      D(2)=0.
      DO 173 LZ=1,IZ
      D(1)=D(1)*(AKL(LZ)**XCP(LZ))
 173  D(2)=D(2)+(AKL(LZ)*XCP(LZ))
      WRITE(3,670) D(1),D(2)
 670  FORMAT(/20X,'ŠŽŒ‹…Š‘›‰ Šˆ’…ˆ‰ Š€—…‘’‚€ (0,6922) '/
     1 10•,'„‹Ÿ Œ“‹’ˆ‹ˆŠ€’ˆ‚ŽŠ ŒŽ„…‹ˆ =',F11.5/
     2 10•,'„‹Ÿ   €„„ˆ’ˆ‚Ž‰     ŒŽ„…‹ˆ =',F11.5)
C
C     „ˆ‘…‘ˆŸ ŠŽŒ‹…Š‘ŽƒŽ Šˆ’…ˆŸ Š€—…‘’‚€
C
      AN=NM-1
      DO 909 I=1,2
      DIS(I)=0.
      DO 909 LZ=1,IZ
 909  DIS(I)=DIS(I)+(CX(LZ,I)-D(I))*(CX(LZ,I)-D(I))
      DO 666 I=1,2
      T=DIS(I)/AN
 666  S(I)=SQRT(T)
      WRITE(3,675) DIS(1),DIS(2)
 675  FORMAT(/5X,'„ˆ‘…‘ˆŸ ŠŽŒ‹…Š‘ŽƒŽ Šˆ’…ˆŸ Š€—…‘’‚€'
     1 /10•,'„‹Ÿ Œ“‹’ˆ‹ˆŠ€’ˆ‚ŽŠ ŒŽ„…‹ˆ = '.F11.5
     2 /10•,'„‹Ÿ    €„„ˆ’ˆ‚Ž‰    ŒŽ„…‹ˆ = ',F11.5)
      WRITE(3,677) S(1),S(2)
 677  FORMAT(/2•,'‘…„…Š‚€„€’ˆ—…‘ŠŽ… Ž’Š‹Ž…ˆ… ŠŽŒ‹…Š‘ŽƒŽ
     1 Šˆ’…ˆŸ Š€—…‘’‚€ '/20•,'€€‹Žƒˆ—Ž  ',F11.5/32X,F11.5)
      WRITE(5,5020)
 5020 FORMAT(' ‚‚…„ˆ’… Šˆ’…ˆ‰ ‘’œž„…’€ ’‘’= ')
      READ(*,5000) ’‘’
 5000 FORMAT(F11.5)
      WRITE(3,5111) ’‘’
 5111 FORMAT(/2X,'Šˆ’…ˆ‰ ‘’œž„…’€   ’‘’= ',F11.5)
      AN=NM
      DO 5007 I=1,2
      DIS(I)=TCT*S(I)/SQRT(AN)
 5007 D(I)=DIS(I)/D(I)
      WRITE(3,680)
 680  FORMAT(/10X,'€‘Ž‹ž’›… ˆ Ž’Ž‘ˆ’…‹œ›… Žƒ…˜Ž‘’ˆ
     1 ŠŽŒ‹…Š‘ŽƒŽ Šˆ’…ˆŸ Š€—…‘‚€')
      DO 7412 I=1,2
 7412 WRITE(3,681) DIS(I),D(I)
 681  FORMAT(10X, '€‘Ž‹ž’. ',F11.5,'    Ž’Ž‘ˆ’. ',F11.5/)
      GOTO 910
 164  WRITE(3,911)
 911  FORMAT(/'ŽŒˆ€‹œ›… ˆ „Ž“‘’ˆŒ›… ‡€—…ˆŸ ŽŽ™…›•
     * ŽŠ€‡€’…‹…‰ €‚›')
 910  WRITE(*,*) '…‡“‹œ’€’› ‘ŒŽ’ˆ ‚ ”€‰‹…   ',FREZ
      CLOSE(3)
      STOP 'ŠŽ…– €Ž’›'
      END
C
      SUBROUTINE MATR(AO,BO,N,AA,X,DET)
      DIMENSION AO(N,N),BO(N),AA(N,N),X(N)
      DO 147 I=1,N
      DO 147 J=1,N
 147  AA(I,J)=AO(I,J)
      CALL OPER(AA,N,O)
      DET=O
      IF(DET.EQ.0.) RETURN
      DO 7 K=1,N
      DO 1 I=1,N
      DO 1 J=1,N
 1    AA(I,J)=AO(I,J)
      DO 2 I=1,N
 2    AA(I,K)=BO(I)
      CALL OPER(AA,N,T)
      DET=T
      IF(DET.EQ.0.) RETURN
      X(K)=T
 7    CONTINUE
      DO 6 K=1,N
 6    X(K)=X(K)/O
      DO 142 I=1,N
 142  BO(I)=X(I)
      RETURN
      END
C
      SUBROUTINE OPER(A,N,DET)
      DIMENSION A(N,N)
      K=1
      DET=1.
 8    L=K
      AMA=ABS(A(L,K))
      LROW=K
 1    L=L+1
      AMB=ABS(A(L,K))
      IF(AMA.GE.AMB) GOTO 2
      LROW=L
      AMA=AMB
 2    IF(L.LT.N) GOTO 1
      IF(LROW.NE.K) GOTO 3
      SIGN=1.
      GOTO 6
 3    SIGN=-1.
      DO 5 J=K,N
      S=A(K,J)
      A(K,J)=A(LROW,J)
 5    A(LROW,J)=S
      IF(A(K,K).NE.0) GOTO 6
      WRITE(*,1315)
 1315 FORMAT(1H0,10X,'Ž…„…‹ˆ’…‹œ ‘ˆ‘’…Œ› €‚… “‹ž')
      DET=0.
      RETURN
 6    DET=SIGN*A(K,K)*DET
      Š1=Š+1
      DO 7 I=K1,N
      DO 7 J=K1,N
 7    A(I,J)=A(I,J)-A(I,K)*A(K,J)/A(K,K)
      K=K1
      IF(K.LT.N) GOTO 8
      DET=DET*A(N,N)
      RETURN
      END
C
      SUBROUTINE DIM(NM,M)
      DIMENSION X(10),XD(10),XX(10),XXD(10),FI(10),F2(10)
      INTEGER G,N1,N2
      REAL T,XX,XXD,UI,U2,X,XD,F1,F2
      WRITE(*,50)
 50   FORMAT(15X,'€‘—…’ „ˆ•Ž’ŽŒˆ—…‘ŠŽ‰ ”“Š–ˆˆ')
      WRITE(*,*) '‚‚…„ˆ’… Ž‹œ˜Ž… —ˆ‘‹Ž "NG"'
      READ(*,*) G
      WRITE(*,*) '‚‚…„ˆ’… Œ€‹Ž… —ˆ‘‹Ž "’'"
      READ(*,*) ’
      WRITE(*,*) '‚‚…„ˆ’… ŠŽ‹-‚Ž €€Œ…’Ž‚ 1 ƒ“›'
      READ(*,*) N1
      WRITE(*,*) N1
      WRITE(*,*) '‚‚…„ˆ’… ŠŽ‹-‚Ž €€Œ…’Ž‚ 2 ƒ“›'
      READ(*,*) N2
      WRITE(*,*) N2
      WRITE(*,*) '‚‚Ž„ €€Œ…’Ž‚ 1 ƒ“›'
      WRITE(*,*)' '
      WRITE(*,*) '‚‚…„ˆ’… ŽŒˆ€‹œ›… ‡€—…ˆŸ €€Œ…’Ž‚'
      READ(*,*)(X(I),I=1,N1)
      WRITE(*,*)(X(I),I=1,NI)
      WRITE(*,*) '‚‚…„ˆ’… „Ž“‘’ˆŒ›… ‡€—…ˆŸ €€Œ…’Ž‚'
      READ(*,*)(XD(I),I=1,N1)
      WRITE(*,*)(XD(I),I=1,N1)
      WRITE(*,*) '‚‚Ž„ €€Œ…’Ž‚ 2 ƒ“›'
      WRITE(*,*)' '
      WRITE(*,*) '‚‚…„ˆ’… ŽŒˆ€‹œ›… ‡€—…ˆŸ €€Œ…’Ž‚'
      READ(*,*)(XX(I),I=1,N2)
      WRITE(*,*)(XX(I),I=1,N2)
      WRITE(*,*) '‚‚…„ˆ’… „Ž“‘’ˆŒ›… ‡€—…ˆŸ €€Œ…’Ž‚'
      READ(*,*) (XXD(I),I=1,N2)
      WRITE(*,*)(XXD(I),I=1,N2)
C
C €‘—…’ …‚Ž‰ ”“Š–ˆˆ F1
C
      DO 30 I=1,N1
      U1=(X(I)/(XD(I)*T))**G
      WRITE(*,*)'F1= ',U1
      WRITE(*,*) '‚‚…„ˆ’… ”“Š–ˆž ”(F1) 	'
      READ(*,*) F1(I)
 30   CONTINUE
C
C €‘—…’ ‚’ŽŽ‰ ”“Š–ˆˆ F2
C
      DO 40 I=1,N2
      U2=(XXD(I)/(XX(I)*T))**G
      WRITE(*,*)'F2= ',U2
      WRITE(*,*) '‚‚…„ˆ’… ”“Š–ˆž ”(F2) 	'
      READ(*,*) F2(I)
 40   CONTINUE
      P1=1.
      P2=1.
      DO 4 I=1,N1
      P1=P1*F1(I)
 4    CONTINUE
      DO 3 I=1,N2
      P2=P2*F2(I)
 3    CONTINUE
      DK=P1*P2
      WRITE(3,*) '„ˆ•Ž’ŽŒˆ—…‘Š€Ÿ ”“Š–ˆŸ = ',DK
      RETURN
      END