PROBLEM FORMULATION:
Considering the region \Omega bounded as a square box within x \in [0,1], y \in [0,1] . We wish to solve the 2D, stationary, advection-diffusion equation,
0 = D\nabla^2 \rho(x,y) + \vec{V} \cdot \nabla \rho(x,y)
where D is a scalar constant, and \vec{V} =...