bem_solve

代码说明：

  This function solves the boundary integral equation (2.33) in the notes by the simplest boundary element method in the 2D case (though the code applies in 3D if the function G is redefined to be the 3D Green s function).