H'(t) = a*(H(t)/(H(t) + b) - p), H(0)=c