Newtonian analytical studies of spiral galaxy rotation have most often considered the gravitational force imparted to a test body (whose rotational velocity is measured) to be determined by the summation of all mass contained within its orbit, diminished by the inverse-square of its radial distance.
This extends, to planar disk mass distributions, shell theorem's justification for treating discrete spherical masses as though their mass was concentrated at a central point - thereby allowing the application of the inverse-square, two body point-mass law of universal gravitation. Obviously however, a spiral galaxy is not a spherically symmetrical distribution of mass.
http://en.wikipedia.org/wiki/Law_of_universal_gravitation#Bodies_with_spatial_extent states:
"If the bodies in question have spatial extent (rather than being theoretical point masses), then the gravitational force between them is calculated by summing the contributions of the notional point masses which constitute the bodies. In the limit, as the component point masses become "infinitely small", this entails integrating the force (in vector form, see below) over the extents of the two bodies.
"In this way it can be shown that an object with a spherically-symmetric distribution of mass exerts the same gravitational attraction on external bodies as if all the object's mass were concentrated at a point at its centre. (This is not generally true for non-spherically-symmetrical bodies.)"
This would seem to indicate that solutions to a representative distribution of point masses should be obtained, then combined by vector summation to determine a more representative result. This may be necessary because the inverse-square force contributed by millions of discrete masses, each with individual separation distances much less than the subject's radial distance from the galactic center, would exceed the force estimation using a single separation distance term.
Also see http://en.wikipedia.org/wiki/Shell_theorem - for points outside the 'shell'.