Cubic-coordinate reduction of a singular Sturm-Liouville operator ($L=-x^{-4}{D}^2+2x^{-5}D-a⟼-D_z^2-a,z=x^3/3$)

For $0<x<1$ and zero endpoint values, write $L=w^{-1}[-(p{y}^{\prime }{)}^{\prime }+qy]$ with $p=x^{-2}$, $w=x^2$ and $q=-a{x}^2$. The change of variable $z=x^3/3$ has $dz=x^{2}dx$ and sends $y^{\prime }=x^2{Y}_z$, $y^{\prime \prime }=2x{Y}_z+x^4{Y}_{zz}$. Thus it is an isometry from the weighted inner product to the ordinary integral on $0<z<1/3$, and the apparent singular operator becomes $-Y_{zz}-aY$. Its Dirichlet eigenvalues are $9n^2{\pi }^2-a$ and its normalized eigenfunctions are $\sqrt{6}\sin (n\pi x^3)$, $n\ge 1$. The singular endpoint is handled through this transformed regular domain; a zero eigenvalue occurs exactly when $a=9n^2{\pi }^2$.