Self-propelled colloidal swimmers move by pushing the adjacent fluid backwards. The resulting motion of an asymmetric body depends on the profile of pushing velocity over its surface. We describe a method of predicting the motion arising from arbitrary velocity profiles over a given body shape, using a discrete-source "stokeslet" representation. The net velocity and angular velocity is a sum of contributions from each point on the surface. The contributions from a given point depend only on the shape. We give a numerical method to find these contributions in terms of the stokeslet positions defining the shape. Each contribution is determined by linear operations on the Oseen interaction matrix between pairs of stokeslets. We first adapt the Lorentz Reciprocal Theorem to discrete sources. We then use the theorem to implement the method of Teubner to determine electrophoretic mobilities of nonuniformly charged bodies.