Intense debate at work place around the solution for this:
Let $M \in C[a,b]$ be a continuous function on the closed interval $[a,b]$ that satisfies $$\int_{a}^{b}M(x)\eta^{\prime\prime}dx = 0,$$ for all $\eta \in C^{2}\left[a,b\right]$ satisfying $\eta(a) = \eta(b) = \eta^{\prime}(a) = \eta^{\prime}(b) = 0$. Prove that $M(x) = c_{0}+c_{1}x$ for suitable $c_{0}$ $c_{1}$. What can you say about $c_{0}$, $c_{1}$?
I tried to use integration by parts, and use the fundamental lemma of the calculus of variations, and the lemma of Du Bois and Reynolds to prove it, but that requires $M$ to be $C^1([a, b])$.
Define $M_\epsilon(x) = (M \ast \phi_\epsilon)(x)$ where $\phi_\epsilon$ is a standard mollifier supported on $[-\epsilon,\epsilon]$. Then for all $\eta \in C^\infty$ supported on $[a+\epsilon, b-\epsilon],$
$$\begin{align*} \int_{a+\epsilon}^{b-\epsilon} M_\epsilon''(x)\eta(x)\,dx &= \int_{a+\epsilon}^{b-\epsilon} M_\epsilon(x)\eta''(x)\,dx \\ &= \int_{a+\epsilon}^{b-\epsilon}\int_{-\epsilon}^\epsilon M(x-t)\phi_\epsilon(t)\eta''(x)\,dt\,dx \\ &= \int_{-\epsilon}^\epsilon \phi_\epsilon(t)\int_{a+\epsilon-t}^{b-\epsilon-t} M(x)\eta''(x+t)\,dx\,dt \\ &= \int_{-\epsilon}^\epsilon \phi_\epsilon(t)\int_a^b M(x)\eta''(x+t)\,dx\,dt \\ &= 0 \end{align*}$$
In particular, this shows that $M_\epsilon(x) = c_{0,\epsilon} + c_{1,\epsilon}x$ for $x\in [a+\epsilon, b-\epsilon].$
Since $M_\epsilon\to M$ pointwise on $(a,b)$ as $\epsilon \to 0^+,$ $M(x) = c_0 + c_1x$ for $a< x < b.$