esim.m

来自「一些常被用于教学或者参考的概率论的实例的源代码」· M 代码 · 共 20 行

%esim.m/created by PJNahin for "Duelling Idiots"(12/20/98)
%This m-file estimates e by performing the simulation described
%in the statement of Probblem 5.
%
%
rand('state',100*sum(clock));       %new seed for generator;
N=input('What is N? ')
bin=zeros(1,N);
for j=1:N
   number=floor(N*rand)+1;          %floor rounds towards minus infinity;
   bin(number)=bin(number)+1;
end
empty=0;                            %empty = # of empty bins;
for j=1:N
   if bin(j)==0
   empty=empty+1;
   end
end
e_estimate=N/empty