First, let me state Cauchy's Integral formula:
Let $U$ be an open region in the complex plane and $D = \{z : |z-z_0| \leq R\}$ a disk in $U$. If $f : U \to \mathbb C$ is holomorphic and $\gamma$ is the boundary of $D$, then for all points $a$ inside the disk $$ f(a) = \frac 1{2πi}\oint_\gamma \frac{f(z)}{z-a}\ dz$$
and Cauchy's Integral Formula for the derivatives:
$$f^{(n)}(a) = \frac{n!}{2πi}\oint_\gamma \frac{f(z)}{(z-a)^{n+1}}\ dz$$
The majority of books prove the second formula by differentiation of the first formula and then proceeding by induction. However, it seems like Wikipedia wants to suggest a different proof that I couldn't find anywhere else, but they only give hints:
Since ${\displaystyle 1/(z-a)}$ can be expanded as a power series in the variable ${\displaystyle a}$: ${\displaystyle {\frac {1}{z-a}}={\frac {1+{\frac {a}{z}}+\left({\frac {a}{z}}\right)^{2}+\cdots }{z}}}$ — it follows that holomorphic functions are analytic, i.e. they can be expanded as convergent power series. In particular f is actually infinitely differentiable
So I took the time to try to understand what they are actually suggesting. What follows is my attempt to turn Wikipedia's hint into a proof, however, it's incomplete because I don't know how to justify the red-marked interchange of integral and series:
Let $a$ be a point inside $D$ and $w$ a point in the largest open disk around $a$ contained in $D$. Then, $|w-a| < |z-a|$ for all $z$ on $\gamma$, so $\left|\frac{w-a}{z-a}\right|<1$ and then $$\sum_{n=0}^\infty \frac{(w-a)^n}{(z-a)^{n+1}} = \frac{1}{z-a} \frac{1}{1-\frac{w-a}{z-a}} = \frac 1{(z-a)-(w-a)} = \frac 1{z-w} $$
Applying Cauchy's Integral formula: $$\begin{align} f(w) &= \frac 1{2πi}\oint_\gamma \frac{f(z)}{z-w}\ dz \\ &= \frac 1{2πi}\oint_\gamma \sum_{n=0}^\infty f(z)\frac{(w-a)^n}{(z-a)^{n+1}}\ dz \\ &\color{red}{=} \sum_{n=0}^\infty \frac 1{2πi}\oint_\gamma f(z)\frac{(w-a)^n}{(z-a)^{n+1}}\ dz \\ &= \sum_{n=0}^\infty \left(\frac 1{2πi}\oint_\gamma \frac{f(z)}{(z-a)^{n+1}}\ dz\right)(w-a)^n \end{align}$$ but this is a power series in $w$ around $a$. And we know that power series are infinitely-differentiable in their radius of convergence with $$ f(w) = \sum_{n=0}^\infty \frac{f^{(n)}(a)}{n!}(w-a)^n $$ therefore $$ f^{(n)}(a) = \frac{n!}{2πi}\oint_\gamma \frac{f(z)}{(z-a)^{n+1}}\ dz $$
Is the red-marked interchange valid? Apart from the interchange, is everything alright?