(0.5(u+v)-a)^2/b^2+(0.5(u-v)-c)^2/d^2-2p((0.5(u+v)-a)(0.5(u-v)-c))/bd