PRINT
PRINT "成分も大きさも自然数で表される３次元ベクトル。"
PRINT
LET n=0
LET m=0
FOR i=1 TO 100
   FOR j=i+1 TO 100
      FOR k=j+1 TO 100
         LET R = i^2 + j^2 + k^2
         IF R > 10000 THEN GOTO 10
         IF SQR(R) = INT(SQR(R)) THEN
            LET n=n+1
            LET m=m+1
            LET X=gcd(i,j)
            LET Y=gcd(X,k)
            IF Y>1 THEN
               LET m=m-1
               PRINT "* ";
            END IF
            PRINT "(";i;",";j;",";k;") → ";SQR(R);
            PRINT " ： ";i^2;"+";j^2;"+";k^2;" → ";R;" ← ";SQR(R);"^2"
         END IF
10     NEXT k
   NEXT j
NEXT i
PRINT
PRINT "全部で";n;"個ありました。"
PRINT "基本的なカルテットは全部で";m;"個です。"
PRINT
PRINT "　以上が、 成分も大きさもすべて自然数の３次元ベクトルです。"
PRINT "ただし、大きさが１００以下のものです。"
PRINT
END

EXTERNAL FUNCTION gcd(a,b)     !最大公約数
   DO UNTIL b=0
      LET t=b
      LET b=MOD(a,b)
      LET a=t
   LOOP
   LET gcd=a
END FUNCTION