ベル数

[0001 some one 2023/08/03(木) 17:41:58.16 ID:r+lEYW0A]

REM 最初の100個のベル数を出力します。十進basicを使っています。1000桁モードで実行してください。
OPTION BASE 0
DIM B(105)
LET B(0)=1
LET B(1)=1

LET s=0

FOR n=1 TO 100

FOR k=0 TO n
LET s=s+comb(n,k)*B(k)
NEXT K
PRINT n;s

LET B(n+1)=s
LET s=0
NEXT N

END