% 
% Written by:
% -- 
% John L. Weatherwax                2007-09-10
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

alpha = 10; 

s = alpha; 
b = 1/sqrt(2*pi*10); 
fprintf( 'b = %10.6f\n', b );

% we need to increase s:
while( b > 0.05 )
  s = s+1;
  b = (1/sqrt(2*pi*s))*((exp(1)*alpha/s)^s)*exp(-alpha); 
  fprintf( 's = %10d; b = %10.6f\n', s, b );
end