Let $f_1, f_2:\mathbb{R}^2\to\mathbb{R}$ be the functions $$f_1(x,y)=\frac{3}{2}x-\frac{1}{2}y-\sin{\left(\frac{x+y}{4}\right)}\\ f_2(x,y)=(x-y)(2(x-y)^2+3(x-y)+6)+x+y$$ I want to use the gradient descent to get the solution of $\left(\begin{matrix} f_1(x,y)\\f_2(x,y) \end{matrix}\right)=0.$
I think there should be a way transforming the coordinates by a bijective function $\varphi:\mathbb{R}^2\to\mathbb{R}^2$, s.t. $\tilde f_i(z)=f_i(\varphi^{-1}(z))$ and s.t. there exists a function $F:\mathbb{R}^2\to\mathbb{R}$ where $\nabla F=\left(\begin{matrix} \tilde f_1(x,y)\\\tilde f_2(x,y) \end{matrix}\right)$.
Then we can apply the method of gradient descent to get the zeros of $\nabla F$ and hence the zeros of our system of equations. The method of gradient descent is obvious, but I cannot determine the function $\varphi$.
I am grateful for any hint.