An algorithm for describing operator structures and calculating energy levels for hydrogen-like atoms in the quantum mechanics with nonnegative distribution function (Kuryshkin's quantum mechanics) is developed. The algorithm is implemented in Maple. An original library of necessary functions and quantities (Coulomb wave functions, Sturm functions, their Fourier images and kernels of integral transformations, Clebsch-Gordan coefficients, products of spherical harmonics, etc.) has been created. In accordance with the quantization rules in Kuryshkin's quantum mechanics, operators of observable quantities are computed. Based on the Hamiltonian obtained, energy levels of a hydrogen-like atom are calculated by the Ritz method (in the basis of Sturm functions). Numerical values of model parameters are determined by comparing results obtained with experimental data for hydrogen and alkali metals taken from the NIST Atomic Spectra Database Levels Data library. © Nauka/Interperiodica 2007.