In this work we apply a meshless method based on Local Boundary Integral Equations (LBIE) to find the solution to boundary value problems. We discretize the weak form through the use of special basis functions that, unlike the Finite Element Method (FEM), are not confined to an element and do not need the support of an underlying mesh. The approach developed can be applied to general 3D scalar boundary value problems that arise in areas such as electrostatics and acoustic scattering, among others.