  In the work, the numerical methods are designed for the
Bogoliubov-Tolmachev-Shirkov model in superconductivity theory. The numerical
methods are novel and effective to determine the critical transition
temperature and approximate to the energy gap function of the above model.
Finally, a numerical example confirming the theoretical results is presented.
