I showed that the relation $$f(x,y)=e^x - e^y + xy = 0$$
defines near (0,0) an implicit function y=$\phi (x)$, since the $1x1$ block, $\frac{df}{dy}$, evaluated at (0,0) gives -1, which is non-zero - and this is also the determinant of the $1x1$ block, which makes this block invertible.
I then Taylor expanded $\phi (x)$ at the origin, to order 3, but I'm not sure whether I am computing the derivatives correctly.
The derivative of $\phi(x)$ is given by - $\frac{f_x}{f_y}$, evaluated at (0,0).
My question is: When differentiating f(x,y) w.r.t. the variable x, can I treat y as constant, as we usually do, when computing partial derivatives in multivariable calculus, or should I treat y as a function of x, so that for f(x,y) = $e^x - e^y + xy$ replace the y's with $\phi$(x), before computing the derivatives for the implied function, $\phi$(x)?
The third derivative is very messy, and if I can treat y as constant, when computing $f_x$, that would be great. I'm guessing that treating y as constant is perfectly fine, since even if I replaced the y's with $\phi$(x), then when computing $f_x$, I'd get some terms in the form of $\phi$'(x), which aren't that helpful, if I don't have an explicit function for $\phi$ - and that we only know the existence of $\phi$ by the Implicit Function Theorem.
Thanks in advance,
with $y=y(x)$ we get $e^x-e^yy'+y+xy'=0$ thus we obtain $$y'=-\frac{e^x+y}{x-e^y}$$ if $x\neq e^y$