PRINT
PRINT "    c=A∩B         d=NOT(A∪B)"
PRINT "    a=A∩NOT(B)    b=NOT(A)∩B"
PRINT "    K=a+c+b+d"
PRINT "    4*(a+c)=K      3*(b+c)=K"
PRINT
PRINT "    a    c    b    d    K  カイ二乗 （ 3.84 以上で関係性あり ）"
LET X=O
LET Y=0
LET Z=0
FOR c=1 TO 30
   FOR a=1 TO 30
      FOR b=a TO 30
         LET d=-13*a+11*b-c
         LET K=12*(b-a)
         IF a+c=K/4 AND b+c=K/3 THEN
            CALL KOUY(a,c,b,d,W)
            LET R=KAI(c,a,b,d)
            IF R >= 3.34 THEN
               LET Z=Z+1
            END IF
            PRINT USING "#####":a;
            PRINT USING "#####":c;
            PRINT USING "#####":b;
            PRINT USING "#####":d;
            PRINT USING "#####":K;
            PRINT "    ";R;
            LET X=X+1
            IF 2*d=K AND 12*c=K THEN
               PRINT "    ←  完全独立事象のケース";
               LET Y=Y+1
            END IF
            IF W=1 THEN
               PRINT "  （ 公倍数あり ）";
            END IF
            PRINT
         END IF
      NEXT b
   NEXT a
NEXT c
PRINT
PRINT "以上 a,b,c 30 以下の";X;"ケースを見たうち、"
PRINT "     完全独立事象のケースは";Y;"ケースでした。"
PRINT "     明らかに関係性のあるケースは";Z;"ケースでした。"
PRINT "     c>15 では、すべて明らかに関係性のあるケースでした。"
PRINT
END

EXTERNAL SUB KOUY(q1,q2,q3,q4,Q)
LET Q=0
FOR i=2 TO q2
   IF (INT(q1/i)-q1/i)+(INT(q2/i)-q2/i)+(INT(q3/i)-q3/i)+(INT(q4/i)-q4/i)=0 THEN
      LET Q=1
      EXIT FOR
   END IF
NEXT i
END SUB

EXTERNAL FUNCTION KAI(q1,q2,q3,q4)
LET U=(q1+q2+q3+q4)*(q1*q4-q2*q3)^2/((q1+q3)*(q2+q4)*(q1+q2)*(q3+q4))
LET KAI=ROUND(U,2)
END FUNCTION