% LaTex
\documentstyle[12pt]{article}
%  #############################
\topmargin=-0.5cm
\textheight=23cm
\oddsidemargin=-0.1cm
\textwidth=16cm
%  ##############################
\pagestyle{plain}
%%%%%%%%%%
\newtheorem{th}{Theorem}
\newtheorem{prop}{Proposition}
%%%%%%%%%%
\newfont{\bg}{cmr12 scaled\magstep4}
\newcommand{\bigzerou}{\smash{\lower1.7ex\hbox{\bg 0}}}
\begin{document}

\title{�����Υ��������ݷ������Υա��ꥨ�����ο�����}
\author{ð�� ����\\
̾�Ų���Ω��� �ݽѹ�����
}

\maketitle


%%%%%%%%%%%%%%%%%%  abstract  %%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%   Section 1   %%%%%%%%%%%%%%%%%%%%
\section{Introduction}
��0ǯ���餤���˶���ε��ư˺��פ��� tarbo C ��Ȥä�
 degree 2 �� weight 12, weight 14 ��
 Siegel cusp form ��Fourier���������ץ������ cusp12.exe,
 cusp14.exe ������������礭�ʡ׷����ޤǵ����ǿ�p�˴ؤ���
 Hecke polynomial ��ʬ�򤹤뤳�Ȥ�ѻ�����ͽ�ۤ�Ω�Ƥ�줿
��weight 12 �ξ��p=2,3�� weight 14 ��p=2�ˡ�ɮ�ԤϤ⤦����
������䤷���ۤ����ɤ��ȴ������Τ� weight 12,p=5 ��
�������ͽ�ۤ�����Ω�Ĥ��Ȥ�Τ��᤿�����Τ���ˤ�
���ˤ�äȡ��礭�ʡ�Fourier���������ɬ�פ����ä���
�ϡ��ɥ������������եȥ��������ο���ˤ��ʸ���˺ܤ���
ñ��ʥץ������
�˰ͤäƤǤ��ǽ�Ǥ��ä������ƥ��ե� cusp12.exe,
 cusp14.exe �����˹��פ��줿�����̤ʥץ������
�Ǥ��뤬��ʬ������640�����Х���MS-DOS�Υץ������
�ʤΤ� weight 12,p=5 ��
���ˤϥ����С��ե�������̵���Ǥ��ä�������Ǥ����ø���
ʪ�ǤϤ��뤬ʸ���Υץ��������äƤߤ���
degree 2 �� weight 12, weight 14 ��
 Siegel cusp form �ϳơ����Ĥ����ʤ��Τ�
���줾��$F_{12}(Z),F_{14}(Z)$�Ƚ�([1])��
�ơ����Ĥ����ʤ��Ƥ�����ܤΰ㤤������櫓������
���ΰ㤤��ͭ�äƤ�ơ���Ʊ������ǽ񤯤��Ȥˤ��롣
���г���ʬ��a,b,c����¾����ʬ��e/2,f/2,g/2��
�����оι����[a,b,c,e,f,g]�Ƚ񤯤�$F(Z)=F_{12}(Z),F_{14}(Z)$
��
$$F(Z)=\sum _{a,b,c,e,f,g \in {\bf Z},[a,b,c,e,f,g]>0}
A([a,b,c,e,f,g])
\exp (2\pi i{\rm tr}([a,b,c,e,f,g]Z))$$ 
��ɽ����롣���ƥ��ե� cusp12.exe,
 cusp14.exe ��ư����ȣ��Ĥο������Ϥ������롣
a b c e f g �����Ϥ���Ȥ���A([a,b,c,e,f,g])�����Ϥ���롣
 ${\rm det}([a,b,c,e,f,g])\ge 81/2$�ޤǤ�A([a,b,c,e,f,g])��
ɽ��[1],[2]�ˤ��롣ɮ�Ԥ���᤿�Τ�
${\rm det}([a,b,c,e,f,g])\ge 5^4/2$��A([a,b,c,e,f,g])��
����ͽ�ۤ�Τ����Τ�ɬ�פʣ��Ĥ�A([a,b,c,e,f,g])����
������Ʊ���褦��a,b,c,e,f,g��A([a,b,c,e,f,g])�ϰ��
��ޤ�Τǡ����������Ω�Ĥ��⤷��ʤ���â���ץ������
��ɬ�פʽ�򤽤�����ľ���ʤ��Ȥ����ʤ���
\section{Result}
$\Gamma =\Gamma _3=Sp(3,{\bf Z}),r=1,2,3$�ΤȤ�
$I_r$��$r$��ñ�̹���Ȥ�
$$
T_r(p^3)=\Gamma \left( \begin{array}{cccc}
pI_r & & &\bigzerou\\
 & I_{3-r} & & \\
 & & pI_r & \\
 \bigzerou & & & p^2I_{3-r}
\end{array}\right) \Gamma
$$
�Ȥ��롣�ޤ�$F(Z)$���̾�Τ褦��$T_r(p^2)$����Ѥ�����
$F \vert T_r(p^2)(Z)$��FourierŸ����
\begin{eqnarray}
& & F \vert T_r(p^2)(Z)=\nonumber\\
& &\sum _{a,b,c,e,f,g \in {\bf Z},[a,b,c,e,f,g]>0}
A(T_r(p^2),[a,b,c,e,f,g])
\exp (2\pi i{\rm tr}([a,b,c,e,f,g]Z))\nonumber
\end{eqnarray}
�Ȥ��롣�ʲ� weight k ��k=12��$F$��$F=F_{12}$�Ȥ��롣
\begin{eqnarray}
&A_0=A([1,1,1,1,1,1]),
&A_1=A([1,1,13,-1,-1,0]),\nonumber\\
&A_2=A([1,1,17,1,1,1]),
&A_3=A([1,2,7,-2,-1,0]),\nonumber\\
&A_4=A([1,3,5,2,1,1]),
&A_6=A([2,3,3,1,2,2]),\nonumber\\
&A_{20}=A([1,19,19,13,1,1]),
&A_{27}=A([2,13,13,1,2,2]),\nonumber\\
&A_{35}=A([3,7,17,1,2,3]),
&A_{46}=A([5,5,15,-5,-5,0]),\nonumber\\
&A_{59}=A([6,9,9,-7,-4,-4])&\nonumber
\end{eqnarray}
�Ȥ����
\begin{prop}
\begin{eqnarray}
A(T_1(5^2),[1,1,1,1,1,1])&=&75\cdot 5^{3 k-12} A_0
+
5^{2 k-8} (
12 A_1-24 A_2-36 A_3+48 A_4-6 A_6)\nonumber\\
&+&
5^{k-4} 
(6 A_{20}+12 A_{35}+4 A_{59}+3 A_{27}+6 A_{46})
,\nonumber\\
A(T_2(5^2),[1,1,1,1,1,1])&=&-5^{3 k-12} A_0
+5^{2 k-8} 
(3 A_1+4 A_2+6 A_3+12 A_4+6 A_6).\nonumber
\end{eqnarray}
\end{prop}
�����롣$A_0=1$������������ȵ��ƥ��եȤ�
\begin{eqnarray}
& &A_1=-9826005,
A_2=1675705,
A_3=-17476295,\nonumber\\
& &A_4=16852395,
A_6=-28569550\nonumber
\end{eqnarray}
�Ϥ����狼��
ʸ���Υץ�����फ��
\begin{eqnarray}
& &A_{20}=-30307897267625,
A_{59}=-129587829571625,\nonumber\\
& &A_{27}=264736711948375,
A_{35}=25547315740375,\nonumber\\
& &A_{46}=266011961221750\nonumber
\end{eqnarray}
�Ȥʤ�Τ�$T_1(5^2),T_2(5^2)$�θ�ͭ�ͤϤ��줾��
$$ 1044613580296289062500, -14833396911621093750 $$
�Ǥ��뤳�Ȥ��狼�롣weight 12,weight 20 �Σ��ѿ��� 
normalized Hecke eigen cusp form �� n ����
��Fourier�����򤽤줾��$a_n,b_n$�Ȥ���ȡ�[1]�ˤ���褦��
$T_1(p^2),T_2(p^2)$�θ�ͭ�ͤ����줾��
\begin{eqnarray}
& &a_p^2 b_p p^8 (p+1)+
a_p^2 p^{16} (p^2-1)
+
b_p^2 p^{10}-
b_p p^{17} (p^2+1) (p+1)+
p^{26} (p^3+1) (p-1)
,\nonumber\\
& &a_p^2 p^{16}+b_p p^{17} (p+1)-p^{24} (p^3+1)\nonumber
\end{eqnarray}
�Ȥʤ�Ȥ����Τ�����ͽ�ۤʤΤ�
$$
a_5=4830,
b_5=-2377410
$$
���鼡�����롣
\begin{th}
p=5 �ΤȤ�����ͽ�ۤ���������
\end{th}
\section{Program}
��[1]��p.314��$F_12(Z)$��Fourier����A([1,19,19,13,1,1]�ˤ�
����ץ�������[2],[3],[4]�˽��äƺ�ä���Τ�ʸ����
��Ͽ�ˤ��λ���ˡ�򤳤��˺ܤ��롣������$E_8$-lAttice�θ�����ʬ
�������ˤʤ�褦���ܤ����ΤǳƷ�����2�ο����ܤˤʤäƤ��롣
Ʊ���褦��[a,b,c,e,f,g]�ξ�硢ɬ�פ��ѹ���ܤ���
A([a,b,c,e,f,g])�ϵ�ޤ롣ɬ�פ��ѹ��Ϥ��������Ҥ٤Ƥ��롣
���Υץ������ξ��
$$A([1,1,1,1,1,1])=13915694039040,$$
$$
A([1,19,19,13,1,1])=-421755425342925916078080000
$$
�ˤʤ�Τ�A([1,1,1,1,1,1])=1���������������
$$A([1,1,1,1,1,1])=-30307897267625$$
�Ȥʤ롣
\newline
��Pentium 200 �ᥬ�إ�ġ��� Dell Optiplex Gxmt5200�ˤ�
Windows95�ξ�Microsoft��Fortran Power Station��
���󥹥ȡ��뤵��Ƥ����ΤȤ��ƥץ������μ¹�ˡ��
���⤹�롣
\newline
��n20.for��double click����Power station��ư���롣
������project���뤫ʹ���Ƥ���Τ�yes�Ȥ��롣�����ȥ�С���
����Build��˥塼��Rebuild all�����򡣥���ѥ���ȥ��
�����顼�ʤ��˽�λ������Build��桼��Execute n20.exe������
�������˶��Υե�����n20.kek���äƤ������ޤ������ȥ�С���
����Tools��˥塼��Options���֥�˥塼��Work space������
Reloaded last project at startup��check��Ϥ����Ƥ������Ĥޤ�
��ư���뤿�ӿ�����project�ˤʤ�褦�ˤ��Ƥ�����n20.dat��
�Կ�16466400�����󥽡���˽��Ϥ����n20.exe������롣
�����ǰ�üMSDS(Microsoft Developer Studio)��λ���롣
\newline
��n203.for��double click����Power station��ư���롣
n20.for�ΤȤ���Ʊ�ͤ�n203.kek���äƤ������ʲ�MSDS��
��ư���뤿�ӥե�����̾.kek���ե�����̾.ke2���ե�����̾.ke3��
�ե�����̾.ke4�Ȥ����ե�������äƤ�����240�����󥽡���
�˽��Ϥ���� n203.exe ������롣n20.for, n203.for ��Ĺ�� 2, 38 ��
�������ƽ񤭽Ф���ΤʤΤ� 13, 39 ���ܤ� KON(38, KON(2 ��
 38, 2 �� [1,19,19,13,1,1] �Ǥʤ������б������ͤ�ľ��ɬ�פ����롣
��mf2033.for��double click����Power station��ư���롣
���Υե������NCASE=240��240��n203.exe�Υ��󥽡������
�ǡ�C(200000,8)��F(200000,8)��2000000��NCASE����礭�����
���Ǥ��ɤ�������ѥ���Ƿٹ�Ǥ뤬̵�롣5   240   240  ��
���󥽡���˽��Ϥ����mf2033.exe������롣
\newline
������ 5 �� mf204.for �� NCASE �ˤ���뤿�� mf204.for �� NCASE=5 ��
�ʤäƤ��롣C(200,8) �� 200 �� NCASE ����礭�����Ȥ������ȤǤ��롣
mf204.for��double click����Power station��ư���롣
\newline
����NCASE������μ�����ɤ����뤫�ǤäƤ���ʪ�ˤĤ��Ƥ�
[1,19,19,13,1,1]�Ǥʤ�¾�ξ�硢���ƽ�ľ��ɬ�פ�����
���λ����ο���ɤ����뤫����ꤹ���ΤǤ��롣��
\newline
��mf205.for��double click����Power station��ư���롣
mf205.for�Ǥ�n20.exe�Υ��󥽡�����Ϥ���NCASE=16466400��
mf2033.exe�Υ��󥽡�����Ϥ���MCASE=5����Ƥ��롣
24���� .EQ. 4 ��4��[1,19,19,13,1,1]�θ��1��4�ܤǤ��롣
C(16466400,8)��16466400��n20.dat�ιԿ��Ǥ��ꡢn20.exe��
���󥽡�����ϤǤ��롣
\newline
��mf2061.for��double click����Power station��ư���롣
C(882280,8)��F(882280,8)��882280��mf205.kek����2��κǽ��ԡ�
NCASE=176456��mf205.kek����1�󤫤��ޤ롣�夫��9����NCASE*5
��5��夫��21���ܤ�5��mf2033.exe�Υ��󥽡�����Ϥ�5�Ǥ��롣
A(100,3)��100��5����礭�����Ȥ������Ȥ�I(5000,14)��
5000��mf204.kek����2��κ�����4095�򲼤ޤ��ʤ����Ȥ���
���Ƥ��롣��������6167��������882280�������������󥽡������
�����mf2061.exe�Ͻ���롣�ʣ��ݣ����֤����롣��
\newline
��n207.for��double click����Power station��ư���롣
��괺����Rebuild all��¹Ԥ���ȥ��顼���Ǥ롣�虜��
�������Ƥ��饿���ȥ�С��β���Insert��˥塼��
Files into project������gg.for, sa.for, wa.for, wss.for,
 syoukei.for��project�˲ä��롣�Ƥ�Rebuild all��¹ԡ�
����n207.kek��˺�줺�˺�äƤ�����n207.exe��¹ԡ�
��λ�ˤϣ��ݣ�����ɬ�פǤ��롣�ġ���С��β������Debug 
option�Ǥʤ�Release option�ˤ��Ƥ��������������餫®����
1-30�Ԥ�176456��mf205.kek���裱��6167��mf2061.kek���裱��
�κǽ��ԡ������ֹ��5��������Ƥ�5��mf2033.exe��
���󥽡�����Ϥ�5�Ǥ���Τ�¾�ξ��ˤ�����ľ��ɬ�פ����롣
�裶���Ԥ�52��4*13��[1,19,19,13,1,1]��13������Ƥ��롣
        -1      0       0       421     7554    2534
        2925    9160    0
���󥽡���˽��Ϥ����n207.exe������롣�����
������$$-421755425342925916078080000$$�Ȥʤ뤳�Ȥ��̣���롣
\newline
��mm207.for��n207.for��Z11,Z12,Z21,Z22,Z31,Z32,Z41,Z42��
1000�ʲ��Ǥ��뤳�Ȥ�Τ����Τ˻Ȥ��������ʤäƤʤ�����
���������ΤǤʤ����⤷��ʤ���
\begin{thebibliography}{99}
\bibitem{B-1}\ {\sc I.Miyawaki},{Numerical examples of Siegel cusp 
forms
 of degree 3 and their zeta-functions},{\sl Mem.Fac.Sci.Kyushu Univ.,
Ser.A},{\bf 46}\ (1992),\ 307--339.
\bibitem{B-2}\ { ���ư˺���},{���� Siegel cusp form ��zeta�ؿ��μ����
�Ĥ���},
{ �������ֵ�Ͽ�ʡ�)}

\bibitem{B-3}\ { ���ư˺���},��Unpublished work��
Part 1. Theta series attached to $E_8$ -lattice,
Part 2. Theta series with spherical-function.

\bibitem{B-4}\ { ���ư˺���},��cusp14.exe�Υ������ե������.

\end{thebibliography}

\begin{tiny}
\begin{verbatim}
n20.for

      PROGRAM TEST
      INTEGER I1,I2,I3,I4,I5,I6,I7,I8,A0
      INTEGER KON2
      OPEN (6,FILE='N20.dat',STATUS='NEW')
      A0=0
      DO 80000 I1=-6,6
      DO 75000 I2=-6,6
      DO 70000 I3=-6,6
      DO 60000 I4=-6,6
      DO 50000 I5=-6,6
      DO 40000 I6=-6,6
      DO 30000 I7=-6,6
      I8=KON2(38-I1**2-I2**2-I3**2-I4**2-I5**2-I6**2-I7**2)
      IF ((I8 .NE. -1) .AND. 
     .(MOD(I1+I2+I3+I4+I5+I6+I7+I8,2) .EQ. 0)) THEN
      A0=A0+1
      WRITE(6,90000) 2*I1,2*I2,2*I3,2*I4,2*I5,2*I6,2*I7,2*I8
      ENDIF
      IF ((I8 .GT. 0) .AND. 
     .(MOD(I1+I2+I3+I4+I5+I6+I7-I8,2) .EQ. 0)) THEN
      A0=A0+1
      WRITE(6,90000) 2*I1,2*I2,2*I3,2*I4,2*I5,2*I6,2*I7,-2*I8
      ENDIF
30000 CONTINUE
40000 CONTINUE
50000 CONTINUE
60000 CONTINUE
70000 CONTINUE
75000 CONTINUE
80000 CONTINUE
90000 FORMAT(8I4)
      DO 98000 I1=-11,11,2
      DO 97000 I2=-11,11,2
      DO 96000 I3=-11,11,2
      DO 95000 I4=-11,11,2
      DO 94000 I5=-11,11,2
      DO 93000 I6=-11,11,2
      DO 92000 I7=-11,11,2
      I8=KON2(38*4-I1**2-I2**2-I3**2-I4**2-I5**2-I6**2-I7**2)
      IF ((I8 .NE. -1) .AND. (MOD(I8,2) .EQ. 1) .AND. 
     .(MOD(I1+I2+I3+I4+I5+I6+I7+I8,4) .EQ. 0)) THEN
      A0=A0+1
      WRITE(6,90000) I1,I2,I3,I4,I5,I6,I7,I8
      ENDIF
      IF ((I8 .GT. 0) .AND. (MOD(I8,2) .EQ. 1) .AND. 
     .(MOD(I1+I2+I3+I4+I5+I6+I7-I8,4) .EQ. 0)) THEN
      A0=A0+1
      WRITE(6,90000) I1,I2,I3,I4,I5,I6,I7,-I8
      ENDIF
92000 CONTINUE
93000 CONTINUE
94000 CONTINUE
95000 CONTINUE
96000 CONTINUE
97000 CONTINUE
98000 CONTINUE
      WRITE(*,*) A0
      OPEN (7,FILE='N20.KEK',STATUS='OLD')
      WRITE(7,*) A0
      CLOSE (7)
      CLOSE (6)
      END
      INTEGER FUNCTION KON2(BB0)
      INTEGER BB0
      IF (BB0 .LT. 0) THEN
         KON2=-1
        ELSEIF (BB0 .EQ. 0) THEN 
           KON2=0
           RETURN
          ELSEIF (BB0 .EQ. 1) THEN 
             KON2=1
             RETURN
            ELSEIF (BB0 .EQ. 4) THEN 
               KON2=2
               RETURN
              ELSEIF (BB0 .EQ. 9) THEN 
                 KON2=3
                 RETURN
                ELSEIF (BB0 .EQ. 16) THEN 
                   KON2=4
                   RETURN
                  ELSEIF (BB0 .EQ. 25) THEN 
                     KON2=5
                     RETURN
                    ELSEIF (BB0 .EQ. 36) THEN 
                       KON2=6
                       RETURN
                      ELSEIF (BB0 .EQ. 49) THEN 
                         KON2=7
                         RETURN
                        ELSEIF (BB0 .EQ. 64) THEN 
                           KON2=8
                           RETURN
                          ELSEIF (BB0 .EQ. 81) THEN 
                             KON2=9
                             RETURN
                            ELSEIF (BB0 .EQ. 100) THEN 
                               KON2=10
                               RETURN
                              ELSEIF (BB0 .EQ. 121) THEN 
                                 KON2=11
                                 RETURN
      ELSE 
         KON2=-1
         RETURN
      ENDIF
      END

\end{verbatim}
\begin{verbatim}

n203.for

      PROGRAM TEST
      INTEGER I1,I2,I3,I4,I5,I6,I7,I8,A0
      INTEGER KON2
      OPEN (6,FILE='N203.dat',STATUS='NEW')
      A0=0
      DO 80000 I1=-4,4
      DO 75000 I2=-4,4
      DO 70000 I3=-4,4
      DO 60000 I4=-4,4
      DO 50000 I5=-4,4
      DO 40000 I6=-4,4
      DO 30000 I7=-4,4
      I8=KON2(2-I1**2-I2**2-I3**2-I4**2-I5**2-I6**2-I7**2)
      IF ((I8 .NE. -1) .AND. 
     .(MOD(I1+I2+I3+I4+I5+I6+I7+I8,2) .EQ. 0)) THEN
      A0=A0+1
      WRITE(6,90000) 2*I1,2*I2,2*I3,2*I4,2*I5,2*I6,2*I7,2*I8
      ENDIF
      IF ((I8 .GT. 0) .AND. 
     .(MOD(I1+I2+I3+I4+I5+I6+I7-I8,2) .EQ. 0)) THEN
      A0=A0+1
      WRITE(6,90000) 2*I1,2*I2,2*I3,2*I4,2*I5,2*I6,2*I7,-2*I8
      ENDIF
30000 CONTINUE
40000 CONTINUE
50000 CONTINUE
60000 CONTINUE
70000 CONTINUE
75000 CONTINUE
80000 CONTINUE
90000 FORMAT(8I4)
      DO 98000 I1=-7,7,2
      DO 97000 I2=-7,7,2
      DO 96000 I3=-7,7,2
      DO 95000 I4=-7,7,2
      DO 94000 I5=-7,7,2
      DO 93000 I6=-7,7,2
      DO 92000 I7=-7,7,2
      I8=KON2(2*4-I1**2-I2**2-I3**2-I4**2-I5**2-I6**2-I7**2)
      IF ((I8 .NE. -1) .AND. (MOD(I8,2) .EQ. 1) .AND. 
     .(MOD(I1+I2+I3+I4+I5+I6+I7+I8,4) .EQ. 0)) THEN
      A0=A0+1
      WRITE(6,90000) I1,I2,I3,I4,I5,I6,I7,I8
      ENDIF
      IF ((I8 .GT. 0) .AND. (MOD(I8,2) .EQ. 1) .AND. 
     .(MOD(I1+I2+I3+I4+I5+I6+I7-I8,4) .EQ. 0)) THEN
      A0=A0+1
      WRITE(6,90000) I1,I2,I3,I4,I5,I6,I7,-I8
      ENDIF
92000 CONTINUE
93000 CONTINUE
94000 CONTINUE
95000 CONTINUE
96000 CONTINUE
97000 CONTINUE
98000 CONTINUE
      WRITE(*,*) A0
      OPEN (7,FILE='N203.KEK',STATUS='OLD')
      WRITE(7,*) A0
      CLOSE (7)
      CLOSE (6)
      END
      INTEGER FUNCTION KON2(BB0)
      INTEGER BB0
      IF (BB0 .LT. 0) THEN
         KON2=-1
        ELSEIF (BB0 .EQ. 0) THEN 
           KON2=0
           RETURN
          ELSEIF (BB0 .EQ. 1) THEN 
             KON2=1
             RETURN
            ELSEIF (BB0 .EQ. 4) THEN 
               KON2=2
               RETURN
              ELSEIF (BB0 .EQ. 9) THEN 
                 KON2=3
                 RETURN
                ELSEIF (BB0 .EQ. 16) THEN 
                   KON2=4
                   RETURN
                  ELSEIF (BB0 .EQ. 25) THEN 
                     KON2=5
                     RETURN
                    ELSEIF (BB0 .EQ. 36) THEN 
                       KON2=6
                       RETURN
                      ELSEIF (BB0 .EQ. 49) THEN 
                         KON2=7
                         RETURN
                        ELSEIF (BB0 .EQ. 64) THEN 
                           KON2=8
                           RETURN
                          ELSEIF (BB0 .EQ. 81) THEN 
                             KON2=9
                             RETURN
                            ELSEIF (BB0 .EQ. 100) THEN 
                               KON2=10
                               RETURN
                              ELSEIF (BB0 .EQ. 121) THEN 
                                 KON2=11
                                 RETURN
      ELSE 
         KON2=-1
         RETURN
      ENDIF
      END

\end{verbatim}
\begin{verbatim}

mf2033.for

      PROGRAM TEST
      INTEGER I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13,A0
      INTEGER C1,C2,C3,C4,C5,C6,C7,C8,C9,C0
      INTEGER X1,X2,X3,X4,X5,X6,X7,X8,X9
      INTEGER K,J,A1,B0,L,M,L0,N,M0,N0
      INTEGER C(200000,8),NCASE,MCASE
      INTEGER B(2**14,8),F(200000,8)
      NCASE=240
      OPEN (5,FILE='N203.DAT',STATUS='OLD')
      READ(5,*) ((C(K,J),J=1,8),K=1,NCASE)
      CLOSE (5)
      M0=0
      A1=1+0
      OPEN (6,FILE='MF2033.KEK',STATUS='OLD')
      OPEN (7,FILE='MF2033.KE2',STATUS='OLD')
      OPEN (8,FILE='MF2033.KE3',STATUS='OLD')
      OPEN (9,FILE='MF2033.KE4',STATUS='OLD')
      M0=0
      B0=0
      DO 18500 K=1,NCASE
      DO 10020 M=1,M0
      IF ((C(K,1) .EQ. F(M,1)) .AND. (C(K,2) .EQ. F(M,2))
     . .AND. (C(K,3) .EQ. F(M,3)) .AND. (C(K,4) .EQ. F(M,4))
     . .AND. (C(K,5) .EQ. F(M,5)) .AND. (C(K,6) .EQ. F(M,6))
     . .AND. (C(K,7) .EQ. F(M,7)) .AND. (C(K,8) .EQ. F(M,8))) 
     . GOTO 18500
10020 CONTINUE
      DO 18206 I13=0,1
      DO 18205 I12=0,1
      DO 18204 I11=0,1
      DO 18203 I10=0,1
      DO 18200 I0=0,1
      DO 18100 I9=0,1
      DO 18000 I8=0,1
      DO 17000 I7=0,1
      DO 16000 I6=0,1
      DO 15000 I5=0,1
      DO 14000 I4=0,1
      DO 13000 I3=0,1
      DO 12000 I2=0,1
      DO 11000 I1=0,1
      X1=C(K,1)
      X2=C(K,2)
      X3=C(K,3)
      X4=C(K,4)
      X5=C(K,5)
      X6=C(K,6)
      X7=C(K,7)
      X8=C(K,8)
      CALL FF(I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13,
     .X1,X2,X3,X4,X5,X6,X7,X8)
      A1=1+I0*2**9+I9*2**8+I8*2**7+I7*2**6+I6*2**5+I5*2**4+
     .I4*2**3+I3*2**2+I2*2+I1+I10*2**10+I11*2**11+I12*2**12+I13*2**13
      B(A1,1)=X1
      B(A1,2)=X2
      B(A1,3)=X3
      B(A1,4)=X4
      B(A1,5)=X5
      B(A1,6)=X6
      B(A1,7)=X7
      B(A1,8)=X8
11000 CONTINUE
12000 CONTINUE
13000 CONTINUE
14000 CONTINUE
15000 CONTINUE
16000 CONTINUE
17000 CONTINUE
18000 CONTINUE
18100 CONTINUE
18200 CONTINUE
18203 CONTINUE
18204 CONTINUE
18205 CONTINUE
18206 CONTINUE
      L0=0
      DO 18300 L=1,2**14
      DO 18250 M=1,L-1
      IF ((B(L,1) .EQ. B(M,1)) .AND. (B(L,2) .EQ. B(M,2))
     . .AND. (B(L,3) .EQ. B(M,3)) .AND. (B(L,4) .EQ. B(M,4))
     . .AND. (B(L,5) .EQ. B(M,5)) .AND. (B(L,6) .EQ. B(M,6))
     . .AND. (B(L,7) .EQ. B(M,7)) .AND. (B(L,8) .EQ. B(M,8))) 
     . GOTO 18300
18250 CONTINUE
      L0=L0+1
      M0=M0+1
      F(M0,1)=B(L,1)
      F(M0,2)=B(L,2)
      F(M0,3)=B(L,3)
      F(M0,4)=B(L,4)
      F(M0,5)=B(L,5)
      F(M0,6)=B(L,6)
      F(M0,7)=B(L,7)
      F(M0,8)=B(L,8)
18300 CONTINUE
      B0=B0+1
      WRITE(7,90000) C(K,1),C(K,2),C(K,3),C(K,4)
     .,C(K,5),C(K,6),C(K,7),C(K,8)
      WRITE(8,*) B0,K,L0,M0
18500 CONTINUE
      WRITE(6,*) B0,M0
      DO 20000 L=1,M0
      WRITE(9,90000) F(L,1),F(L,2),F(L,3),F(L,4)
     .,F(L,5),F(L,6),F(L,7),F(L,8)
20000 CONTINUE
      CLOSE(6)
      CLOSE(8)
      CLOSE(7)
      CLOSE(9)
      WRITE(*,*) B0,M0,NCASE
90000 FORMAT(8I4)
      END
      SUBROUTINE FF(I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13,
     .X1,X2,X3,X4,X5,X6,X7,X8)
      INTEGER I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13
      INTEGER X1,X2,X3,X4,X5,X6,X7,X8,X9
      INTEGER Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9
      INTEGER K,J,A1
      DO 91000 K=1,I1
      Y1=X1
      X1=X2
      X2=Y1
91000 CONTINUE
      DO 92000 K=1,I2
      Y4=X4
      X4=X3
      X3=Y4
92000 CONTINUE
      DO 93000 K=1,I3
      Y5=X5
      X5=X6
      X6=Y5
93000 CONTINUE
      DO 93100 K=1,I4
      Y7=X7
      X7=X8
      X8=Y7
93100 CONTINUE
      DO 93200 K=1,I5
      Y3=X3
      X3=X1
      X1=Y3
      Y4=X4
      X4=X2
      X2=Y4
93200 CONTINUE
      DO 94000 K=1,I6
      Y5=X5
      X5=X3
      X3=Y5
      Y6=X6
      X6=X4
      X4=Y6
94000 CONTINUE
      DO 94100 K=1,I7
      Y5=X5
      X5=X1
      X1=Y5
      Y6=X6
      X6=X2
      X2=Y6
94100 CONTINUE
      DO 94150 K=1,I8
      X1=-X1
      X8=-X8
94150 CONTINUE
      DO 95000 K=1,I9
      X2=-X2
      X8=-X8
95000 CONTINUE
      DO 95100 K=1,I0
      X3=-X3
      X8=-X8
95100 CONTINUE
      DO 95150 K=1,I10
      X4=-X4
      X8=-X8
95150 CONTINUE
      DO 95200 K=1,I11
      X5=-X5
      X8=-X8
95200 CONTINUE
      DO 95500 K=1,I12
      X6=-X6
      X8=-X8
95500 CONTINUE
      DO 95600 K=1,I13
      X7=-X7
      X8=-X8
95600 CONTINUE
      RETURN
      END

\end{verbatim}
\begin{verbatim}

mf204.for

      PROGRAM TEST
      INTEGER I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13
      INTEGER X1,X2,X3,X4,X5,X6,X7,X8,X9
      INTEGER K,J,A1,B0,L,M,L0,N,M0,N0
      INTEGER C(200,8),NCASE,MCASE
      NCASE=5
      OPEN (5,FILE='MF2033.KE2',STATUS='OLD')
      READ(5,*) ((C(K,J),J=1,8),K=1,NCASE)
      CLOSE (5)
      OPEN (6,FILE='MF204.KEK',STATUS='OLD')
      OPEN (7,FILE='MF204.KE2',STATUS='OLD')
      B0=0
      DO 18500 K=1,NCASE
      L0=0
10020 CONTINUE
      DO 18206 I13=0,1
      DO 18205 I12=0,1
      DO 18204 I11=0,1
      DO 18203 I10=0,1
      DO 18200 I0=0,1
      DO 18100 I9=0,1
      DO 18000 I8=0,1
      DO 17000 I7=0,1
      DO 16000 I6=0,1
      DO 15000 I5=0,1
      DO 14000 I4=0,1
      DO 13000 I3=0,1
      DO 12000 I2=0,1
      DO 11000 I1=0,1
      X1=C(K,1)
      X2=C(K,2)
      X3=C(K,3)
      X4=C(K,4)
      X5=C(K,5)
      X6=C(K,6)
      X7=C(K,7)
      X8=C(K,8)
      CALL FF(I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13,
     .X1,X2,X3,X4,X5,X6,X7,X8)
      IF ((X1 .EQ. C(K,1)) .AND. (X2 .EQ. C(K,2))
     . .AND. (X3 .EQ. C(K,3)) .AND. (X4 .EQ. C(K,4))
     . .AND. (X5 .EQ. C(K,5)) .AND. (X6 .EQ. C(K,6))
     . .AND. (X7 .EQ. C(K,7)) .AND. (X8 .EQ. C(K,8))) THEN
      L0=L0+1
      WRITE(7,80000) I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13
      ENDIF
11000 CONTINUE
12000 CONTINUE
13000 CONTINUE
14000 CONTINUE
15000 CONTINUE
16000 CONTINUE
17000 CONTINUE
18000 CONTINUE
18100 CONTINUE
18200 CONTINUE
18203 CONTINUE
18204 CONTINUE
18205 CONTINUE
18206 CONTINUE
      B0=B0+L0
      WRITE(6,*) K,L0,B0
18500 CONTINUE
      CLOSE(6)
      CLOSE(7)
80000 FORMAT(14I4)
90000 FORMAT(8I4)
      END
      SUBROUTINE FF(I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13,
     .X1,X2,X3,X4,X5,X6,X7,X8)
      INTEGER I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13
      INTEGER X1,X2,X3,X4,X5,X6,X7,X8,X9
      INTEGER Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9
      INTEGER K,J,A1
      DO 91000 K=1,I1
      Y1=X1
      X1=X2
      X2=Y1
91000 CONTINUE
      DO 92000 K=1,I2
      Y4=X4
      X4=X3
      X3=Y4
92000 CONTINUE
      DO 93000 K=1,I3
      Y5=X5
      X5=X6
      X6=Y5
93000 CONTINUE
      DO 93100 K=1,I4
      Y7=X7
      X7=X8
      X8=Y7
93100 CONTINUE
      DO 93200 K=1,I5
      Y3=X3
      X3=X1
      X1=Y3
      Y4=X4
      X4=X2
      X2=Y4
93200 CONTINUE
      DO 94000 K=1,I6
      Y5=X5
      X5=X3
      X3=Y5
      Y6=X6
      X6=X4
      X4=Y6
94000 CONTINUE
      DO 94100 K=1,I7
      Y5=X5
      X5=X1
      X1=Y5
      Y6=X6
      X6=X2
      X2=Y6
94100 CONTINUE
      DO 94150 K=1,I8
      X1=-X1
      X8=-X8
94150 CONTINUE
      DO 95000 K=1,I9
      X2=-X2
      X8=-X8
95000 CONTINUE
      DO 95100 K=1,I0
      X3=-X3
      X8=-X8
95100 CONTINUE
      DO 95150 K=1,I10
      X4=-X4
      X8=-X8
95150 CONTINUE
      DO 95200 K=1,I11
      X5=-X5
      X8=-X8
95200 CONTINUE
      DO 95500 K=1,I12
      X6=-X6
      X8=-X8
95500 CONTINUE
      DO 95600 K=1,I13
      X7=-X7
      X8=-X8
95600 CONTINUE
      RETURN
      END

\end{verbatim}
\begin{verbatim}

mf205.for

      PROGRAM TEST
      INTEGER I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13,A0
      INTEGER C1,C2,C3,C4,C5,C6,C7,C8,C9,C0
      INTEGER X1,X2,X3,X4,X5,X6,X7,X8,X9
      INTEGER K,J,A1,B0,L,M,L0,N,M0,N0
      INTEGER C(1646400,8),NCASE,MCASE
      INTEGER A(8)
      NCASE=1646400
      MCASE=5
      OPEN (5,FILE='N20.DAT',STATUS='OLD')
      READ(5,*) ((C(K,J),J=1,8),K=1,NCASE)
      CLOSE (5)
      OPEN (8,FILE='MF2033.KE2',STATUS='OLD')
      OPEN (7,FILE='MF205.KE2',STATUS='OLD')
      OPEN (6,FILE='MF205.KEK',STATUS='OLD')
      B0=0
      DO 18500 L=1,MCASE
      READ(8,*) (A(J),J=1,8)
      L0=0
      DO 10000 K=1,NCASE
      IF ((C(K,1)*A(1)+C(K,2)*A(2))
     .+(C(K,3)*A(3)+C(K,4)*A(4))
     .+(C(K,5)*A(5)+C(K,6)*A(6))
     .+(C(K,7)*A(7)+C(K,8)*A(8)) .EQ. 4) THEN 
      L0=L0+1
      WRITE(7,90000) C(K,1),C(K,2),C(K,3),C(K,4)
     .,C(K,5),C(K,6),C(K,7),C(K,8)
      ENDIF
10000 CONTINUE
      B0=B0+L0
      WRITE(6,*) L0,B0
18500 CONTINUE
      CLOSE(6)
      CLOSE(8)
      CLOSE(7)
80000 FORMAT(14I4)
90000 FORMAT(8I4)
      END
      SUBROUTINE FF(I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13,
     .X1,X2,X3,X4,X5,X6,X7,X8)
      INTEGER I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13
      INTEGER X1,X2,X3,X4,X5,X6,X7,X8,X9
      INTEGER Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9
      INTEGER K,J,A1
      DO 91000 K=1,I1
      Y1=X1
      X1=X2
      X2=Y1
91000 CONTINUE
      DO 92000 K=1,I2
      Y4=X4
      X4=X3
      X3=Y4
92000 CONTINUE
      DO 93000 K=1,I3
      Y5=X5
      X5=X6
      X6=Y5
93000 CONTINUE
      DO 93100 K=1,I4
      Y7=X7
      X7=X8
      X8=Y7
93100 CONTINUE
      DO 93200 K=1,I5
      Y3=X3
      X3=X1
      X1=Y3
      Y4=X4
      X4=X2
      X2=Y4
93200 CONTINUE
      DO 94000 K=1,I6
      Y5=X5
      X5=X3
      X3=Y5
      Y6=X6
      X6=X4
      X4=Y6
94000 CONTINUE
      DO 94100 K=1,I7
      Y5=X5
      X5=X1
      X1=Y5
      Y6=X6
      X6=X2
      X2=Y6
94100 CONTINUE
      DO 94150 K=1,I8
      X1=-X1
      X8=-X8
94150 CONTINUE
      DO 95000 K=1,I9
      X2=-X2
      X8=-X8
95000 CONTINUE
      DO 95100 K=1,I0
      X3=-X3
      X8=-X8
95100 CONTINUE
      DO 95150 K=1,I10
      X4=-X4
      X8=-X8
95150 CONTINUE
      DO 95200 K=1,I11
      X5=-X5
      X8=-X8
95200 CONTINUE
      DO 95500 K=1,I12
      X6=-X6
      X8=-X8
95500 CONTINUE
      DO 95600 K=1,I13
      X7=-X7
      X8=-X8
95600 CONTINUE
      RETURN
      END

\end{verbatim}
\begin{verbatim}

mf2061.for

      PROGRAM TEST
      INTEGER X1,X2,X3,X4,X5,X6,X7,X8,X9
      INTEGER K,J,A0,A1,B0,L,M,L0,N,M0,N0
      INTEGER C(882280,8),NCASE,MCASE,A(100,3)
      INTEGER B(2**14,8),F(882280,8),I(5000,14)
C      NCASE=
      NCASE=176456
      OPEN (5,FILE='MF205.KE2',STATUS='OLD')
      READ(5,*) ((C(K,J),J=1,8),K=1,NCASE*5)
      CLOSE (5)
      OPEN (5,FILE='MF204.KEK',STATUS='OLD')
      READ(5,*) ((A(K,J),J=1,3),K=1,5)
      CLOSE (5)
      OPEN (6,FILE='MF2061.KEK',STATUS='OLD')
      OPEN (7,FILE='MF2061.KE2',STATUS='OLD')
      OPEN (8,FILE='MF2061.KE3',STATUS='OLD')
      OPEN (9,FILE='MF2061.KE4',STATUS='OLD')
      OPEN (5,FILE='MF204.KE2',STATUS='OLD')
      A0=0
      N0=0
      DO 30000 J=1,5
      READ(5,*) ((I(L,K),K=1,14),L=1,A(J,2))
C      READ(5,*) I(L,1),I(L,2),I(L,3),I(L,4),I(L,5),I(L,6),
C       .I(L,7),I(L,8),I(L,9),I(L,10),I(L,11),I(L,12),I(L,13),I(L,14)
      M0=0
      B0=0
      DO 18500 K=NCASE*(J-1)+1,NCASE*J
      DO 10020 M=1,M0
      IF ((C(K,1) .EQ. F(M,1)) .AND. (C(K,2) .EQ. F(M,2))
     . .AND. (C(K,3) .EQ. F(M,3)) .AND. (C(K,4) .EQ. F(M,4))
     . .AND. (C(K,5) .EQ. F(M,5)) .AND. (C(K,6) .EQ. F(M,6))
     . .AND. (C(K,7) .EQ. F(M,7)) .AND. (C(K,8) .EQ. F(M,8))) 
     . GOTO 18500
10020 CONTINUE
      A1=0
C      DO 11000 L=1,A(J,2)
      DO 11000 L=1,A(J,2)
      X1=C(K,1)
      X2=C(K,2)
      X3=C(K,3)
      X4=C(K,4)
      X5=C(K,5)
      X6=C(K,6)
      X7=C(K,7)
      X8=C(K,8)
      CALL FF(I(L,1),I(L,2),I(L,3),I(L,4),I(L,5),I(L,6),
     .I(L,7),I(L,8),I(L,9),I(L,10),I(L,11),I(L,12),I(L,13),I(L,14),
     .X1,X2,X3,X4,X5,X6,X7,X8)
      A1=A1+1
      B(A1,1)=X1
      B(A1,2)=X2
      B(A1,3)=X3
      B(A1,4)=X4
      B(A1,5)=X5
      B(A1,6)=X6
      B(A1,7)=X7
      B(A1,8)=X8
11000 CONTINUE
      L0=0
C      DO 18300 L=1,A(J,2)
      DO 18300 L=1,A(J,2)
      DO 18250 M=1,L-1
      IF ((B(L,1) .EQ. B(M,1)) .AND. (B(L,2) .EQ. B(M,2))
     . .AND. (B(L,3) .EQ. B(M,3)) .AND. (B(L,4) .EQ. B(M,4))
     . .AND. (B(L,5) .EQ. B(M,5)) .AND. (B(L,6) .EQ. B(M,6))
     . .AND. (B(L,7) .EQ. B(M,7)) .AND. (B(L,8) .EQ. B(M,8))) 
     . GOTO 18300
18250 CONTINUE
      L0=L0+1
      M0=M0+1
      F(M0,1)=B(L,1)
      F(M0,2)=B(L,2)
      F(M0,3)=B(L,3)
      F(M0,4)=B(L,4)
      F(M0,5)=B(L,5)
      F(M0,6)=B(L,6)
      F(M0,7)=B(L,7)
      F(M0,8)=B(L,8)
18300 CONTINUE
      B0=B0+1
      WRITE(7,90000) C(K,1),C(K,2),C(K,3),C(K,4)
     .,C(K,5),C(K,6),C(K,7),C(K,8)
      WRITE(8,*) B0,K,L0,M0
18500 CONTINUE
      WRITE(6,*) B0,M0
      DO 20000 L=1,M0
      WRITE(9,90000) F(L,1),F(L,2),F(L,3),F(L,4)
     .,F(L,5),F(L,6),F(L,7),F(L,8)
20000 CONTINUE
      A0=A0+B0
      N0=N0+M0
30000 CONTINUE
      WRITE(6,*) A0,N0
      CLOSE(5)
      CLOSE(6)
      CLOSE(8)
      CLOSE(7)
      CLOSE(9)
      WRITE(*,*) A0,N0
80000 FORMAT(14I4)
90000 FORMAT(8I4)
      END
      SUBROUTINE FF(I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13,
     .X1,X2,X3,X4,X5,X6,X7,X8)
      INTEGER I1,I2,I3,I4,I5,I6,I7,I8,I9,I0,I10,I11,I12,I13
      INTEGER X1,X2,X3,X4,X5,X6,X7,X8,X9
      INTEGER Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9
      INTEGER K,J,A1
      DO 91000 K=1,I1
      Y1=X1
      X1=X2
      X2=Y1
91000 CONTINUE
      DO 92000 K=1,I2
      Y4=X4
      X4=X3
      X3=Y4
92000 CONTINUE
      DO 93000 K=1,I3
      Y5=X5
      X5=X6
      X6=Y5
93000 CONTINUE
      DO 93100 K=1,I4
      Y7=X7
      X7=X8
      X8=Y7
93100 CONTINUE
      DO 93200 K=1,I5
      Y3=X3
      X3=X1
      X1=Y3
      Y4=X4
      X4=X2
      X2=Y4
93200 CONTINUE
      DO 94000 K=1,I6
      Y5=X5
      X5=X3
      X3=Y5
      Y6=X6
      X6=X4
      X4=Y6
94000 CONTINUE
      DO 94100 K=1,I7
      Y5=X5
      X5=X1
      X1=Y5
      Y6=X6
      X6=X2
      X2=Y6
94100 CONTINUE
      DO 94150 K=1,I8
      X1=-X1
      X8=-X8
94150 CONTINUE
      DO 95000 K=1,I9
      X2=-X2
      X8=-X8
95000 CONTINUE
      DO 95100 K=1,I0
      X3=-X3
      X8=-X8
95100 CONTINUE
      DO 95150 K=1,I10
      X4=-X4
      X8=-X8
95150 CONTINUE
      DO 95200 K=1,I11
      X5=-X5
      X8=-X8
95200 CONTINUE
      DO 95500 K=1,I12
      X6=-X6
      X8=-X8
95500 CONTINUE
      DO 95600 K=1,I13
      X7=-X7
      X8=-X8
95600 CONTINUE
      RETURN
      END

\end{verbatim}
\begin{verbatim}

n207.for

      PROGRAM TEST
      INTEGER J,K,I,BB,L
      INTEGER Z11,Z12,Z21,Z22,Z31,Z32,Z41,Z42
      INTEGER DS,D8,D7,D6,D5,D4,D3,D2,D1,D0
      INTEGER CS,C8,C7,C6,C5,C4,C3,C2,C1,C0
      INTEGER BS,B8,B7,B6,B5,B4,B3,B2,B1,B0
      INTEGER C(176456*5,8),B(6167,8),D(6167,4),E(5,4),MCASE(5,2)
      INTEGER LCASE(5),A(5,8)
      NCASE=176456
      OPEN (5,FILE='MF2033.KE2',STATUS='OLD')
      READ(5,*) ((A(K,J),J=1,8),K=1,5)
      CLOSE (5)
      OPEN (5,FILE='MF205.KE2',STATUS='OLD')
      READ(5,*) ((C(K,J),J=1,8),K=1,NCASE*5)
      CLOSE (5)
      OPEN (5,FILE='MF2061.KE2',STATUS='OLD')
      READ(5,*) ((B(K,J),J=1,8),K=1,6167)
      CLOSE (5)
      OPEN (5,FILE='MF2061.KE3',STATUS='OLD')
      READ(5,*) ((D(K,J),J=1,4),K=1,6167)
      CLOSE (5)
      OPEN (5,FILE='MF2033.KE3',STATUS='OLD')
      READ(5,*) ((E(K,J),J=1,4),K=1,5)
      CLOSE (5)
      OPEN (5,FILE='MF2061.KEK',STATUS='OLD')
      READ(5,*) ((MCASE(K,J),J=1,2),K=1,5)
      CLOSE (5)
      LCASE(1)=0
      DO 12000 L=1,5-1
      LCASE(L+1)=LCASE(L)+MCASE(L,1)
12000 CONTINUE
      DS=1
      D8=0
      D7=0
      D6=0
      D5=0
      D4=0
      D3=0
      D2=0
      D1=0
      D0=0
      DO 81500 L=1,5
      CS=1
      C8=0
      C7=0
      C6=0
      C5=0
      C4=0
      C3=0
      C2=0
      C1=0
      C0=0
      DO 80500 K=1+LCASE(L),LCASE(L)+MCASE(L,1)
      BS=1
      B8=0
      B7=0
      B6=0
      B5=0
      B4=0
      B3=0
      B2=0
      B1=0
      B0=0
      DO 80000 J=1+NCASE*(L-1),NCASE*L
      BB=0
      DO 18100 I=1,8
      BB=BB+B(K,I)*C(J,I)
18100 CONTINUE
      IF (BB .EQ. 52) THEN
      Z11=-C(J,6)*B(K,4)*A(L,1)-C(J,6)*B(K,3)*
     .A(L,2)+C(J,6)*B(K,2)*A(L,3)+C(J,6)*B(K,1)*
     .A(L,4)-C(J,5)*B(K,4)*A(L,2)+C(J,5)*B(K,
     .3)*A(L,1)+C(J,5)*B(K,2)*A(L,4)-C(J,5)*
     .B(K,1)*A(L,3)+C(J,4)*B(K,6)*A(L,1)+C(J,4)*
     .B(K,5)*A(L,2)-C(J,4)*B(K,2)*A(L,5)-C(J,
     .4)*B(K,1)*A(L,6)+C(J,3)*B(K,6)*A(L,2)-
     .C(J,3)*B(K,5)*A(L,1)-C(J,3)*B(K,2)*A(L,6)+
     .C(J,3)*B(K,1)*A(L,5)-C(J,2)*B(K,6)*A(L,
     .3)-C(J,2)*B(K,5)*A(L,4)+C(J,2)*B(K,4)*
     .A(L,5)+C(J,2)*B(K,3)*A(L,6)-C(J,1)*B(K,6)*
     .A(L,4)+C(J,1)*B(K,5)*A(L,3)+C(J,1)*B(K,
     .4)*A(L,6)-C(J,1)*B(K,3)*A(L,5)
      Z12=-C(J,6)*B(K,4)*A(L,2)+C(J,6)*B(K,3)*
     .A(L,1)+C(J,6)*B(K,2)*A(L,4)-C(J,6)*B(K,1)*
     .A(L,3)+C(J,5)*B(K,4)*A(L,1)+C(J,5)*B(K,
     .3)*A(L,2)-C(J,5)*B(K,2)*A(L,3)-C(J,5)*
     .B(K,1)*A(L,4)+C(J,4)*B(K,6)*A(L,2)-C(J,4)*
     .B(K,5)*A(L,1)-C(J,4)*B(K,2)*A(L,6)+C(J,
     .4)*B(K,1)*A(L,5)-C(J,3)*B(K,6)*A(L,1)-
     .C(J,3)*B(K,5)*A(L,2)+C(J,3)*B(K,2)*A(L,5)+
     .C(J,3)*B(K,1)*A(L,6)-C(J,2)*B(K,6)*A(L,
     .4)+C(J,2)*B(K,5)*A(L,3)+C(J,2)*B(K,4)*
     .A(L,6)-C(J,2)*B(K,3)*A(L,5)+C(J,1)*B(K,6)*
     .A(L,3)+C(J,1)*B(K,5)*A(L,4)-C(J,1)*B(K,
     .4)*A(L,5)-C(J,1)*B(K,3)*A(L,6)
      Z21=C(J,6)*B(K,4)*A(L,1)+C(J,6)*B(K,3)*
     .A(L,2)-C(J,6)*B(K,2)*A(L,3)-C(J,6)*B(K,1)*
     .A(L,4)-C(J,5)*B(K,4)*A(L,2)+C(J,5)*B(K,
     .3)*A(L,1)+C(J,5)*B(K,2)*A(L,4)-C(J,5)*
     .B(K,1)*A(L,3)-C(J,4)*B(K,6)*A(L,1)+C(J,4)*
     .B(K,5)*A(L,2)-C(J,4)*B(K,2)*A(L,5)+C(J,
     .4)*B(K,1)*A(L,6)-C(J,3)*B(K,6)*A(L,2)-
     .C(J,3)*B(K,5)*A(L,1)+C(J,3)*B(K,2)*A(L,6)+
     .C(J,3)*B(K,1)*A(L,5)+C(J,2)*B(K,6)*A(L,
     .3)-C(J,2)*B(K,5)*A(L,4)+C(J,2)*B(K,4)*
     .A(L,5)-C(J,2)*B(K,3)*A(L,6)+C(J,1)*B(K,6)*
     .A(L,4)+C(J,1)*B(K,5)*A(L,3)-C(J,1)*B(K,
     .4)*A(L,6)-C(J,1)*B(K,3)*A(L,5)
      Z22=C(J,6)*B(K,4)*A(L,2)-C(J,6)*B(K,3)*
     .A(L,1)-C(J,6)*B(K,2)*A(L,4)+C(J,6)*B(K,1)*
     .A(L,3)+C(J,5)*B(K,4)*A(L,1)+C(J,5)*B(K,
     .3)*A(L,2)-C(J,5)*B(K,2)*A(L,3)-C(J,5)*
     .B(K,1)*A(L,4)-C(J,4)*B(K,6)*A(L,2)-C(J,4)*
     .B(K,5)*A(L,1)+C(J,4)*B(K,2)*A(L,6)+C(J,
     .4)*B(K,1)*A(L,5)+C(J,3)*B(K,6)*A(L,1)-
     .C(J,3)*B(K,5)*A(L,2)+C(J,3)*B(K,2)*A(L,5)-
     .C(J,3)*B(K,1)*A(L,6)+C(J,2)*B(K,6)*A(L,
     .4)+C(J,2)*B(K,5)*A(L,3)-C(J,2)*B(K,4)*
     .A(L,6)-C(J,2)*B(K,3)*A(L,5)-C(J,1)*B(K,6)*
     .A(L,3)+C(J,1)*B(K,5)*A(L,4)-C(J,1)*B(K,
     .4)*A(L,5)+C(J,1)*B(K,3)*A(L,6)
      Z31=C(J,6)*B(K,4)*A(L,1)-C(J,6)*B(K,3)*
     .A(L,2)+C(J,6)*B(K,2)*A(L,3)-C(J,6)*B(K,1)*
     .A(L,4)+C(J,5)*B(K,4)*A(L,2)+C(J,5)*B(K,
     .3)*A(L,1)-C(J,5)*B(K,2)*A(L,4)-C(J,5)*
     .B(K,1)*A(L,3)-C(J,4)*B(K,6)*A(L,1)-C(J,4)*
     .B(K,5)*A(L,2)+C(J,4)*B(K,2)*A(L,5)+C(J,
     .4)*B(K,1)*A(L,6)+C(J,3)*B(K,6)*A(L,2)-
     .C(J,3)*B(K,5)*A(L,1)-C(J,3)*B(K,2)*A(L,6)+
     .C(J,3)*B(K,1)*A(L,5)-C(J,2)*B(K,6)*A(L,
     .3)+C(J,2)*B(K,5)*A(L,4)-C(J,2)*B(K,4)*
     .A(L,5)+C(J,2)*B(K,3)*A(L,6)+C(J,1)*B(K,6)*
     .A(L,4)+C(J,1)*B(K,5)*A(L,3)-C(J,1)*B(K,
     .4)*A(L,6)-C(J,1)*B(K,3)*A(L,5)
      Z32=C(J,6)*B(K,4)*A(L,2)+C(J,6)*B(K,3)*
     .A(L,1)-C(J,6)*B(K,2)*A(L,4)-C(J,6)*B(K,1)*
     .A(L,3)-C(J,5)*B(K,4)*A(L,1)+C(J,5)*B(K,
     .3)*A(L,2)-C(J,5)*B(K,2)*A(L,3)+C(J,5)*
     .B(K,1)*A(L,4)-C(J,4)*B(K,6)*A(L,2)+C(J,4)*
     .B(K,5)*A(L,1)+C(J,4)*B(K,2)*A(L,6)-C(J,
     .4)*B(K,1)*A(L,5)-C(J,3)*B(K,6)*A(L,1)-
     .C(J,3)*B(K,5)*A(L,2)+C(J,3)*B(K,2)*A(L,5)+
     .C(J,3)*B(K,1)*A(L,6)+C(J,2)*B(K,6)*A(L,
     .4)+C(J,2)*B(K,5)*A(L,3)-C(J,2)*B(K,4)*
     .A(L,6)-C(J,2)*B(K,3)*A(L,5)+C(J,1)*B(K,6)*
     .A(L,3)-C(J,1)*B(K,5)*A(L,4)+C(J,1)*B(K,
     .4)*A(L,5)-C(J,1)*B(K,3)*A(L,6)
      Z41=-C(J,6)*B(K,4)*A(L,1)+C(J,6)*B(K,3)*
     .A(L,2)-C(J,6)*B(K,2)*A(L,3)+C(J,6)*B(K,1)*
     .A(L,4)+C(J,5)*B(K,4)*A(L,2)+C(J,5)*B(K,
     .3)*A(L,1)-C(J,5)*B(K,2)*A(L,4)-C(J,5)*
     .B(K,1)*A(L,3)+C(J,4)*B(K,6)*A(L,1)-C(J,4)*
     .B(K,5)*A(L,2)+C(J,4)*B(K,2)*A(L,5)-C(J,
     .4)*B(K,1)*A(L,6)-C(J,3)*B(K,6)*A(L,2)-
     .C(J,3)*B(K,5)*A(L,1)+C(J,3)*B(K,2)*A(L,6)+
     .C(J,3)*B(K,1)*A(L,5)+C(J,2)*B(K,6)*A(L,
     .3)+C(J,2)*B(K,5)*A(L,4)-C(J,2)*B(K,4)*
     .A(L,5)-C(J,2)*B(K,3)*A(L,6)-C(J,1)*B(K,6)*
     .A(L,4)+C(J,1)*B(K,5)*A(L,3)+C(J,1)*B(K,
     .4)*A(L,6)-C(J,1)*B(K,3)*A(L,5)
      Z42=-C(J,6)*B(K,4)*A(L,2)-C(J,6)*B(K,3)*
     .A(L,1)+C(J,6)*B(K,2)*A(L,4)+C(J,6)*B(K,1)*
     .A(L,3)-C(J,5)*B(K,4)*A(L,1)+C(J,5)*B(K,
     .3)*A(L,2)-C(J,5)*B(K,2)*A(L,3)+C(J,5)*
     .B(K,1)*A(L,4)+C(J,4)*B(K,6)*A(L,2)+C(J,4)*
     .B(K,5)*A(L,1)-C(J,4)*B(K,2)*A(L,6)-C(J,
     .4)*B(K,1)*A(L,5)+C(J,3)*B(K,6)*A(L,1)-
     .C(J,3)*B(K,5)*A(L,2)+C(J,3)*B(K,2)*A(L,5)-
     .C(J,3)*B(K,1)*A(L,6)-C(J,2)*B(K,6)*A(L,
     .4)+C(J,2)*B(K,5)*A(L,3)+C(J,2)*B(K,4)*
     .A(L,6)-C(J,2)*B(K,3)*A(L,5)-C(J,1)*B(K,6)*
     .A(L,3)-C(J,1)*B(K,5)*A(L,4)+C(J,1)*B(K,
     .4)*A(L,5)+C(J,1)*B(K,3)*A(L,6)
      CALL SYOUKEI(Z11,Z12,BS,B8,B7,B6,B5,B4,B3,B2,B1,B0)
      CALL SYOUKEI(Z21,Z22,BS,B8,B7,B6,B5,B4,B3,B2,B1,B0)
      CALL SYOUKEI(Z31,Z32,BS,B8,B7,B6,B5,B4,B3,B2,B1,B0)
      CALL SYOUKEI(Z41,Z42,BS,B8,B7,B6,B5,B4,B3,B2,B1,B0)
      ENDIF
80000 CONTINUE
      CALL BAI(D(K,3),B8,B7,B6,B5,B4,B3,B2,B1,B0)
      WRITE(*,*) K,L
      CALL WSS(CS,C8,C7,C6,C5,C4,C3,C2,C1,C0,
     .BS,B8,B7,B6,B5,B4,B3,B2,B1,B0)
80500 CONTINUE
      CALL BAI(E(L,3),C8,C7,C6,C5,C4,C3,C2,C1,C0)
      CALL WSS(DS,D8,D7,D6,D5,D4,D3,D2,D1,D0,
     .CS,C8,C7,C6,C5,C4,C3,C2,C1,C0)
81500 CONTINUE
      WRITE(*,*) DS,D8,D7,D6,D5,D4,D3,D2,D1,D0
      OPEN(6,FILE='N207.KEK',STATUS='OLD')
      WRITE(6,90000) DS,D8,D7,D6,D5,D4,D3,D2,D1,D0
90000 FORMAT(10I4)
      CLOSE(6)
      END

\end{verbatim}
\begin{verbatim}

wa.for

      SUBROUTINE WA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      INTEGER S1,S2,S3,S4,S0,CC0
      INTEGER B0,B1,B2,B3,B4,B5,B6,B7,B8,A0,A1,A2,A3,A4,A5,A6,A7,A8
      INTEGER C0,C1,C2,C3,C4,C5,C6,C7,C8
      INTEGER C00,C10,C20,C30,C40,C50,C60,C70,C80
      INTEGER C01,C11,C21,C31,C41,C51,C61,C71,C81
      INTEGER Q0,Q1,RQ1
      INTEGER XS,X8,X7,X6,X5,X4,X3,X2,X1,X0
      INTEGER CC
      C0=A0+B0
      X0=Q0(C0)
      C01=Q1(C0)
      C1=A1+B1+C01
      X1=Q0(C1)
      C11=Q1(C1)
      C2=A2+B2+C11
      X2=Q0(C2)
      C21=Q1(C2)
      C3=A3+B3+C21
      X3=Q0(C3)
      C31=Q1(C3)
      C4=A4+B4+C31
      X4=Q0(C4)
      C41=Q1(C4)
      C5=A5+B5+C41
      X5=Q0(C5)
      C51=Q1(C5)
      C6=A6+B6+C51
      X6=Q0(C6)
      C61=Q1(C6)
      C7=A7+B7+C61
      X7=Q0(C7)
      C71=Q1(C7)
      C8=A8+B8+C71
      X8=Q0(C8)
      RETURN
      END


sa.for

      SUBROUTINE SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      INTEGER S1,S2,S3,S4,S5,S6,S7,S8,S0,CC0,CC1
      INTEGER C0,C1,C2,C3,C4,C5,C6,C7,C8
      INTEGER B0,B1,B2,B3,B4,B5,B6,B7,B8
      INTEGER A0,A1,A2,A3,A4,A5,A6,A7,A8
      INTEGER XS,X8,X7,X6,X5,X4,X3,X2,X1,X0
      INTEGER Q0,Q1,RQ1
      CC1=9999
      S0=CC1-B0
      S1=CC1-B1
      S2=CC1-B2
      S3=CC1-B3
      S4=CC1-B4
      S5=CC1-B5
      S6=CC1-B6
      S7=CC1-B7
      S8=CC1-B8
      CALL WA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .S8,S7,S6,S5,S4,S3,S2,S1,S0,
     .C8,C7,C6,C5,C4,C3,C2,C1,C0)
      CALL WA(C8,C7,C6,C5,C4,C3,C2,C1,C0,
     .0,0,0,0,0,0,0,0,1,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      END

\end{verbatim}
\begin{verbatim}

wss.for

      SUBROUTINE WSS(AS,A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .BS,B8,B7,B6,B5,B4,B3,B2,B1,B0)
      INTEGER AS,A8,A7,A6,A5,A4,A3,A2,A1,A0
      INTEGER BS,B8,B7,B6,B5,B4,B3,B2,B1,B0
      INTEGER XS,X8,X7,X6,X5,X4,X3,X2,X1,X0
      CALL WS(AS,A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .BS,B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .XS,X8,X7,X6,X5,X4,X3,X2,X1,X0)
      AS=XS
      A0=X0
      A1=X1
      A2=X2
      A3=X3
      A4=X4
      A5=X5
      A6=X6
      A7=X7
      A8=X8
      RETURN
      END

\end{verbatim}
\begin{verbatim}

ws.for

      SUBROUTINE WS(AS,A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .BS,B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .XS,X8,X7,X6,X5,X4,X3,X2,X1,X0)
      INTEGER AS,A8,A7,A6,A5,A4,A3,A2,A1,A0
      INTEGER BS,B8,B7,B6,B5,B4,B3,B2,B1,B0
      INTEGER CS,C8,C7,C6,C5,C4,C3,C2,C1,C0
      INTEGER XS,X8,X7,X6,X5,X4,X3,X2,X1,X0
      INTEGER CC
      IF (AS .NE. BS) THEN
      IF (A8 .NE. B8) THEN
      IF (A8 .GT. B8) THEN
      XS=AS
      CALL SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ELSE
      XS=BS
      CALL SA(B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      ELSE IF (A7 .NE. B7) THEN
      IF (A7 .GT. B7) THEN
      XS=AS
      CALL SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ELSE
      XS=BS
      CALL SA(B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      ELSE IF (A6 .NE. B6) THEN
      IF (A6 .GT. B6) THEN
      XS=AS
      CALL SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ELSE
      XS=BS
      CALL SA(B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      ELSE IF (A5 .NE. B5) THEN
      IF (A5 .GT. B5) THEN
      XS=AS
      CALL SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ELSE
      XS=BS
      CALL SA(B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      ELSE IF (A4 .NE. B4) THEN
      IF (A4 .GT. B4) THEN
      XS=AS
      CALL SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ELSE
      XS=BS
      CALL SA(B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      ELSE IF (A3 .NE. B3) THEN
      IF (A3 .GT. B3) THEN
      XS=AS
      CALL SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ELSE
      XS=BS
      CALL SA(B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      ELSE IF (A2 .NE. B2) THEN
      IF (A2 .GT. B2) THEN
      XS=AS
      CALL SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ELSE
      XS=BS
      CALL SA(B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      ELSE IF (A1 .NE. B1) THEN
      IF (A1 .GT. B1) THEN
      XS=AS
      CALL SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ELSE
      XS=BS
      CALL SA(B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      ELSE IF (A0 .NE. B0) THEN
      IF (A0 .GT. B0) THEN
      XS=AS
      CALL SA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ELSE
      XS=BS
      CALL SA(B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      ELSE
      XS=1
      X0=0
      X1=0
      X2=0
      X3=0
      X4=0
      X5=0
      X6=0
      X7=0
      X8=0
      RETURN
      ENDIF
      ELSE
      XS=AS
      XS=AS
      CALL WA(A8,A7,A6,A5,A4,A3,A2,A1,A0,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .X8,X7,X6,X5,X4,X3,X2,X1,X0)
      RETURN
      ENDIF
      END

\end{verbatim}
\begin{verbatim}

gg.for

      SUBROUTINE G8(RA,C70,C60,C50,C40,C30,C20,C10,C00)
      INTEGER RA
      INTEGER S1,S2,S3,S4,S0,CC0
      INTEGER B0,B1,B2,B3,B4,B5,B6,B7,B8
      INTEGER C0,C1,C2,C3,C4,C5,C6,C7,C8
      INTEGER C00,C10,C20,C30,C40,C50,C60,C70,C80
      INTEGER C01,C11,C21,C31,C41,C51,C61,C71,C81
      INTEGER Q0,Q1,RQ1,G20,G40,G41,G42,G43
      INTEGER R0
      R0=RA**2
      S0=MOD(R0,10000)
      S1=(R0-S0)/10000
      CALL G4(RA,B3,B2,B1,B0)
      C00=Q0(B0**2)
      C1=Q1(B0**2)+2*Q0(B0*B1)
      C10=Q0(C1)
      C11=Q1(C1)
      C2=C11+2*Q1(B0*B1)+2*Q0(B2*B0)+Q0(B1**2)
      C20=Q0(C2)
      C21=Q1(C2)
      C3=C21+2*Q1(B2*B0)+Q1(B1**2)+2*Q0(B0*B3)+2*Q0(B1*B2)
      C30=Q0(C3)
      C31=Q1(C3)
      C4=C31+2*Q1(B0*B3)+2*Q1(B1*B2)+2*Q0(B1*B3)+Q0(B2**2)
      C40=Q0(C4)
      C41=Q1(C4)
      C5=C41+2*Q1(B1*B3)+Q1(B2**2)+2*Q0(B2*B3)
      C50=Q0(C5)
      C51=Q1(C5)
      C6=C51+2*Q1(B2*B3)+Q0(B3**2)
      C60=Q0(C6)
      C61=Q1(C6)
      C7=C61+Q1(B3**2)
      C70=Q0(C7)
      RETURN
      END
      INTEGER FUNCTION Q0(AA0)
      INTEGER CC0,AA0
      Q0=MOD(AA0,10000)
      RETURN
      END
      INTEGER FUNCTION RQ1(AA0)
      INTEGER CC0,AA0,BB0,BB1
      BB0=MOD(AA0,10000)
      RQ1=(AA0-BB0)/10000
      RETURN
      END
      INTEGER FUNCTION Q1(AA0)
      INTEGER CC0,AA0,BB0,BB1,RQ1
      BB0=MOD(AA0,10000)
      BB1=(AA0-BB0)/10000
      Q1=MOD(BB1,10000)
      RETURN
      END
      SUBROUTINE G2(QA,C10,C00)
      INTEGER QA,C00,C10
      INTEGER Q0,Q1,RQ1
      C00=Q0(QA**2)
      C10=Q1(QA**2)
      RETURN
      END
      SUBROUTINE G4(QA,C30,C20,C10,C00)
      INTEGER QA
      INTEGER Q0,Q1,Q2,RQ1
      INTEGER S1,S2,S3,S4,S0,CC0
      INTEGER C0,C1,C2,C3,C4,C5
      INTEGER C00,C10,C20,C30,C40,C50
      INTEGER C01,C11,C21,C31,C41,C51
      INTEGER R0
      R0=QA**2
      S0=MOD(R0,10000)
      S1=(R0-S0)/10000
      C00=Q0(S0**2)
      C1=Q1(S0**2)+Q0(S0*S1)+Q0(S0*S1)
      C10=Q0(C1)
      C11=Q1(C1)
      C2=C11+Q1(S0*S1)+Q1(S0*S1)+Q0(S1**2)
      C20=Q0(C2)
      C21=Q1(C2)
      C3=C21+Q1(S1**2)
      C30=Q0(C3)
      RETURN
      END
      SUBROUTINE G6(RA,C50,C40,C30,C20,C10,C00)
      INTEGER RA
      INTEGER S1,S2,S3,S4,S0,CC0
      INTEGER B0,B1,B2,B3,B4,B5
      INTEGER C0,C1,C2,C3,C4,C5
      INTEGER C00,C10,C20,C30,C40,C50
      INTEGER C01,C11,C21,C31,C41,C51
      INTEGER Q0,Q1,RQ1
      INTEGER R0
      R0=RA**2
      S0=MOD(R0,10000)
      S1=(R0-S0)/10000
      CALL G4(RA,B3,B2,B1,B0)
      C00=Q0(S0*B0)
      C1=Q1(S0*B0)+Q0(S0*B1)+Q0(S1*B0)
      C10=Q0(C1)
      C11=Q1(C1)
      C2=C11+Q1(S0*B1)+Q1(S1*B0)+Q0(S0*B2)+Q0(S1*B1)
      C20=Q0(C2)
      C21=Q1(C2)
      C3=C21+Q1(S0*B2)+Q1(S1*B1)+Q0(S0*B3)+Q0(S1*B2)
      C30=Q0(C3)
      C31=Q1(C3)
      C4=C31+Q1(S0*B3)+Q1(S1*B2)+Q0(S1*B3)
      C40=Q0(C4)
      C41=Q1(C4)
      C5=C41+Q1(S1*B3)
      C50=Q0(C5)
      RETURN
      END
      SUBROUTINE G62(RA,RB,C70,C60,C50,C40,C30,C20,C10,C00)
      INTEGER RA,RB
      INTEGER S1,S2,S3,S4,S0,CC0
      INTEGER B0,B1,B2,B3,B4,B5
      INTEGER C0,C1,C2,C3,C4,C5,C6,C7,C8
      INTEGER C00,C10,C20,C30,C40,C50,C60,C70,C80
      INTEGER C01,C11,C21,C31,C41,C51,C61,C71,C81
      INTEGER Q0,Q1,RQ1
      INTEGER R0
      R0=RB**2
      S0=MOD(R0,10000)
      S1=(R0-S0)/10000
      CALL G6(RA,B5,B4,B3,B2,B1,B0)
      C00=Q0(S0*B0)
      C1=Q1(S0*B0)+Q0(S0*B1)+Q0(S1*B0)
      C10=Q0(C1)
      C11=Q1(C1)
      C2=C11+Q1(S0*B1)+Q1(S1*B0)+Q0(S0*B2)+Q0(S1*B1)
      C20=Q0(C2)
      C21=Q1(C2)
      C3=C21+Q1(S0*B2)+Q1(S1*B1)+Q0(S0*B3)+Q0(S1*B2)
      C30=Q0(C3)
      C31=Q1(C3)
      C4=C31+Q1(S0*B3)+Q1(S1*B2)+Q0(S0*B4)+Q0(S1*B3)
      C40=Q0(C4)
      C41=Q1(C4)
      C5=C41+Q1(S0*B4)+Q1(S1*B3)+Q0(S0*B5)+Q0(S1*B4)
      C50=Q0(C5)
      C51=Q1(C5)
      C6=C51+Q1(S0*B5)+Q1(S1*B4)+Q0(S1*B5)
      C60=Q0(C6)
      C61=Q1(C6)
      C7=C61+Q1(S1*B5)
      C70=Q0(C7)
      RETURN
      END
      SUBROUTINE G44(RA,RB,C70,C60,C50,C40,C30,C20,C10,C00)
      INTEGER RA,RB
      INTEGER S1,S2,S3,S4,S0,CC0
      INTEGER B0,B1,B2,B3,B4,B5
      INTEGER C0,C1,C2,C3,C4,C5,C6,C7,C8
      INTEGER C00,C10,C20,C30,C40,C50,C60,C70,C80
      INTEGER C01,C11,C21,C31,C41,C51,C61,C71,C81
      INTEGER Q0,Q1,RQ1
      INTEGER R0
      CALL G4(RA,S3,S2,S1,S0)
      CALL G4(RB,B3,B2,B1,B0)
      C00=Q0(S0*B0)
      C1=Q1(S0*B0)+Q0(S0*B1)+Q0(S1*B0)
      C10=Q0(C1)
      C11=Q1(C1)
      C2=C11+Q1(S0*B1)+Q1(S1*B0)+Q0(S0*B2)+Q0(S1*B1)+Q0(S2*B0)
      C20=Q0(C2)
      C21=Q1(C2)
      C3=C21+Q1(S0*B2)+Q1(S1*B1)+Q1(S2*B0)+
     .Q0(S0*B3)+Q0(S1*B2)+Q0(S2*B1)+Q0(S3*B0)
      C30=Q0(C3)
      C31=Q1(C3)
      C4=C31+Q1(S0*B3)+Q1(S1*B2)+Q1(S2*B1)+Q1(S3*B0)+
     .Q0(S1*B3)+Q0(S2*B2)+Q0(S3*B1)
      C40=Q0(C4)
      C41=Q1(C4)
      C5=C41+Q1(S1*B3)+Q1(S2*B2)+Q1(S3*B1)+Q0(S2*B3)+Q0(S3*B2)
      C50=Q0(C5)
      C51=Q1(C5)
      C6=C51+Q1(S2*B3)+Q1(S3*B2)+Q0(S3*B3)
      C60=Q0(C6)
      C61=Q1(C6)
      C7=C61+Q1(S3*B3)
      C70=Q0(C7)
      RETURN
      END
      SUBROUTINE BAI(RA,
     .B8,B7,B6,B5,B4,B3,B2,B1,B0)
      INTEGER RA
      INTEGER C0,C1,C2,C3,C4,C5,C6,C7,C8
      INTEGER B0,B1,B2,B3,B4,B5,B6,B7,B8
      INTEGER XS,X8,X7,X6,X5,X4,X3,X2,X1,X0
      INTEGER Q0,Q1,RQ1
      C0=B0*RA
      X0=Q0(C0)
      C1=Q1(C0)+Q0(B1*RA)
      X1=Q0(C1)
      C2=Q1(C1)+Q1(B1*RA)+Q0(B2*RA)
      X2=Q0(C2)
      C3=Q1(C2)+Q1(B2*RA)+Q0(B3*RA)
      X3=Q0(C3)
      C4=Q1(C3)+Q1(B3*RA)+Q0(B4*RA)
      X4=Q0(C4)
      C5=Q1(C4)+Q1(B4*RA)+Q0(B5*RA)
      X5=Q0(C5)
      C6=Q1(C5)+Q1(B5*RA)+Q0(B6*RA)
      X6=Q0(C6)
      C7=Q1(C6)+Q1(B6*RA)+Q0(B7*RA)
      X7=Q0(C7)
      C8=Q1(C7)+Q1(B7*RA)+Q0(B8*RA)
      X8=Q0(C8)
      B0=X0
      B1=X1
      B2=X2
      B3=X3
      B4=X4
      B5=X5
      B6=X6
      B7=X7
      B8=X8
      RETURN
      END

\end{verbatim}
\begin{verbatim}

syoukei.for

      SUBROUTINE SYOUKEI(RA,RB,BS,B8,B7,B6,B5,B4,B3,B2,B1,B0)
      INTEGER RA,RB
      INTEGER CS,C8,C7,C6,C5,C4,C3,C2,C1,C0
      INTEGER BS,B8,B7,B6,B5,B4,B3,B2,B1,B0
      CALL G8(RA,C7,C6,C5,C4,C3,C2,C1,C0)
      CS=1
      C8=0
      CALL WSS(BS,B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .CS,C8,C7,C6,C5,C4,C3,C2,C1,C0)
      CALL G62(RB,RA,C7,C6,C5,C4,C3,C2,C1,C0)
      C8=0
      CS=-1
      CALL BAI(28,C8,C7,C6,C5,C4,C3,C2,C1,C0)
      CALL WSS(BS,B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .CS,C8,C7,C6,C5,C4,C3,C2,C1,C0)
      CALL G44(RB,RA,C7,C6,C5,C4,C3,C2,C1,C0)
      C8=0
      CS=1
      CALL BAI(70,C8,C7,C6,C5,C4,C3,C2,C1,C0)
      CALL WSS(BS,B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .CS,C8,C7,C6,C5,C4,C3,C2,C1,C0)
      CALL G62(RA,RB,C7,C6,C5,C4,C3,C2,C1,C0)
      C8=0
      CS=-1
      CALL BAI(28,C8,C7,C6,C5,C4,C3,C2,C1,C0)
      CALL WSS(BS,B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .CS,C8,C7,C6,C5,C4,C3,C2,C1,C0)
      CALL G8(RB,C7,C6,C5,C4,C3,C2,C1,C0)
      C8=0
      CS=1
      CALL WSS(BS,B8,B7,B6,B5,B4,B3,B2,B1,B0,
     .CS,C8,C7,C6,C5,C4,C3,C2,C1,C0)
      RETURN
      END

\end{verbatim}
\begin{verbatim}

mm207.for

      PROGRAM TEST
      INTEGER J,K,I,BB,L
      INTEGER Z11,Z12,Z21,Z22,Z31,Z32,Z41,Z42
      INTEGER DS,D8,D7,D6,D5,D4,D3,D2,D1,D0
      INTEGER CS,C8,C7,C6,C5,C4,C3,C2,C1,C0
      INTEGER BS,B8,B7,B6,B5,B4,B3,B2,B1,B0
      INTEGER C(6864*48,8),B(6864*48,8),D(45377,4),E(48,4),MCASE(48,2)
      INTEGER LCASE(48),A(48,8)
      NCASE=6864
      OPEN (5,FILE='MF2033.KE2',STATUS='OLD')
      READ(5,*) ((A(K,J),J=1,8),K=1,48)
      CLOSE (5)
      OPEN (5,FILE='MF205.KE2',STATUS='OLD')
      READ(5,*) ((C(K,J),J=1,8),K=1,NCASE*48)
      CLOSE (5)
      OPEN (5,FILE='MF2061.KE2',STATUS='OLD')
      READ(5,*) ((B(K,J),J=1,8),K=1,45377)
      CLOSE (5)
      OPEN (5,FILE='MF2061.KE3',STATUS='OLD')
      READ(5,*) ((D(K,J),J=1,4),K=1,45377)
      CLOSE (5)
      OPEN (5,FILE='MF2033.KE3',STATUS='OLD')
      READ(5,*) ((E(K,J),J=1,4),K=1,48)
      CLOSE (5)
      OPEN (5,FILE='MF2061.KEK',STATUS='OLD')
      READ(5,*) ((MCASE(K,J),J=1,2),K=1,48)
      CLOSE (5)
      LCASE(1)=0
      DO 12000 L=1,48-1
      LCASE(L+1)=LCASE(L)+MCASE(L,1)
12000 CONTINUE
      DS=1
      D8=0
      D7=0
      D6=0
      D5=0
      D4=0
      D3=0
      D2=0
      D1=0
      D0=0
      DO 81500 L=1,48
      CS=1
      C8=0
      C7=0
      C6=0
      C5=0
      C4=0
      C3=0
      C2=0
      C1=0
      C0=0
      DO 80500 K=1+LCASE(L),LCASE(L)+MCASE(L,1)
      BS=1
      B8=0
      B7=0
      B6=0
      B5=0
      B4=0
      B3=0
      B2=0
      B1=0
      B0=0
      DO 80000 J=1+NCASE*(L-1),NCASE*L
      BB=0
      DO 18100 I=1,8
      BB=BB+B(K,I)*C(J,I)
18100 CONTINUE
      IF (BB .EQ. 52) THEN
      Z11=-C(J,6)*B(K,4)*A(L,1)-C(J,6)*B(K,3)*
     .A(L,2)+C(J,6)*B(K,2)*A(L,3)+C(J,6)*B(K,1)*
     .A(L,4)-C(J,5)*B(K,4)*A(L,2)+C(J,5)*B(K,
     .3)*A(L,1)+C(J,5)*B(K,2)*A(L,4)-C(J,5)*
     .B(K,1)*A(L,3)+C(J,4)*B(K,6)*A(L,1)+C(J,4)*
     .B(K,5)*A(L,2)-C(J,4)*B(K,2)*A(L,5)-C(J,
     .4)*B(K,1)*A(L,6)+C(J,3)*B(K,6)*A(L,2)-
     .C(J,3)*B(K,5)*A(L,1)-C(J,3)*B(K,2)*A(L,6)+
     .C(J,3)*B(K,1)*A(L,5)-C(J,2)*B(K,6)*A(L,
     .3)-C(J,2)*B(K,5)*A(L,4)+C(J,2)*B(K,4)*
     .A(L,5)+C(J,2)*B(K,3)*A(L,6)-C(J,1)*B(K,6)*
     .A(L,4)+C(J,1)*B(K,5)*A(L,3)+C(J,1)*B(K,
     .4)*A(L,6)-C(J,1)*B(K,3)*A(L,5)
      Z12=-C(J,6)*B(K,4)*A(L,2)+C(J,6)*B(K,3)*
     .A(L,1)+C(J,6)*B(K,2)*A(L,4)-C(J,6)*B(K,1)*
     .A(L,3)+C(J,5)*B(K,4)*A(L,1)+C(J,5)*B(K,
     .3)*A(L,2)-C(J,5)*B(K,2)*A(L,3)-C(J,5)*
     .B(K,1)*A(L,4)+C(J,4)*B(K,6)*A(L,2)-C(J,4)*
     .B(K,5)*A(L,1)-C(J,4)*B(K,2)*A(L,6)+C(J,
     .4)*B(K,1)*A(L,5)-C(J,3)*B(K,6)*A(L,1)-
     .C(J,3)*B(K,5)*A(L,2)+C(J,3)*B(K,2)*A(L,5)+
     .C(J,3)*B(K,1)*A(L,6)-C(J,2)*B(K,6)*A(L,
     .4)+C(J,2)*B(K,5)*A(L,3)+C(J,2)*B(K,4)*
     .A(L,6)-C(J,2)*B(K,3)*A(L,5)+C(J,1)*B(K,6)*
     .A(L,3)+C(J,1)*B(K,5)*A(L,4)-C(J,1)*B(K,
     .4)*A(L,5)-C(J,1)*B(K,3)*A(L,6)
      Z21=C(J,6)*B(K,4)*A(L,1)+C(J,6)*B(K,3)*
     .A(L,2)-C(J,6)*B(K,2)*A(L,3)-C(J,6)*B(K,1)*
     .A(L,4)-C(J,5)*B(K,4)*A(L,2)+C(J,5)*B(K,
     .3)*A(L,1)+C(J,5)*B(K,2)*A(L,4)-C(J,5)*
     .B(K,1)*A(L,3)-C(J,4)*B(K,6)*A(L,1)+C(J,4)*
     .B(K,5)*A(L,2)-C(J,4)*B(K,2)*A(L,5)+C(J,
     .4)*B(K,1)*A(L,6)-C(J,3)*B(K,6)*A(L,2)-
     .C(J,3)*B(K,5)*A(L,1)+C(J,3)*B(K,2)*A(L,6)+
     .C(J,3)*B(K,1)*A(L,5)+C(J,2)*B(K,6)*A(L,
     .3)-C(J,2)*B(K,5)*A(L,4)+C(J,2)*B(K,4)*
     .A(L,5)-C(J,2)*B(K,3)*A(L,6)+C(J,1)*B(K,6)*
     .A(L,4)+C(J,1)*B(K,5)*A(L,3)-C(J,1)*B(K,
     .4)*A(L,6)-C(J,1)*B(K,3)*A(L,5)
      Z22=C(J,6)*B(K,4)*A(L,2)-C(J,6)*B(K,3)*
     .A(L,1)-C(J,6)*B(K,2)*A(L,4)+C(J,6)*B(K,1)*
     .A(L,3)+C(J,5)*B(K,4)*A(L,1)+C(J,5)*B(K,
     .3)*A(L,2)-C(J,5)*B(K,2)*A(L,3)-C(J,5)*
     .B(K,1)*A(L,4)-C(J,4)*B(K,6)*A(L,2)-C(J,4)*
     .B(K,5)*A(L,1)+C(J,4)*B(K,2)*A(L,6)+C(J,
     .4)*B(K,1)*A(L,5)+C(J,3)*B(K,6)*A(L,1)-
     .C(J,3)*B(K,5)*A(L,2)+C(J,3)*B(K,2)*A(L,5)-
     .C(J,3)*B(K,1)*A(L,6)+C(J,2)*B(K,6)*A(L,
     .4)+C(J,2)*B(K,5)*A(L,3)-C(J,2)*B(K,4)*
     .A(L,6)-C(J,2)*B(K,3)*A(L,5)-C(J,1)*B(K,6)*
     .A(L,3)+C(J,1)*B(K,5)*A(L,4)-C(J,1)*B(K,
     .4)*A(L,5)+C(J,1)*B(K,3)*A(L,6)
      Z31=C(J,6)*B(K,4)*A(L,1)-C(J,6)*B(K,3)*
     .A(L,2)+C(J,6)*B(K,2)*A(L,3)-C(J,6)*B(K,1)*
     .A(L,4)+C(J,5)*B(K,4)*A(L,2)+C(J,5)*B(K,
     .3)*A(L,1)-C(J,5)*B(K,2)*A(L,4)-C(J,5)*
     .B(K,1)*A(L,3)-C(J,4)*B(K,6)*A(L,1)-C(J,4)*
     .B(K,5)*A(L,2)+C(J,4)*B(K,2)*A(L,5)+C(J,
     .4)*B(K,1)*A(L,6)+C(J,3)*B(K,6)*A(L,2)-
     .C(J,3)*B(K,5)*A(L,1)-C(J,3)*B(K,2)*A(L,6)+
     .C(J,3)*B(K,1)*A(L,5)-C(J,2)*B(K,6)*A(L,
     .3)+C(J,2)*B(K,5)*A(L,4)-C(J,2)*B(K,4)*
     .A(L,5)+C(J,2)*B(K,3)*A(L,6)+C(J,1)*B(K,6)*
     .A(L,4)+C(J,1)*B(K,5)*A(L,3)-C(J,1)*B(K,
     .4)*A(L,6)-C(J,1)*B(K,3)*A(L,5)
      Z32=C(J,6)*B(K,4)*A(L,2)+C(J,6)*B(K,3)*
     .A(L,1)-C(J,6)*B(K,2)*A(L,4)-C(J,6)*B(K,1)*
     .A(L,3)-C(J,5)*B(K,4)*A(L,1)+C(J,5)*B(K,
     .3)*A(L,2)-C(J,5)*B(K,2)*A(L,3)+C(J,5)*
     .B(K,1)*A(L,4)-C(J,4)*B(K,6)*A(L,2)+C(J,4)*
     .B(K,5)*A(L,1)+C(J,4)*B(K,2)*A(L,6)-C(J,
     .4)*B(K,1)*A(L,5)-C(J,3)*B(K,6)*A(L,1)-
     .C(J,3)*B(K,5)*A(L,2)+C(J,3)*B(K,2)*A(L,5)+
     .C(J,3)*B(K,1)*A(L,6)+C(J,2)*B(K,6)*A(L,
     .4)+C(J,2)*B(K,5)*A(L,3)-C(J,2)*B(K,4)*
     .A(L,6)-C(J,2)*B(K,3)*A(L,5)+C(J,1)*B(K,6)*
     .A(L,3)-C(J,1)*B(K,5)*A(L,4)+C(J,1)*B(K,
     .4)*A(L,5)-C(J,1)*B(K,3)*A(L,6)
      Z41=-C(J,6)*B(K,4)*A(L,1)+C(J,6)*B(K,3)*
     .A(L,2)-C(J,6)*B(K,2)*A(L,3)+C(J,6)*B(K,1)*
     .A(L,4)+C(J,5)*B(K,4)*A(L,2)+C(J,5)*B(K,
     .3)*A(L,1)-C(J,5)*B(K,2)*A(L,4)-C(J,5)*
     .B(K,1)*A(L,3)+C(J,4)*B(K,6)*A(L,1)-C(J,4)*
     .B(K,5)*A(L,2)+C(J,4)*B(K,2)*A(L,5)-C(J,
     .4)*B(K,1)*A(L,6)-C(J,3)*B(K,6)*A(L,2)-
     .C(J,3)*B(K,5)*A(L,1)+C(J,3)*B(K,2)*A(L,6)+
     .C(J,3)*B(K,1)*A(L,5)+C(J,2)*B(K,6)*A(L,
     .3)+C(J,2)*B(K,5)*A(L,4)-C(J,2)*B(K,4)*
     .A(L,5)-C(J,2)*B(K,3)*A(L,6)-C(J,1)*B(K,6)*
     .A(L,4)+C(J,1)*B(K,5)*A(L,3)+C(J,1)*B(K,
     .4)*A(L,6)-C(J,1)*B(K,3)*A(L,5)
      Z42=-C(J,6)*B(K,4)*A(L,2)-C(J,6)*B(K,3)*
     .A(L,1)+C(J,6)*B(K,2)*A(L,4)+C(J,6)*B(K,1)*
     .A(L,3)-C(J,5)*B(K,4)*A(L,1)+C(J,5)*B(K,
     .3)*A(L,2)-C(J,5)*B(K,2)*A(L,3)+C(J,5)*
     .B(K,1)*A(L,4)+C(J,4)*B(K,6)*A(L,2)+C(J,4)*
     .B(K,5)*A(L,1)-C(J,4)*B(K,2)*A(L,6)-C(J,
     .4)*B(K,1)*A(L,5)+C(J,3)*B(K,6)*A(L,1)-
     .C(J,3)*B(K,5)*A(L,2)+C(J,3)*B(K,2)*A(L,5)-
     .C(J,3)*B(K,1)*A(L,6)-C(J,2)*B(K,6)*A(L,
     .4)+C(J,2)*B(K,5)*A(L,3)+C(J,2)*B(K,4)*
     .A(L,6)-C(J,2)*B(K,3)*A(L,5)-C(J,1)*B(K,6)*
     .A(L,3)-C(J,1)*B(K,5)*A(L,4)+C(J,1)*B(K,
     .4)*A(L,5)+C(J,1)*B(K,3)*A(L,6)
     . .OR. (ABS(Z21) .GE. 1000) .OR. (ABS(Z22) .GE. 1000)
     . .OR. (ABS(Z31) .GE. 1000) .OR. (ABS(Z32) .GE. 1000)
     . .OR. (ABS(Z41) .GE. 1000) .OR. (ABS(Z42) .GE. 1000)) THEN
      WRITE(*,*) "OVERFLOW"
      OPEN(6,FILE='M207.KEK',STATUS='OLD')
      WRITE(6,*) "OVERFLOW"
      CLOSE(6)
      GOTO 99000
      ENDIF
      ENDIF
80000 CONTINUE
      WRITE(*,*) K,L
80500 CONTINUE
81500 CONTINUE
90000 FORMAT(10I4)
      WRITE(*,*) "OK"
99000 STOP
      END

\end{verbatim}
\end{tiny}

\end{document}
