function dxdt=webservcon(t,x, flags, r)
mu=1;
t1=20;t2=25;
if t<t1
   lambda=0.5;
elseif t<t2
   lambda=4;
else
   lambda=0.5;
end

u=1;
dxdt=lambda*u-mu*x/(1+x);