This paper deals with the numerical solution of second kind integral equations with fixed singularities of Mellin convolution type. The main difficulty in solving such equations is the proof of the stability of the chosen numerical method, being the noncompactness of the Mellin integral operator the chief theoretical barrier. Here, we propose a Nyström method suitably modified in order to achieve the theoretical stability under proper assumptions on the Mellin kernel. We also provide an error estimate in weighted uniform norm and prove the well-conditioning of the involved linear systems. Some numerical tests which confirm the efficiency of the method are shown
