OPTION ARITHMETIC DECIMAL_HIGH
OPTION BASE 0
PRINT 
LET x=SQR(31)
PRINT
PRINT "SQR(31) の連分数表示： ";
DIM R(20)
LET y=1
LET s=0
FOR j=9 TO 1 STEP -1
   LET R(j)=INT(x/y)
   PRINT R(j);
   LET s=s+1
   LET z=MOD(x,y)
   LET x=y
   LET y=z
NEXT j
PRINT " ・ ・ ・ ・ ・"
PRINT

FOR t=1 TO 9
   LET a=R(t)
   LET b=R(t+1)
   LET c=R(t+2)
   LET d=R(t+3)
   LET e=R(t+4)
   LET f=R(t+5)
   LET g=R(t+6)
   LET h=R(t+7)
   LET i=R(t+8)
   LET K=a*b+1
   LET L=K*c+a
   LET M=L*d+K
   LET N=M*e+L
   LET O=N*f+M
   LET P=O*g+N
   LET Q=P*h+O
   IF (Q*i+P)^2-31*Q^2=1 THEN 
      LET V=(Q*i+P)
      LET W=Q
   ELSEIF (P*h+O)^2-31*P^2=1 THEN
      LET V=(P*h+O)
      LET W=P
   ELSEIF (O*g+N)^2-31*O^2=1 THEN
      LET V=(O*g+N)
      LET W=O
   ELSEIF (N*f+M)^2-31*N^2=1 THEN
      LET V=(N*f+M)
      LET W=N
   ELSEIF (M*e+L)^2-31*M^2=1 THEN
      LET V=(M*e+L)
      LET W=M
   ELSEIF (L*d+K)^2-31*L^2=1 THEN
      LET V=(L*d+K)
      LET W=L
   ELSEIF (K*c+a)^2-31*K^2=1 THEN
      LET V=(K*c+a)
      LET W=K
   ELSEIF (a*b+1)^2-31*a^2=1 THEN
      LET V=(a*b+1)
      LET W=a
   END IF
NEXT t
PRINT "χ ＝";V;"，  ｙ ＝";W
PRINT
PRINT "SRQ(31) ＝";ROUND(SQR(31),7);" ≒ ";ROUND(V/W,7);"＝ χ/y ＝";V;"/";W
PRINT
END