I am trying to solve the equation $f(x+y)-f(x-y)=2f'(x)f'(y)$ for all $f:\mathbb{R}\to\mathbb{R}$ non-constant, differentiable functions.
Here is my progress:
Any solution must be an even function and therefore zero is a critical point. Further if $f(x)$ is a solution, then $f(x)+c$ is also a solution for any constant $c.$
In addition $\dfrac{x^2}2, -\cos x$ and $\cosh x$ are some solutions of our functional equation.
How can I solve this equation completely?
OP has already derived
$f(2x)-f(0)~=~2f^{\prime}(x)^2~\geq~0$.
$f^{\prime}(0)~=~0$.
$f(y)-f(-y)~=~0$ even function.
$x\mapsto f(x)$ solution $\quad\Rightarrow\quad$ $x\mapsto f(x)+c$ solution.
Assume that $f$ is not a constant function. Then there exists a point $y_0\in\mathbb{R}$ such that $f^{\prime}(y_0)\neq 0$. From the bootstrap/recurrence relation $$ f^{(n+1)}(x)~=~\frac{f^{(n)}(x+y_0)-f^{(n)}(x-y_0)}{2f^{\prime}(y_0)}, \qquad n~\in~\mathbb{N}_0,$$ we can deduce inductively that $f\in C^{\infty}(\mathbb{R})$ is infinitely many times differentiable.
Under the assumption that $f$ is 4 times differentiable, one can derive the ODE $$f^{(3)}(x)~=~f^{\prime}(x)f^{(4)}(0).$$
Combined with what OP already has derived, it follows that the full solution is