27 lines
497 B
Text
27 lines
497 B
Text
PROGRAM BINOMIAL
|
|
|
|
!$DOUBLE
|
|
|
|
PROCEDURE BINOMIAL(N,K->BIN)
|
|
LOCAL R,D
|
|
R=1 D=N-K
|
|
IF D>K THEN K=D D=N-K END IF
|
|
WHILE N>K DO
|
|
R*=N
|
|
N-=1
|
|
WHILE D>1 AND (R-D*INT(R/D))=0 DO
|
|
R/=D
|
|
D-=1
|
|
END WHILE
|
|
END WHILE
|
|
BIN=R
|
|
END PROCEDURE
|
|
|
|
BEGIN
|
|
BINOMIAL(5,3->BIN)
|
|
PRINT("Binomial (5,3) = ";BIN)
|
|
BINOMIAL(100,2->BIN)
|
|
PRINT("Binomial (100,2) = ";BIN)
|
|
BINOMIAL(33,17->BIN)
|
|
PRINT("Binomial (33,17) = ";BIN)
|
|
END PROGRAM
|