OPTION ARITHMETIC RATIONAL !有理数モード
PUBLIC NUMERIC N !上限
LET N=28
DIM F(N),A(N),B(N)
MAT F=ZER
PRINT 
PRINT "＜ １の単位分数分解 ＞"
PRINT
LET A(N)=1/N
LET B(N)=A(N)
FOR i=N-1 TO 1 STEP -1
   LET A(i)=1/i
   LET B(i)=A(i)+B(i+1)
NEXT i
CALL try(1,1,A,B,F)
PRINT
END

EXTERNAL SUB try(M,P,A(),B(),F())
OPTION ARITHMETIC RATIONAL
FOR i=P TO N 
   LET T=M-A(i)
   IF T=0 THEN
      FOR j=1 TO P-1
         IF F(j)=1 THEN
            PRINT "1/";STR$(j);" + ";
         END IF
      NEXT j
      PRINT  "1/";STR$(i)
   ELSEIF T>0 AND i<N THEN
      IF B(i+1)<T THEN EXIT FOR 
      LET F(i)=1
      CALL try(T,i+1,A,B,F)
      LET F(i)=0
   END IF
NEXT i
END SUB