function yp = drugmetab(t,y)

a21 = log(2)/0.5;
a02 = log(2)/5;
subt = t-floor(t/6)*6;
D = zeros(size(t));
D(subt<0.5) = 4;
yp = [-a21*y(1)+D; a21*y(1)-a02*y(2)];