Nonlinear reaction-diffusion equations are used to model many nonlinear phenomena. In this work we present a numerical method to the solve two-dimentional nonlinear reaction-diffusion equations. Meshless and collocation techniques using radial basis functions (RBFs) with the help of radial point Hermite interpolation (RPHI) method are used to construct basis functions so called shape functions. Due to the use of meshless method, no mesh generation is done in the spatial domain. Time discretization is performed using the finite difference method and Taylor expansion is used for the nonlinear part. The accuracy and efficiency of the method are checked with a numerical example. Results show that this procedure is stable through the time