% 
% 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 = 3; 
mu     = 4;
rho    = lambda/mu; 

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

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

PT = gamcdf(1/2,2,1/mu);
fprintf( '1 - F_T(t=1/2) = %10.6f\n', 1-PT );