In this paper, a local meshless method (LMM) based on radial basis functions (RBFs) is utilized for the numerical solution of various types of PDEs, due to the flexibility with respect to geometry and high order of convergence rate. In case of global meshless methods, the two major deficiencies are the computational cost and the optimum value of shape parameter. Therefore, research is currently focus towards localized radial basis function approximations, as the local meshless method is proposed here. The local meshless procedures is used for spatial discretization whereas for temporal discretization different time integrators are employed. The proposed local meshless method is testified in terms of efficiency, accuracy and ease of implementation using regular and irregular domains.