function v = prob6f(t,y,Z)
global a b d d1 T
v = b * exp(-a*Z)*Z*exp(-d1*T) - d*y;