@%=&1010

      PRINT "Binomial (5,3) = "; FNbinomial(5, 3)
      PRINT "Binomial (100,2) = "; FNbinomial(100, 2)
      PRINT "Binomial (33,17) = "; FNbinomial(33, 17)
      END

      DEF FNbinomial(N%, K%)
      LOCAL R%, D%
      R% = 1 : D% = N% - K%
      IF D% > K% THEN K% = D% : D% = N% - K%
      WHILE N% > K%
        R% *= N%
        N% -= 1
        WHILE D% > 1 AND (R% MOD D%) = 0
          R% /= D%
          D% -= 1
        ENDWHILE
      ENDWHILE
      = R%