A method is presented to reduce the singular Lippmann-Schwinger integral equation to a simple matrix equation. This method is applied to calculate the matrix elements of the reaction and transition operators, respectively, on the real axis and on the complex plane. The phase shifts and the differential scattering amplitudes are computable as well as the differential cross sections if the R- and/or T-matrix elements on the energy-shell are known. The method is applicable by using the Gaussian quadratures based on the Legenre, Laguer Chebyshev and shifted Chebyshev polynomials. Choosing the nodal points and weight functions depends on the aspects of the problem.