Mean value theorem functional equation: $ f ' \left ( \frac { x + y } 2 \right ) = \frac { f ( x ) - f ( y ) } { x - y } $ I need help solving the following functional equation. I found this question in the book "Elementary Real Analysis", right after the section on the mean value theorem.

Find all differentiable $ f : \mathbb R \to \mathbb R $ such that
$$ f ' \left ( \frac { x + y } 2 \right ) = \frac { f ( x ) - f ( y ) } { x - y } $$
holds for all $ x \ne y $.

I can clearly see that any function of the form $ f ( x ) = \alpha x $ solves this, but I don't know how to find any other solutions or show that this is the only one.
 A: The functional equation
$$
f'\left(\frac{x + y}{2}\right) = \frac{f(y) - f(x)}{y - x}\quad\text{for all real $x$ and $y$,}
$$
guarantees $f'$ is continuous: Fixing $x$ and letting $y = x + 2h \to x$, we have
$$
\lim_{h \to 0} f'(x + h) = \lim_{h \to 0} \frac{f(x + 2h) - f(x)}{2h} = f'(x).
$$
This technical fact is needed below.
For each real $z$ and for all $h \neq 0$, the functional equation says
$$
2h\, f'(z) = f(z + h) - f(z - h).
$$
Fixing $z$ and differentiating with respect to $h$ (each side is differentiable in $h$ by hypothesis),
$$
2f'(z) = f'(z + h) + f'(z - h)\quad\text{for all real $z$ and $h$.}
$$
In other words, $2f'(\frac{x + y}{2}) = f'(x) + f'(y)$ for all real $x$ and $y$. Since $f'$ is continuous, it follows that $f'$ is a linear polynomial, i.e., there exist real numbers $a$ and $b$ such that
$$
f'(x) = 2ax + b\quad\text{for all real $x$,}
$$
and consequently that $f$ is a quadratic polynomial.
This fact is well-known (and surely answered elsewhere on site), but in the spirit of being self-contained, here's an ad hoc sketch: If $f'(x)$ and $f'(y)$ are known for some $x < y$, then


*

*$f'(z) = f'(x) + \dfrac{z - x}{y - x} \bigl(f'(y) - f'(x)\bigr)$ on $[x, y]$ by successive bisection and continuity;

*$f'(2y - x) = f'(y + (y - x)) = 2f'(y) - f'(x)$ and $f'(y - 2x) = f'(y) - 2f'(x)$ are uniquely determined, so $f'$ is uniquely determined (hence linear) on the larger interval $[y - 2x, 2y - x]$. Inductively, $f'$ is linear on (arbitrary closed, bounded subintervals of) $(-\infty, \infty)$.

(Once the end result is known, more elegant arguments can doubtless be given, based on the fact that the functional equation is invariant if a quadratic polynomial is subtracted from $f$; interpolate $f$ at three points with a quadratic $ax^{2} + bx + c$, put $g(x) = f(x) - (ax^{2} + bx + c)$, and show $g \equiv 0$. I haven't carried out the details, however.)
A: It's easy to show that any function $f(x)=ax^2+bx+c$ fulfills the equation.  I remember that this result was discovered by Archimede.
A: You can maybe try restricting yourself to certain families. For example the following drafts an informal proof that, of the family of functions that can be expressed as Maclaurin series,  only  the polynomials up to 2nd order have this property. This family is quite big as it includes all polynomials, exp, sin, cos, log(1+x), ... 
Write $f(x) = \sum_{i=0}^{\infty} \alpha_i x^i $. Replacing in your equation you should obtain
$$\sum_{i=1}^{\infty} \frac {i \alpha_i}{2^{i-1}} (x+y)^{i-1} = 
\frac{ \sum_{i=0}^{\infty} \alpha_i x^i -  \sum_{i=0}^{\infty} \alpha_i y^i}{x-y}.$$
Then multiply both sides by $(x-y) $:
$$\sum_{i=1}^{\infty} \frac {i \alpha_i}{2^{i-1}} (x+y)^{i-1} (x-y) = \sum_{i=0}^{\infty} \alpha_i (x^i - y^i).$$
I'll now give the gist of the proof. We need to find values of $\alpha_i$ for which the above identity holds for all $x$ and $y$. This means that the coefficients of all monomials $x^jy^k$ must be the same on the lhs and hrs. We notice that on both sides $\alpha_i$ multiplies all and only monomials $x^jy^k$ such that $i=j+k$.
We have that $\alpha_0$ can be any real number since it doesn't appear on the lhs and it cancels out in the rhs. Similar conclusions can be drawn for $\alpha_1$ and $\alpha_2$ after expanding some products.
For $\alpha_i\; i\ge 3$, we notice that on the lhs this parameter multiples mixed monomials ($x^jy^k\; j, k > 0$) that are not present in the rhs. Hence $\alpha_i = 0\; i\ge 3$.
This proof is quite informal but I hope you got the main idea. 
A: There's another way of looking at the problem, resulting in a characterization of all the possible solutions, with almost no appeal to differentiation and derivatives. Consider the following more general problem:

Find all pairs of functions $ f , g : \mathbb R \to \mathbb R $ satisfying
$$ f ( x ) - f ( y ) = ( x - y ) g ( x + y ) \tag 1 \label 1 $$
for all $ x , y \in \mathbb R $.

Putting $ y = 0 $ in \eqref{1} we get $ f ( x ) = x g ( x ) + f ( 0 ) $ for all $ x \in \mathbb R $. Using this, we can rewrite \eqref{1} only in terms of $ g $, as
$$ x g ( x ) - y g ( y ) = ( x - y ) g ( x + y ) \tag 2 \label 2 $$
for all $ x , y \in \mathbb R $. The functional equation \eqref{2} has already appeared on Mathematics Stack Exchange, for example in the posts "Determining all $ f : \mathbb R \to \mathbb R $ that satisfy $ x f ( x ) - y f ( y ) = ( x - y ) f ( x + y ) $" and "How can I find $ f ' ( 2 x ) $ given that $ x f ( x ) - y f ( y ) = ( x - y ) f ( x + y ) $?", and you can follow the methods presented there to find out $ g ( x ) = \bigl ( g ( 1 ) - g ( 0 ) \bigr ) x + g ( 0 ) $ (focus on those methods that use nothing about differentiability). Therefore, we've found that every pair of functions satisfying \eqref{1} must be of the form
$$ \begin {cases} g ( x ) = a x + b \\ f ( x ) = a x ^ 2 + b x + c \end {cases} \tag 3 \label 3 $$
for some constants $ a , b , c \in \mathbb R $. It's straightforward to verify that indeed any pair of functions of this form gives a solution of \eqref{1}.
Getting back to the original equation
$$ f ' \left ( \frac { x + y } 2 \right ) = \frac { f ( x ) - f ( y ) } { x - y } \text , \tag 0 \label 0 $$
defining $ g : \mathbb R \to \mathbb R $ with $ g ( x ) = f ' \left ( \frac x 2 \right ) $ for all $ x \in \mathbb R $, we get the equation \eqref{1}. Since any $ f $ of the form in \eqref{3} is differentiable and for $ g $ of the form in \eqref{3} we have $ f ' ( x ) = 2 a x + b = g ( 2 x ) $ for all $ x \in \mathbb R $, the extra required relation between $ g $ and $ f ' $ is automatically ensured, and thus \eqref{3} gives the general form of the solutions of \eqref{0}.
