A new purely meshless method for solving elliptic partial differential equations (PDEs) is presented. The method is based on the principle of the analog equation of Katsikadelis, hence its name meshless analog equation method (MAEM), which converts the original equation into a simple solvable substitute one of the same order under a fictitious source. The fictitious source is represented by multiquadric radial basis functions (MQ-RBFs). The integration of the analog equation yields new RBFs, which are used to approximate the sought solution. Then inserting the approximate solution into the PDE and boundary conditions (BCs) and collocating at the mesh-free nodal points results in a system of linear equations, which permit the evaluation of the expansion coefficients of the RBFs series. The method exhibits key advantages over other RBF collocation methods as it is highly accurate and the coefficient matrix of the resulting linear equations is always invertible. The accuracy is achieved using optimal values for the shape parameters and the centers of the multiquadrics as well as of the integration constants of the analog equation, which are obtained by minimizing the functional that produces the PDE. Without restricting its generality, the method is illustrated by applying it to the general second order 2D and 3D elliptic PDEs. The studied examples demonstrate the efficiency and high accuracy of the developed method.