Abstract.A meshless local method of approximated particular solutions (LMAPS) is used to analyze problem described by the convectiondiffusion equation. The method solves the steady convection-diffusion equation with reaction term. The discretized system of equations is derived via interpolation procedure using integrated multiquadrics (MQ) radial basis functions (RBF). The solution of the equation is performed over simple geometry with non-uniform velocity field and results are presented in the article. The LMAPS method is capable to produce stable solutions with results comparable to the analytical solutions.