In numerical mathematics, the boundary knot method (BKM) is proposed as an alternative boundary-type meshfree distance function collocation scheme. Recent decades have witnessed a research boom on the meshfree numerical PDE techniques since the construction of a mesh in the standard finite element method and boundary element method is not trivial especially for moving boundary, and higher-dimensional problems. The boundary knot method is different from the other methods based on the fundamental solutions, such as boundary element method, method of fundamental solutions and singular boundary method in that the former does not require special techniques to cure the singularity. The BKM is truly meshfree, spectral convergent (numerical observations), symmetric (self-adjoint PDEs), integration