% 
% 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 = 1/10;
mu     = 1/8;
M      = 4; 
s      = 2; 
p_n = machine_repair_model(lambda,mu,M,s); 

fprintf('the average number of machines not in use is %10.6f\n', dot( 0:M, p_n ) ); 
fprintf( 'the coefficient of loss for the machines is given by %10.6f\n', C_machines(p_n,M,s) );
fprintf( 'the proportion of time that both repairmen are busy is %10.6f\n', sum( p_n(3:end) ) ); 
fprintf( 'the coefficient of loss for the repairmen is %10.6f\n', C_repairmen(p_n,s) );