バグ修正して動作確認
６面サイコロを１０回投げた場合

> # N面体のサイコロをn回投げてm種類の目がでている確率を計算
> Dice=\(N,n,m,fraction=FALSE){
+ j=m:1
+ k=(-1)^(0:(m-1))
+ nu=ifelse(m==1,N,choose(N,m)*sum(k*choose(m,j)*(j^n-j)))
+ de=N^n
+ gcd=numbers::GCD(nu,de)
+ if(fraction) cat(nu,'/',de,'=',nu/gcd,'/',de/gcd,'=',nu/de,'\n')
+ return(nu/de)
+ }
> Dice=Vectorize(Dice)
> Dice(6,10,1:6,T)
6 / 60466176 = 1 / 10077696 = 9.922903e-08
15330 / 60466176 = 2555 / 10077696 = 0.0002535302
1119600 / 60466176 = 7775 / 419904 = 0.01851614
12277800 / 60466176 = 170525 / 839808 = 0.2030524
30618000 / 60466176 = 875 / 1728 = 0.5063657
16435440 / 60466176 = 38045 / 139968 = 0.2718121