Abstract. In this work the method of approximate particular solutions using compactly supported kernels is investigated. In the work of [1] the globally supported radial kernels for the construction of solution are used, and it is observed that for large scaled PDEs the differentiation matrix is ill-conditioned. We extended the work of [1] for compactly supported kernels, in order to solve large-scaled PDEs engineering sciences. The numerical scheme of the present method of approximate particular solutions is very accurate and simple in implementation. Three benchmark problems are solved by the present numerical scheme and the results are compared to other methods in the literature.