In high school, we learned to reason like so:
$$(*) \qquad \frac{d}{dx}(x^2+x) = \frac{d}{dx}(x^2)+\frac{d}{dx}(x) = 2x+1$$
Now that I know more, I can "reanalyze" this chain of reasoning using ideas that I have more faith in, like "function" and "lambda abstraction." We begin by defining $\nabla (f)$ as the derivative of $f$. Then the above chain of reasoning becomes:
$$\nabla\mathop{\lambda}_{x:\mathbb{R}}(x^2+x) = \nabla\mathop{\lambda}_{x:\mathbb{R}}(x^2)+\nabla\mathop{\lambda}_{x:\mathbb{R}}(x) = \left(\mathop{\lambda}_{x:\mathbb{R}}2x\right) +\left(\mathop{\lambda}_{x:\mathbb{R}}1\right) = \mathop{\lambda}_{x:\mathbb{R}}(2x+1)$$
So I can confidently say that I understand $(*)$, because I can reanalyze it in terms of ideas that I have a lot of faith in, like "function" and "lambda abstraction."
Onwards.
In high school, we also learned a pattern of reasoning that was referred to as "implicit differentiation." It looks a bit like so:
Suppose $y^2+x = x^2+y.$
Then: $$\frac{d}{dx}(y^2+x) = \frac{d}{dx}(x^2+y).$$
$$\therefore 2y \frac{dy}{dx}+1 = 2x+ \frac{dy}{dx}$$
$$\therefore (2y-1)\frac{dy}{dx} = 2x-1$$
$$\therefore \left(\frac{dy}{dx} = \frac{2x-1}{2y-1}\right) \vee (y = 1/2)$$
Unfortunately, I still have absolute no idea what any of this means.
Question. How can we reanalyze implicit differentiation using respectable concepts like "function" and "lambda abstraction" and "limit", and without using concepts such as as "dependent variable" and "independent variable" and "differential."
Edit. It seems to be unclear what I'm looking for. I'm not looking for handwaiving and intuition. I want something as formal as possible, with the minimum of handwaiving. For example, if you're going to "switch" semantics so that $y^2+x=x+y^2$ is no longer a condition on pairs of points $(x,y) \in \mathbb{R}^2$ but instead becomes a condition on smooth functions $(a,b) \rightarrow \mathbb{R}^2,$ you should make that completely clear, introduce notation for the set of paths, etc., and the rest of your answer should be phrased in terms of this notation. Write definitions. State theorems if relevant. Tell me the domains and codomains of all your functions. I want the absolute technical logical and/or set-theoretic nitty-gritty here.
The context of this is that we're looking at the subset of the plane $$ \{(x,y)\in\mathbb R^2 \mid y^2+x=x^2+y \}$$ which is a curve of a certain shape, and we're trying to find a parameterization of that curve -- or at least part of it -- that is, nontrivial functions smooth functions $x(t)$ and $y(t)$ such that $y(t)^2+x(t)=x(t)^2+y(t)$ for all $t$ in an interval.
In order to simplify things before we start, we decide to choose the parameterization such that $t$ is just the $x$-coordinate -- in other words, $x$ is now the identity function, so we can forget about the $t$ and $y$ is now secretly a function of $x$ (whereas $x$ itself is again just a number)!
With lambdas, the original equation is then $$ \tag{*} (\lambda x.(y(x))^2+x) = (\lambda x.x^2+y(x)) $$ where the unknown is now the function $y:\mathbb [a,b]\to\mathbb R$ on some yet-to-be-identified interval $[a,b]$.
Implicit differentiation now applies your $\nabla$ to both sides of $\text{(*)}$, and after some routine manipulations we find that $\text{(*)}$ implies is $$ \forall x\in[a,b]. \left((\nabla y)(x) = \frac{2\cdot x-1}{2\cdot y(x)-1} ~\lor~ y(x)=\frac12\right) $$ We can then use that to find $y'(x)$ -- and thus the tangent to the curve -- at a particular point that we already know is on it, without needing to derive an explicit expression for $y$. (Doing so would be an easy application of the quadratic formula in this particular case, but the technique also applies in cases where this is not easy).