Consider the following Sturm Liouville problem on an interval $[a,b]$

$$\frac{\mathrm{d}}{\mathrm{d} x}\left[p(x) \frac{\mathrm{d} y}{\mathrm{d} x}\right]+q(x) y=-\lambda w(x) y$$
for given coefficient function $p(x),q(x),w(x)>0$. Here we can assume they are smooth. One can consider two sets eigenvalues
$$0<\lambda_{1}^D<\lambda_{2}^D<\lambda_{3}^D<\cdots<\lambda_{n}^D<\cdots \rightarrow \infty$$
$$\lambda_{1}^N<\lambda_{2}^N<\lambda_{3}^N<\cdots<\lambda_{n}^N<\cdots \rightarrow \infty$$
where $\lambda_i^D$ are Dirichlet boundary condition and $\lambda_i^N$ are Neumann boundary condition. It is quite easy to show $\lambda_i^N\leq \lambda_i^D$ by the variational characterizations of the eigenvalues. However, it seems that under some condition of $p,q,w$, we will have
$$\lambda_{i+1}^N\leq \lambda_i^D$$
I encounter this when reading some papers, and it says this is well known. Does anyone know a good reference for this type of inequality? Who first found this? Do you have a simple proof?

PS. I googled a little bit and found many works on eigenvalue inequality of this type for Laplacian operator on dimension greater than 1. I am just interested in one dimension and here the operator is more general than Laplacian.