Computing the spectral function for singular Sturm-Liouville problems

Algorithms for computing Sturm-Liouville spectral density functions are developed based on several mathematical characterizations. Convergence and error bounds are derived and methods are tested on several examples. The results are compared with those from the existing software package SLEDGE.