Simple, mesh/grid free, explicit and implicit numerical schemes for the solution of linear advection-diffusion problems is developed and validated herein. Unlike the mesh or grid-based methods, these schemes use well distributed quasi-random points and approximate the solution using global radial basis functions. The schemes can be seen as generalized finite differences with random points instead of a regular grid system. This allows the computation of problems with complex-shaped boundaries in higher dimensions with no need for complex mesh/grid structure and with no extra implementation difficulties