The regular integrals that arise are approximated by a quadrature rule defined on a triangle. Laursen and Gellert [21] contains a selection of GaussLegendre quadrature rules for the triangle. The nonregular integrals that arise are computed by the following methods. See Jaswon and Symm [14], Terai [24], Banerjee and Butterfield [2] for the background to these methods.
The M_{0} and M_{0}^{t} operators have regular kernels, hence the aim is to find expressions for:
 (45) 
 (46) 
The integrals (45) and (46) may be written in the form:




In this section the subroutine L3LC.FOR is introduced and demonstrated through the test problem L3LC_T.FOR .
The program in file L3LC_T.FOR is a test for subroutine L3LC. It computes the solution to the Laplace problem interior to a cube. In order to use L3LC, the cube is represented by 24 uniform triangles. The boundary functions are approximated by a constant on each triangular panel.
Simple changes in the program allow the set up of the following data: the choice of quadrature rule, the geometry of the surface (of the approximating polyhedron), the boundary condition and the exact solution, the points in the interior where the solution is sought.
During execution, the program gives the solution at the collocation points (the points at the centre of each element) and the solution at the selected interior point. The program also give the exact solution at the same points so that computed and exact solutions may be compared.
In L3LC_T.FOR results are computed for the test problem of a cube, bounded by x,y,z = ±1. The cube is approximated by 24 uniform elements. The integrals are computed by a 7 point Gaussian quadrature rule for triangles.
The analytic problem that is considered is f(p) = p_{3}. The computed and exact results are compared at the collocation points in the output to the program. The solutions are also compared at (0,0,0.5); the exact solution is 0.5, the computed solution is 0.490392 to six decimal places.