Para 2 variables creo que la siguiente expresión sea correcta. Si es incorrecto, me gustaría saber!
$$f(x,y) = \sum_{m=0}^{\infty} \sum_{n=0}^{\infty} \beta_{m,n} (x-x_0)^m (y-y_0)^n$$
$$\beta_{m,n} = \frac{1}{n! m!} \left[ \frac{\partial^{m+n}}{\partial x^m \partial y^n} \ f(x,y) \right]_{(x=x_0,y=y_0)}$$
Donde la notación de corchetes con coordenadas subíndice significa evaluar.
Yo antes derivados de la aplicación de la siguiente fórmula, sin embargo ahora me sospechar que no es la correcta. Quizás alguien mejor que yo en la prueba por inducción podría comprobar esto?
$$f(x,y) = \sum_{k=0}^{\infty} \frac{1}{k!} \sum_{j=0}^{k} \alpha_{k,j} (x-x_0)^{j} (y-y_0)^{k-j}$$
$$\alpha_{k,j} = \,^{k}C_{j} \left[ \frac{\partial^{k}}{\partial x^j \partial y^{k-j}} \ f(x,y) \right]_{(x=x_0,y=y_0)}$$
Con el fin de demostrar la primera forma, hice lo siguiente:
$$\left[ f(x,y) \right]_{(x0,y0)} = \sum_{m=0}^{\infty} \sum_{n=0}^{\infty} \beta_{m,n} \cdot 0^m \cdot 0^n = \beta_{0,0}$$
La diferenciación parcial con respecto a x nos encontramos con
$$\left[ \frac{\partial f(x,y)}{\partial x} \right]_{(x0,y0)} = \sum_{m=0}^{\infty} \sum_{n=0}^{\infty} \left( \beta_{m,n} \cdot m \cdot 0^{m-1} \cdot 0^n + \left(\frac{\partial \beta_{m,n}}{\partial x} \right) \cdot0^m\cdot0^n\right) = \beta_{1,0} $$ where $\tfrac{\parcial \beta_{m,n}}{\partial x} = 0$.
De continuar con este proceso para todas las combinaciones de las derivadas parciales con respecto a x e y, nos encontramos con que, en general, $$\beta_{m,n}=\frac{1}{m!n!}\cdot \left[ \frac{\partial f(x,y)}{\partial x} \right]_{(x=x0,y=y0)}$$
A menos que yo haya cometido un error que es. (En cuyo caso, por favor, añadir un comentario para que yo pueda corregir mi respuesta.)