Epsilon-delta derivative proof of $x^n$ I'm currently trying to prove the power rule using the epsilon-delta definition of a derivative. I've already done it for the basic limit definition, but I thought it might be a helpful exercise to test my understanding by doing it this way. However, I'm struggling and would really appreciate any help/hints.
My work so far:
The epsilon-delta definition says a function is differentiable if, for every $\epsilon > 0$, there exists a $\delta > 0$ such that $|x-x_0| < \delta$ implies $\frac{f(x)-f(x_0)}{x-x_0} - f'(x_0) < \epsilon$. Thus, I need to find a delta, as function of epsilon, and possibly $x_0$, such that whenever the first inequality is true the second one is as well.
So, fix $x_o \in R$, and let $|x-x_0| < \delta$, then $|\frac{x^n - x_{0}^{n}}{x-x_0} - nx_{0}^{n-1}| = |\frac{x^n - nxx_{0}^{n-1} + (n-1)x_{0}^{n}}{x-x_0}|$.
Upper-bounding this with $x < \delta + x_0$ and simplifying, we get $|\frac{(\delta + x_0)^n - n\delta x_{0}^{n-1}-x_{0}^{n}}{\delta}| = |\frac{(\delta + x_0)^n}{\delta} -nx_{0}^{n-1} - x_{0}^{n}| < \epsilon$. 
Is that right so far, and if so, any advice on how to push through this last step and solve for $\delta$? I'm stumped, to be honest.
-Thank you in advance for your help!

For $x\ne 0$ we have $$(x+d)^n-x^n=-x^n+\sum_{j=0}^n x^{n-j}d^j\binom {n}{j}=nx^{n-1}d +d^2 P(d)$$ where $P(d)$ is a polynomial.   So there exists $M\in \Bbb R^+$ such that $\forall d\; (|d|\leq 1\implies |P(d)|<M).$  Therefore for $0<|d|\leq 1$ we have $$|\frac {(x+d)^n-x^n}{d}-nx^{n-1}| =$$ $$= |\frac {nx^{n-1}d+d^2P(d)}{d}-nx^{n=1}| =$$ $$= |dP(x)|\leq |d|M.$$ And $|d|M\to 0$ as $d\to 0.$
To prove that $M$ exists : If $P(d)=\sum_{j=0}^ma_jd^j$ then $|d|\leq 1\implies |P(d)|\leq$ $\sum_{j=0}^m|a_j|\cdot |d|^j\leq $ $\leq \sum_{j=0}^m|a_j|=M.$ 
For $x=0\ne d$ we have $\frac {(x+d)^n-x^n}{d}=\frac {(0+d)^n-0^n}{d}=d^{n-1}.$
