2*(c-a+2*d*t)*d-(d-b-c*t)*c where t = (-a^2 + a c - 2 b^2 + 2 b d)/(2 (b c - a d))