#!/usr/bin/env python
# coding: utf-8

# In[1]:


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()


# In[2]:


eq1 = x - a0 + b0*u + c0*v + d0*u*v
eq2 = y - a1 + b1*u + c1*v + d1*u*v


# In[3]:


sol_u, sol_v = solve([eq1, eq2], [u,v])


# In[4]:


sol_u


# In[5]:


sol_v