% 
% 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.
% 
%-----

clc; 
close all; drawnow; 

mu_min = (1/12)*( 1 + sqrt( 100/(4800/12) ) ); 
fprintf( 'mu min = %10.6f\n', mu_min ); 

lambda = 1/12;
rho    = lambda/mu_min; 
fprintf( 'p_0 = %10.6f\n', 1-rho); 

t = 24; 
p = rho*exp(-(mu_min-lambda)*t); 
fprintf( 'p = %10.6f\n', p);