OPTION ARITHMETIC decimal_HIGH
PRINT
PRINT "ウィルソンの定理を用いて、１００未満の素数を見つけます。"
PRINT
DIM A(100)
LET C=0
FOR i=1 TO 99
   IF ROUND(MOD(FACT(i),i+1),10)=i THEN
      LET C=C+1
      LET A(C)=i+1
      PRINT USING "##":C;
      PRINT ": ";
      PRINT USING "##":A(C)
   END IF
NEXT i
PRINT
PRINT "１００未満の素数は、次の数たちです。"
FOR i=1 TO C
   PRINT A(i);
NEXT i
PRINT
PRINT
END