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

lambda = 2; 
mu     = 3;
rho    = lambda/mu; 
T      = 100/24;     % <- time in days 

fprintf( 'T*lambda = %10.6f\n', T*lambda );

W_Q = rho/((1-rho)*mu);
fprintf( 'W_Q = %10.6f (hours) = %10.6f (minutes)\n', W_Q, W_Q*60 ); 

L_Q = lambda*W_Q; 
fprintf( 'L_Q = %10.6f\n', L_Q ); 

p0 = 1-rho; 
fprintf( 'p0 = %10.6f; 73/226 = %10.6f\n', p0, 73/226 );

fprintf( '1/mu = %10.6f (hours) = %10.6f (minutes)\n', 1/mu, (1/mu)*60 );

p1 = rho; 
fprintf( 'p1 = %10.6f\n', p1 );