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


% Problem 2: 
%
fprintf( 'Problem 2: Part (b):\n' ); 
mu     = 0.3; 
lambda = 0.2; 

t = 1.0; 
q0 = mu/(mu+lambda) + (lambda/(mu+lambda))*exp(-(mu+lambda)*t);
fprintf( 'q_0(1) = %10.6f\n', q0 ); 

t = 10.0; 
q0 = mu/(mu+lambda) + (lambda/(mu+lambda))*exp(-(mu+lambda)*t);
fprintf( 'q_0(10) = %10.6f\n', q0 ); 


% Problem 3: 
%
fprintf( 'Problem 3: Part (a):\n' ); 
t_min = (1/(mu+lambda))*log( 10*lambda/mu );
fprintf( 't_min = %10.6f\n', t_min ); 

% Problem 4: 
%
fprintf( 'Problem 4:\n' ); 
T = 10*lambda/(mu*(mu+lambda));
fprintf( 'T = %10.6f\n', T );