A classical problem in quantum mechanics involving the Dirac Delta function is given by
$$
y''+(\delta(x)-\lambda^2)y=0.
$$ 
Then, to find ''bound states'', you solve on the right and find the converging solution as $x\rightarrow \infty$, then solve on the left similarly. Assume continuity of the solution $y$. The jump condition on $y'$ is found by integrating from $-\epsilon$ to $\epsilon$ and taking $\epsilon\rightarrow 0$. The only value of $\lambda$ giving a bound state is then found to be $1/2$. The corresponding state is continuous everywhere and differentiable everywhere except at zero.

My issue with this is that to make this solution mathematically correct, one has to make sense of the multiplication of a continuous non-differentiable function by the delta Dirac function. In a previous [question](https://math.stackexchange.com/questions/1317882/multiplication-of-a-distribution-by-a-continuous-function) of mine, Alan U. Kennington suggested to use Radon measures. However, I believe there must be something simpler that makes sense of the problem above by mathematically making sense of multiplying a continuous function by a Dirac delta function. I am not necessarily looking for a complete answer but rather a reference that discusses the problem.

The reason I am asking this question is because I am facing a third order equation with coefficients involving the Dirac delta function. I am able to find a solution but I would like to make my computations more mathematically sound. I already asked this question on Math Stack Exchange but did not get an answer.