% queue_model.m - queing modeling
% KJA, 2006

function dxdt=queue_model(t, x)

% Initialize parameters
mu=1;
t1=20;t2=25.1;

% Determine the loading condition
if t<t1
   lambda=0.5;
elseif t<t2
   lambda=4;
else
   lambda=0.5;
end

% Compute the change in queue size
dxdt=lambda-mu*x/(1+x);
