from __future__ import division
from sympy import *
x, y, u, v = symbols('x y u v')
a0, b0, c0, d0 = symbols('a0 b0 c0 d0')
a1, b1, c1, d1 = symbols('a1 b1 c1 d1')
init_printing()
eq1 = x - a0 + b0*u + c0*v + d0*u*v
eq2 = y - a1 + b1*u + c1*v + d1*u*v
sol_u, sol_v = solve([eq1, eq2], [u,v])
sol_u
sol_v