Why does one need to use $f(1,0)$ a linear approximation, rather than calculating $f(1.1, -0.1)$ directly?
e.g. when $f(x,y)=xe^{xy}$ (used to be more complicated $e^{xy}+xye^{xy}$)
It's easy to see that $f$ is a not a linear function. however I cannot recall, why is it uncomputable on fractional points?
Example given here:
Note: I did this using your original function: $f(x,y)=e^{xy}+xye^{xy}$, the point still stands.
This result is quite computable, but it may be impractical to compute. It may be more worth while to simply get a quite very strong approximation. I am by no means experienced with these sorts of things, but when it comes to some types of physics or engineering, sometimes a linear approximation is good enough for the purpose. Worrying about the specific value can be just a waste of time and computing power in practice.
In this case, we have by approximation that $$f(x,y)\approx f(x_0,y_0)+f_x(x_0,y_0)(x-x_0)+f_y(x_0,y_0)(y-y_0). $$ So, if we approximate at $(1,0)$, we have $$ e^{xy}+xye^{xy}\approx 1+(2ye^{xy}+xy^2e^{xy})\bigg|_{(1,0)}(x-1)+ (2xe^{xy}+x^2ye^{xy})\bigg|_{(1,0)}(y-0)=1+2y.$$ So, we have that $f(1.1,-0.1)\approx 1+2(-0.1)=.8$. If we use Wolfram Alpha to compute the result directly, we find that $f(1.1,-0.1)=.79729...$ so that our error is extremely small.