Numerical solution of the inverse spectral problem for Bessel operators

We consider some inverse spectral problems associated with the singular Sturm-Liouville equation −u″+(q(x)+ℓ(ℓ+1)x2)u=λu0