d/dp [4p(1-p) + p^2 + q(4(1-p) + p)]/(q+1))