This paper presents a new approach for the imposing various boundary conditions on radial basis functions and their application in pseudospectral radial basis function method. The various boundary conditions such as Dirichlet, Neumann, Robin, mixed and multi-point boundary conditions, have been considered. Here we propose a new technique to force the radial basis functions to satisfy the boundary conditions exactly. It can improve the applications of existing methods based on radial basis functions especially the pseudospectral radial basis function method to handling the differential equations with more complicated boundary conditions. Several examples of one, two, and three dimensional problems with various boundary conditions have been considered to show the efficacy and versatility of the proposed method.