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

lambda = 7; 
mu_1   = 8; 
mu_2   = 21; 

rho_1  = lambda/mu_1; 
rho_2  = lambda/mu_2; 

p00 = (1-rho_1)*(1-rho_2); 
fprintf( 'p_{00} = %10.6f\n', p00 );

p10 = (1-rho_1)*rho_1*(1-rho_2);
fprintf( 'p_{10} = %10.6f\n', p10 );

t = (1/2) - (1/mu_2)/(1-(lambda/mu_2)); 
fprintf('(1/2) - (1/mu_2)/(1-(lambda/mu_2)) = %10.6f\n', t); 
fprintf('mu_1 = %10.6f\n', (1+lambda*t)/t);