If the function $f: \mathbb R \to \mathbb R$ is $n+1$ times differentiable at $x_0$ and $f^{(n+1)}(x_0) \neq 0$, then a form of the remainder in Taylor's Formula is supposedly
$$r_n(x_0;x) = \frac{f^{(n)}\big(x_0 +\theta(x-x_0)\big)}{n!}(x-x_0)^n,$$ where $ 0< \theta < 1$ and $\theta = \theta(x)$ approaches $\frac{1}{n+1}$ as $x$ approaches $x_0.$
How would one derive this?
This seems quite strange to me.
For example, taking $f(x) =e^x$, the normal Taylor expansion for $n=2$ at $0$ would have the form (Lagrange Remainder): $$e^x = 1 + x + \frac12x^2 +\frac{1}{3!}e^{\xi}x^3,$$ but in this case it is
$$e^x = 1 + x + \frac12x^2 +\frac{1}{2!}e^{\theta(x)}x^2$$.
The properties of $\theta$ do not make much sense. Am I missing something?
Using Taylor-Lagrange formula we have, for $x$ in the neighborhood of $x_0$: $$ f(x)=\sum_{k=0}^{n-1}\frac{f^{(k)}(x_0)}{k!}(x-x_0)^k+\frac{f^{(n)}(x_0+\theta_x(x-x_0))}{n!}(x-x_0)^n\tag{1} $$ with $\theta_x\in(0,1)$, and we are interested in $\lim_{x\to x_0}\theta_x$. Taking one step further in Taylor-Lagrange formula we have also $$\eqalign{ f(x)&=\sum_{k=0}^{n-1}\frac{f^{(k)}(x_0)}{k!}(x-x_0)^k\cr&\phantom{=}+\frac{f^{(n)}(x_0)}{n!}(x-x_0)^n+\frac{f^{(n+1)}(x_0+\phi_x(x-x_0))}{(n+1)!}(x-x_0)^{n+1}\tag{2}} $$ with $\phi_x\in(0,1)$. Comparing $(1)$ and $(2)$ we get $$ \frac{f^{(n)}(x_0+\theta_x(x-x_0))}{n!} =\frac{f^{(n)}(x_0)}{n!} +\frac{f^{(n+1)}(x_0+\phi_x(x-x_0))}{(n+1)!}(x-x_0) $$ or $$ \frac{f^{(n)}(x_0+\theta_x(x-x_0))-f^{(n)}(x_0)}{x-x_0}=\frac{f^{(n+1)}(x_0+\phi_x(x-x_0))}{ n+1 } $$ that is $$ \theta_x=\frac{1}{n+1}\cdot\frac{f^{(n+1)}(x_0+\phi_x(x-x_0))}{ \dfrac{f^{(n)}(x_0+h_x)-f^{(n)}(x_0)}{h_x} } $$ where $h_x=\theta_x(x-x_0)$. Letting $x$ tend to $x_0$ and noting that $f^{(n+1)}(x_0)\ne0$, we get $$\lim_{x\to x_0} \theta_x=\frac{1}{n+1}.$$ which is the desired conclusion.$\qquad\square$