24/04/14 15:07:11.15 T4z17oY+.net
>>163
仕様書をみながらWolfram言語に移植してみた。推敲歓迎。
sim := (
n0=10000;
n=n0;
n20=20;
n10=10;
count=1;
atari=Total[Boole[Table[RandomInteger[n,n10][[i]] <= n20,{i,n10}]]];
While[atari==0,n-=20;count++;atari=Total[Boole[Table[RandomInteger[n,n10][[i]] <= n20,{i,n10}]]]];
li={count,atari}
)
k=10^5
re=Table[sim,k]
p=2500
a=1000000
balance=Table[-p*n10*re[[i,1]]+re[[i,2]]*a,{i,k}]
N[Mean[balance]]
N[Median[balance]]