In this paper, a computational method is presented to solve potential-type Fredholm integral equations of the first kind, equations in which the unknown functions are singular at the endpoints of the integration domain, in addition to the weakly singular logarithmic kernels. This method provides a numerical solution based on the Newton interpolation technique via the Vandermonde matrix, which can accommodate an approximation of the unknown function, in such a manner that its singularity is easily removed, as well as the removal of kernel singularity. In addition, the Gauss–Legendre formula is adapted and applied for the computations of the obtained convergent integrals. Thus, the obtained numerical solution is equivalent to the solution of an algebraic equation in matrix form without applying the collocation method. The numerical solutions of the illustrated example are strongly converging to the exact solution for all values of 1 x including the end-points 1 whereas the exact solution fails to find the functional values at these end-points; which ensures the powerful and high accuracy of the presented computational technique