Main goal of the work was development of the software aimed to compute eigenvalues and eigenfunctions of regular and singular self-adjoin Sturm-Liouville problems, with both separated and coupled boundary conditions, and to approximate the continuous spectrum in the singular case. Computational algorithms of SLEIGN2 package [1] were taken as a basis for our software. SLEIGN2 code itself is based on ideas and methods of the original SLEIGN code [2].