Where λ controls the weight given to the smoothness
constraint and Ω denotes the image domain. Also we assume
that u,v are zero at the boundaries of the image domain. The
minimization of (5) is achieved using Calculus of Variations
(CV) and the approximation of the Laplacian