24/04/30 22:29:16.36 VcpWQbIP.net
>>726
可読性向上
sim[] :=(
item=0;
L=0;
While[L<10,
p1=RandomReal[]; (* runif(1) *)
p2=RandomReal[1-p1]; (* runif(1,0,1-p1) *)
p3=1-p1-p2;
d=RandomChoice[{p1,p2,p3} -> {-1,0,1}]; (* sample(c(-1,1,1),1,prob=c(p1,p2,p3)) *)
If[!(L==0 && d==-1), L=L+d];
item++];
item
)
sim[]