We plug in the pairs (a,x),(a,2x),(a+x,x) and (a−x,x) to get
f(a+x)+f(a−x)−2f(a)f(a+2x)+f(a−2x)−2f(a)f(a+2x)+f(a)−2f(a+x)f(a−2x)+f(a)−2f(a−x)=g(a)x2=4g(a)x2=g(a+x)x2=g(a−x)x2
respectively. Combining these equations in the form 2E1−E2+E3+E4 the left hand side vanishes, yielding an equation in g:(g(a+x)+g(a−x)−2g(a))x2=0, i.e.
g(a)=2g(a+x)+g(a−x)
Since g is continuous, it must be linear, i.e. g(x)=c1x+c0. However, the original equation for x=y together with the concavity condition now gives us
0⩾f(2x)+f(0)−2f(x)=(xc1+c0)x2
for all x, which is only possible if c1=0. Thus g(x)≡c0=2A is constant and
f(x+y)+f(x−y)−2f(x)=2Ay2
This suggests that f is a quadratic function, so we can set f(x)=Ax2+f1(x). Then (∗) becomes f1(x+y)+f1(x−y)−2f1(x)=0, so an easy induction gives us
f1(nx)−f1(0)=n(f1(x)−f1(0)) for all n∈Z
By setting f1(0)=C and f1(1)=B+C we obtain f1(x)=Bx+C and f(x)= Ax2+Bx+C for all x∈Q. By concavity of f we conclude that f(x)=Ax2+Bx+C for all real x.
## Remark.
In fact, (∗) implies that the second derivative of f is constant by taking y→0 and the problem is solved. The solution presented here avoids use of derivatives.