PurposeThis study aims to discuss the numerical solutions of weakly singular Volterra and Fredholm integral equations, which are used to model the problems like heat conduction in engineering and the electrostatic potential theory, using the modified Lagrange polynomial interpolation technique combined with the biconjugate gradient stabilized method (BiCGSTAB). The framework for the existence of the unique solutions of the integral equations is provided in the context of the Banach contraction principle and Bielecki norm.Design/methodology/approachThe authors have applied the modified Lagrange polynomial method to approximate the numerical solutions of the second kind of weakly singular Volterra and Fredholm integral equations.FindingsApproaching the interpolation of the unknown function using the aforementioned method generates an algebraic system of equations that is solved by an appropriate classical technique. Furthermore, some theorems concerning the convergence of the method and error estimation are proved. Some numerical examples are provided which attest to the application, effectiveness and reliability of the method. Compared to the Fredholm integral equations of weakly singular type, the current technique works better for the Volterra integral equations of weakly singular type. Furthermore, illustrative examples and comparisons are provided to show the approach’s validity and practicality, which demonstrates that the present method works well in contrast to the referenced method. The computations were performed by MATLAB software.Research limitations/implicationsThe convergence of these methods is dependent on the smoothness of the solution, it is challenging to find the solution and approximate it computationally in various applications modelled by integral equations of non-smooth kernels. Traditional analytical techniques, such as projection methods, do not work well in these cases since the produced linear system is unconditioned and hard to address. Also, proving the convergence and estimating error might be difficult. They are frequently also expensive to implement.Practical implicationsThere is a great need for fast, user-friendly numerical techniques for these types of equations. In addition, polynomials are the most frequently used mathematical tools because of their ease of expression, quick computation on modern computers and simple to define. As a result, they made substantial contributions for many years to the theories and analysis like approximation and numerical, respectively.Social implicationsThis work presents a useful method for handling weakly singular integral equations without involving any process of change of variables to eliminate the singularity of the solution.Originality/valueTo the best of the authors’ knowledge, the authors claim the originality and effectiveness of their work, highlighting its successful application in addressing weakly singular Volterra and Fredholm integral equations for the first time. Importantly, the approach acknowledges and preserves the possible singularity of the solution, a novel aspect yet to be explored by researchers in the field.