function v = prob4f(t,y,Z,b,c)
v = b * Z * (1 - y) - c * y;