We have developed easy to use fast multipole method (FMM) libraries for the Laplace, low-frequency Helmholtz, and Stokes equations in two and three dimensions. The codes are based on a new method for applying translation operators and provide reasonable performance on either single core processors, or small multi-core systems using OpenMP.