[djf]

I have regular plain old BASIC somewhere. Suppose if I got ambitious, I could run it through ten thousand trials or so.

Then convert the program to do roulette...

;-)

[Lonesome in Massachussets]

% some lunch time matlab, at your service

>> heads = rand(100, 1e5)> .5; % 100,000 trials of coin 100 tosses

>> total_heads = sum(heads); % number of heads

>> mean(total_heads)

ans =

49.9962 % we expect about 50;

>> std(total_heads)

ans =

4.9998 % we expect about 5

>> lx =length( find(mean_x >=45 & mean_x <= 55));

%% “expected” number between 45 and 55
>> fprintf( ‘%20.0f\n’, sum( binopdf( 45:55, 100, .5))*1e5)

72875

%% experimental number between 45 and 55

>> fprintf( ‘%20.0f\n’, lx)

72523