Abstract: A new three-dimensional fundamental solution to the Stokes flow was proposed by transforming the solid harmonic functions in Lamb’s solution into expressions in terms of the oblate spheroidal coordinates. These fundamental solutions are advantageous in treating flows past an arbitrary number of arbitrarily positioned and oriented oblate spheroids. The least squares technique was adopted herein so that the convergence difficulties often encountered in solving three-dimensional problems were completely avoided. The examples demonstrate that present approach is highly accurate, consistently stable and computationally efficient. The oblate spheroid may be used to model a variety of particle shapes between a circular disk and a sphere. For the first time, the effect of various geometric factors on the forces and torques exerted on two oblate spheroids were systematically studied by using the proposed fundamental solutions. The generality of this approach was illustrated by two problems of three spheroids.