The main aim of this paper is to numerically solve the first kind linear Fredholm and Volterra integral equations by using Modified Bernstein–Kantorovich operators. The unknown function in the first kind integral equation is approximated by using the Modified Bernstein–Kantorovich operators. Hence, by using discretization, the obtained linear equations are transformed into systems of algebraic linear equations. Due to the sensitivity of the solutions on the input data, significant difficulties may be encountered, leading to instabilities in the results during actualization. Consequently, to improve on the stability of the solutions which imply the accuracy of the desired results, regularization features are built into the proposed numerical approach. More stable approximations to the solutions of the Fredholm and Volterra integral equations are obtained especially when high order approximations are used by the Modified Bernstein–Kantorovich operators. Test problems are constructed to show the computational efficiency, applicability and the accuracy of the method. Furthermore, the method is also applied to second kind Volterra integral equations.