Introducción
Análogamente al cálculo univariable, una función $f(x,y)$ es continua en un punto $(a,b)$ si se cumplen las siguientes condiciones: la función está definida en $(a,b)$, es decir $f(a,b)$ existe, el límite existe y además coincide con el valor de la función, sin embargo, en varias variables hay varias maneras de acercarse a un punto, por ejemplo: a lo largo del eje $x$, (manteniendo $y$ fijo); a lo largo del eje $y$, (manteniendo $x$ fijo), a lo largo de líneas diagonales; curvas o incluso camimnos extraños. Si el límite cambia dependiendo de la dirección de aproximación, la función no es continua en ese punto.
Proposición 1 Sea $f:\mathbb{R}^{2}\rightarrow \mathbb{R}$ tal que $$\lim_{(x,y)\rightarrow (a,b)}f(x,y)=L$$
Entonces para una función real y continua $g$ definida en un entorno da $a$ tal que $$\lim_{x\rightarrow a}g(x)=b$$ se tiene que $$\lim_{x\rightarrow a}f(x,g(x))=L$$
Demostración. Por la existencia del límite doble, dado $\epsilon>0$ existe un $\delta>0$, tal que $$|(x,y)-(a,b)|<\delta \Rightarrow |f(x,y)-L|<\epsilon.$$ Ahora $$\lim_{x\rightarrow a}g(x)=b$$ quiere decir que dado $\delta>0$ existe $\sigma>0$, con $0<\sigma<\delta$ tal que: $$|x-a|<\sigma \Rightarrow |g(x)-b|<\delta.$$ Por
tanto, si $|x-a|<\sigma$, se tiene que $|(x,g(x))-(a,b)|<\delta.$ Con lo cual, $$|f(x,g(x))-L|<\epsilon$$ $\square$
Ejemplo: Determinar si existe, el límite de la función definida por
$f(x,y)= \left\{\begin{array}{lcc} \frac{x^{2}y}{x^{2}+y^{2}}&(x,y)\neq (0,0)\\ \\0&(x,y)=(0,0)\end{array}\right.$
Para determinar su límite podemos acercarnos por trayectorias (funciones continuas) al origen.
Pongamos $y=g(x)=0$ se tiene entonces que
$$\lim_{(x,y)\rightarrow(0,0)}f(x,y)=\lim_{(x,y)\rightarrow(0,0)}f(x,g(x))=\lim_{(x,y)\rightarrow(0,0)}f(x,0)=\lim_{x\rightarrow0}\frac{x^{2}0}{x^{2}+0^{2}}=0$$
Pongamos ahora $y=g(x)=x$ se tiene entonces que
$$\lim_{(x,y)\rightarrow(0,0)}f(x,y)=\lim_{(x,y)\rightarrow(0,0)}f(x,g(x))=\lim_{(x,x)\rightarrow(0,0)}f(x,x)=\lim_{x\rightarrow0}\frac{x^{2}x}{x^{2}+x^{2}}=\lim_{x\rightarrow0}\frac{x^{3}}{2x^{2}}=0$$
Lo anterior nos dice que si existe el límite, éste tendría que ser 0, para comprobarlo usaremos la definición, se tiene entonces que debemos hallar un $\delta>0$ tal que
$\left|\frac{x^{2}y}{x^{2}+y^{2}}\right|<\epsilon$ siempre que $|(x,y)-(0,0)|<\delta$. Observamos que
$$\left|\frac{x^{2}y}{x^{2}+y^{2}}\right|=\frac{|x^{2}||y|}{|x^{2}+y^{2}|}=\frac{|x|^{2}|y|}{|x^{2}+y^{2}|}\leq
\frac{|\overline{x}|^{2}|\overline{x}|}{|\overline{x}|^{2}}=|\overline{x}|<\delta.$$
$\therefore$ podemos tomar $\delta=\epsilon$
Ejemplo: Determinar si existe, el límite de la función definida por
Para determinar su límite podemos acercarnos por trayectorias (funciones continuas) al origen.
Pongamos $y=g(x)=x$ se tiene entonces que
$$\lim_{(x,y)\rightarrow(0,0)}f(x,y)=\lim_{(x,y)\rightarrow(0,0)}f(x,g(x))=\lim_{(x,x)\rightarrow(0,0)}f(x,x)=\lim_{x\rightarrow0}\frac{x^{2}}{x^{2}+x^{2}}=\frac{1}{2}$$
Pongamos $y=g(x)=0$ se tiene entonces que
$$\lim_{(x,y)\rightarrow(0,0)}f(x,y)=\lim_{(x,y)\rightarrow(0,0)}f(x,g(x))=\lim_{(x,0)\rightarrow(0,0)}f(x,0)=\lim_{x\rightarrow0}\frac{x (0)}{x^{2}+0^{2}}=0$$
como $\frac{1}{2}\neq 0$ entonces $\cancel{\exists}$ el límite de la función.
Continuidad de Funciones de $\mathbb{R}^{n} \rightarrow \mathbb{R}$
Definición 1. Sea $f:\Omega\subset\mathbb{R}^{n} \rightarrow \mathbb{R}$, y sea $x_{0}$ un punto de acumulación de $\Omega$.Se dice que $f(x_{0})\in\mathbb{R}$ es el límite de $f$ en $x_{0}$, y se denota por: $$\displaystyle\lim_{x\rightarrow x_{0}}f(x)=f(x_{0})$$ Si dado $\varepsilon > 0$, existe $\delta > 0$ tal que $|f(x)-f(x_{0})|<\varepsilon$ cuando $x \in \Omega$, $0<|x-x_{0}|<\delta$
Ejemplo: Demostrar la continuidad en $\mathbb{R}^2$ de la función $f(x,y)=xy$.
p.d. Dado $\epsilon>0$ $\exists$ $\delta>0$ tal que $|xy-ab|\leq \epsilon$ siempre que $0<|x-a|<\delta_1$ y $0<|y-b|\leq \delta_2$ tenemos que:
$|xy-ab|
=|xy-xb+xb-ab|
\leq |x(y-b)|+|b(x-a)|
\leq \left(|x-a|+|a|\right)|y-b|+|b||x-a|\leq \left(\delta+|a|\right)\delta+|b|\delta
=\delta\underset{\underset {\text{Esta la podemos acotar}}{\searrow\ \ \ \ \ \ \ \ \ \ \ \
}}{\left(\left(\delta+|a|\right)+|b|\right)}$
Si $\delta=1$ tenemos que $\delta(1+|a|+|b|)$ y asi tomamos
$$\delta=mín \left\{ 1 , \frac{\epsilon}{1+|a|+|b|}\right\}$$
Diferenciación de funciones $\mathbb{R}^{n}\rightarrow\mathbb{R}$
Sea $f:A \subseteq\mathbb{R}^{n}\rightarrow \mathbb{R}$ y $\overline{a}=(a_{1},\ldots,a_{n})\epsilon$ $A$. Se define la derivada pacial $i$-esima en $\overline{a}$ denotada $f_{x}(\overline{a})$, $D_{x}f(\bar{a})$ ó $\displaystyle\frac{\partial f}{\partial x}(\bar{a})$ de la forma $$f_{x}=\displaystyle\lim_{h \rightarrow 0}\frac{f(a_{1},\ldots,a_{i}+h,\ldots.a_{n})-f(\bar{a})}{h}=\displaystyle\lim_{h \rightarrow 0}\frac{f(a+he_{i})-f(a)}{h}$$ siendo $\bar{e}_{i}=(0,\ldots,\underset{i-esimo}{1},\ldots,0)$.
Si $n=2$ existen 2 derivadas parciales.
Sea $\bar{a}=(x_{0},y_{0})$ un punto del interior del dominio de $f:A \subseteq\mathbb{R}^{2}\rightarrow \mathbb{R}$ las derivas parciales de $f$ en el punto $\bar{a}$ denotada respectivamente por $f_{x}(x_{0},y_{0})$, $f_{y}(x_{0},y_{0})$
son:
$$f_{x}(x_{0},y_{0})=\displaystyle\lim_{h \rightarrow 0}\frac{f(x_{0}+h,y_{0})-f(x_{0},y_{0})}{h}$$
$$f_{y}(x_{0},y_{0})=\displaystyle\lim_{k \rightarrow 0}\frac{f(x_{0},y_{0}+k)-f(x_{0},y_{0})}{k}$$
Ejemplo. Si $f(x,y)=x^{2}+x+1$ entonces
$f_{x}(0,0)=1$ ya que $f_{x}=\displaystyle\lim_{h \rightarrow 0}\frac{f(0+h,0)-f(0,0)}{h}=
\displaystyle\lim_{h \rightarrow 0}\frac{h^{2}+h+1-1}{h}=\displaystyle\lim_{h \rightarrow 0}\frac{h(h+1)}{h}=
\displaystyle\lim_{h \rightarrow 0}h+1 = 1$ y $f_{y}=\displaystyle\lim_{k \rightarrow 0}\frac{f(0,0+k)-f(0,0)}{k}=
\displaystyle\lim_{k \rightarrow 0}\frac{1-1}{k}=0$
Ejemplo. Si $f(x,y)=x^{2}+x+1$ entonces $f_{x}(0,0)=1$ ya que
$f_{x}=\displaystyle\lim_{h \rightarrow 0}\frac{f(0+h,0)-f(0,0)}{h}=
\displaystyle\lim_{h \rightarrow 0}\frac{h^{2}+h+1-1}{h}=\displaystyle\lim_{h \rightarrow 0}\frac{h(h+1)}{h}= \displaystyle\lim_{h \rightarrow 0}h+1 = 1$ y $f_{y}=\displaystyle\lim_{k \rightarrow 0}\frac{f(0,0+k)-f(0,0)}{k}= \displaystyle\lim_{k \rightarrow 0}\frac{1-1}{k}=0$
Observación: La derivada parcial en un punto de una función de varias variables en la derivada de la función de una variable, obtenida haciendo constante todas las variables, menos una. en consecuencia se pueden aplicar con esta interpretación, las reglas de derivación en una variable.
Las derivadas parciales en el punto $(x_{0},y_{0})$ de la función $z=f(x,y)$ representa la pendiente de las curvas intersección $C_{1}$ y $C_{2}$ de la superficie $z=f(x,y)$ con los planos $y=y_{0}$, $x=x_{0}$ respectivamente

Ejemplo. Calcular las derivadas parciales
$a)$ $f(x,y)=a\arcsin(x-y)$
$b)$ $f(t,u)=\displaystyle\frac{\cos(2tu)}{t^{2}+u^{2}}$
$c)$ $f(x,y,z)=\displaystyle\frac{xyz}{x^{2}+y^{2}+z^{2}}$
$d)$ $f(x,y)=\displaystyle\int_{0}^{\sqrt{xy}}e^{-t^{2}}dt \qquad x>0, \, y>0$
Solución.
$a)$ $f_{x}=\displaystyle\frac{x}{\sqrt{1-(x-y)^{2}}}+ \arcsin(x-y)$
$f_{y}=\displaystyle\frac{-x}{\sqrt{1-(x-y)^{2}}}$
$b)$ $f_{t}=\displaystyle\frac{-(t^{2}+u^{2})\sin(2tu)\cdot2u-\cos(2tu)2t}{(t^{2}+u^{2})^{2}}$
$f_{u}=\displaystyle\frac{(t^{2}+u^{2})-\sin(2tu)2u-\cos(2tu)2u}{(t^{2}+u^{2})^{2}}$
$c)$ $f_{x}=\displaystyle\frac{(x^{2}+y^{2}+z^{2})yz-xyz(2x)}{(x^{2}+y^{2}+z^{2})^{2}}$
$f_{y}=\displaystyle\frac{(x^{2}+y^{2}+z^{2})xz-xyz(2y)}{(x^{2}+y^{2}+z^{2})^{2}}$
$f_{z}=\displaystyle\frac{(x^{2}+y^{2}+z^{2})xy-xyz(2z)}{(x^{2}+y^{2}+z^{2})^{2}}$
$d)$ $f_{x}=\displaystyle e^{-xy}\frac{y}{2\sqrt{xy}}$
$f_{y}=\displaystyle e^{-xy}\frac{x}{2\sqrt{xy}}$
Más adelante
Definiremos la derivada parcial y notaremos como es similar a la derivada ordinaria a una dimensión evaluando un límite de un cociente que va incrementando en una direccion.
Tarea Moral
1.- Sea $f: \mathbb{R}^2 \rightarrow \mathbb{R}$, $(x,y) \rightarrow x^2 + y^2 + 5$, calcular el límite $\lim_{(x,y)\rightarrow(0,1)}f(x,y)$
2.- Mostrar que $f: \mathbb{R}^2 \rightarrow \mathbb{R}^2$, $(x,y) \rightarrow (x^2y,(y+x^3)/(1+x^2))$
3.- Considera la función $f(x,y)=\dfrac{sen(x^2+y^2)}{x^2+y^2}$ aunque $f$ no esté definida en $(0,0)$ determina si la función tiende a algún número cuando $(x,y)$ tiende a $(0,0)$
4.-Muestra que $f:\mathbb{R}^n \rightarrow \mathbb{R}^m$ es continua en todos los puntos si y sólo si la imagen inversa de todo abierto es abierta.
5.- Calcula las siguientes derivadas parciales :
$f(x,y)= e^{xy}log(x^2+y^2)$
$f(x,y)=\dfrac{2xy}{(x^2+y^2)^2}$