Meshless methods have some obvious advantages such as they do not require meshes in the domain and on the boundary, only some nodes are needed in the computation. Furthermore, for the boundary-type meshless methods, the nodes are even not needed in the domain and only distributed on the boundary. Practice shows that boundary-type meshless methods are effective for homogeneous problems. But for inhomogeneous problems, the application of these boundary-type meshless methods has some difficulties and need to be studied further. The hybrid boundary node method (HBNM) is a boundary-only meshless method, which is based on the moving least squares (MLS) approximation and the hybrid displacement variational principle. No cell is required either for the interpolation of solution variables or for numerical integration. It has a drawback of ‘boundary layer effect’, so a new regular hybrid boundary node method (RHBNM) has been proposed to avoid this pitfall, in which the source points of the fundamental solutions are located outside the domain. These two methods, however, can only be used for solving homogeneous problems. Combining the dual reciprocity method (DRM) and the HBNM, the dual reciprocity hybrid boundary node method (DRHBNM) has been proposed for the inhomogeneous terms. The DRHBNM requires a substantial number of internal points to interpolate the particular solution by the radial basis function, where approximation based only on boundary nodes may not guarantee sufficient accuracy. Now a further improvement to the RHBNM, i.e., a combination of the RHBNM and the multiple reciprocity method (MRM), is presented and called the multiple reciprocity hybrid boundary node method (MRHBNM). The solution comprises two parts, i.e., the complementary and particular solutions. The complementary solution is solved by the RHBNM. The particular solution is solved by the MRM, i.e., a sum of high-order homogeneous solutions, which can be approximated by the same-order fundamental solutions. Compared with the DRHBNM, the MRHBNM does not require internal points to obtain the particular solution for inhomogeneous problems. Therefore, the present method is a real boundary-only meshless method, and can be used to deal with inhomogeneous problems conveniently. The validity and efficiency of the present method are demonstrated by a series of numerical examples of inhomogeneous potential problems.