function v = prob5f(t,y,Z,lambda)
global a b d d1 e gamma T
I = y(1);
N = y(2);
Nlag = Z(2,1);
dIdt = lambda*(N - I)*(I/N) - ( d + e + gamma)*I;
dNdt = b*exp(-a*Nlag)*Nlag*exp(-d1*T) - d*N - e*I;
v = [dIdt; dNdt];