Astronomers have been using the measured luminosity to estimate the {\em luminous mass} of stars, based on empirically established mass-to-light ratio which seems to be only applicable to a special class of stars---the main-sequence stars---with still considerable uncertainties. Another basic tool to determine the mass of a system of stars or galaxies comes from the study of their motion, as Newton demonstrated with his law of gravitation, which yields the {\em gravitational mass}. Because the luminous mass can at best only represent a portion of the gravitational mass, finding the luminous mass to be different or less than the gravitational mass should not be surprising. Using such an apparent discrepancy as a compelling evidence for the so-called dark matter, which has been believed to possess mysterious nonbaryonic properties and present a dominant amount in galaxies and the universe, seems to be too far a stretch when seriously examining the facts and uncertainties in the measurement techniques. In our opinion, a galaxy with star type distribution varying from its center to edge may have a mass-to-light ratio varying accordingly. With the thin-disk model computations based on measured rotation curves, we found that most galaxies have a typical mass density profile that peaks at the galactic center and decreases rapidly within $\sim 5%$ of the cut-off radius, and then declines nearly exponentially toward the edge. The predicted mass density in the Galactic disk is reasonably within the reported range of that observed in interstellar medium. This leads us to believe that ordinary baryonic matter can be sufficient for supporting the observed galactic rotation curves; speculation of large amount of non-baryonic matter may be based on an ill-conceived discrepancy between gravitational mass and luminous mass which appears to be unjustified.

An accurate computational method is presented to determine the mass distribution in a rotating thin-disk galaxy from given rotation curve by applying Newtonian dynamics for an axisymmetrically rotating thin disk of finite size with or without a central spherical bulge. The governing integral equation for mass distribution, resulting from the balance between the Newtonian gravitational force and centrifugal force due to rotation at every point on the disk, is transformed via a boundary-element method into a linear algebra matrix equation that can be solved numerically for rotation curves with a wide range of shapes. To illustrate the effectiveness of this computational method, mass distributions in several mature spiral galaxies are determined from their measured rotation curves. All the surface mass density profiles predicted by our model exhibit approximately a common exponential law of decay, qualitatively consistent with the observed surface brightness distributions. When a central spherical bulge is present, the total galactic mass increases only slightly but the mass distribution in the galaxy is altered in such a way that the periphery mass density is reduced while more mass appears toward the galactic center. By extending the computational domain beyond the galactic edge, we can determine rotation velocity outside the cut-off radius which appears to continuously decrease and gradually approach the Keplerian rotation velocity out over twice the cut-off radius. Am examination of the circular orbit stability suggests that galaxies with flat or increasing rotation velocities with radius are more stable than those with decreasing rotation velocities especially in the region near the galactic edge. Our results demonstrate the fact that Newtonian dynamics can be adequate for describing the observed rotation behavior of mature spiral galaxies.

By numerically solving the mass distribution in a rotating disk based on Newton's laws of motion and gravitation, we demonstrate that the observed flat rotation curves for most spiral galaxies correspond to exponentially decreasing mass density from galactic center for the most of the part except within the central core and near periphery edge. Hence, we believe the galaxies described with our model are consistent with that seen through the eyes of Newton. Although Newton's laws and Kepler's laws seem to yield the same results when they are applied to the planets in the solar system, they are shown to lead to quite different results when describing the stellar dynamics in disk galaxies. This is because that Keplerian dynamics may be equivalent to Newtonian dynamics for only special circumstances, but not generally for all the cases. Thus, the conclusions drawn from calculations based on Keplerian dynamics are often likely to be erroneous when used to describe rotating disk galaxies.

We present an efficient, robust computational method for modeling the Newtonian dynamics for rotation curve analysis of thin-disk galaxies. For a disk galaxy with a typical flat rotation curve, our modeling results show that the surface mass density monotonically decreases from the galactic center toward periphery, according to Newtonian dynamics. In a large portion of the galaxy, the surface mass density follows an approximately exponential law of decay with respect to the galactic radial coordinate. Yet the radial scale length for the surface mass density seems to be generally larger than that of the measured brightness distribution, suggesting an increasing mass-to-light ratio with the radial distance in a disk galaxy. In a nondimensionalized form, our mathematical system contains a dimensionless parameter which we call the "galactic rotation number" that represents the gross ratio of centrifugal force and gravitational force. The value of this galactic rotation number is determined as part of the numerical solution. Through a systematic computational analysis, we have illustrated that this galactic rotation number remains within $\pm 10%$ of 1.70 for a wide variety of rotation curves. This implies that the total mass in a disk galaxy is proportional to $V_0^2\,R_g$, with $V_0$ denoting the characteristic rotation velocity and $R_g$ the radius of the galactic disk. The predicted total galactic mass of the Milky Way is in good agreement with the star-count data.

Observations reveal that mature spiral galaxies consist of stars, gases and plasma approximately distributed in a thin disk of circular shape, usually with a central bulge. The rotation velocities quickly increase from the galactic center and then achieve a constant velocity from the core to the periphery. The basic dynamic behavior of a mature spiral galaxy, such as the Milky Way, is well described by simple models balancing Newtonian gravitational forces against the centrifugal forces associated with a rotating thin axisymmetric disk. In this research, we investigate the effects of adding central bulges to thin disk gravitational models. Even with the addition of substantial central bulges, all the critical essential features of our thin disk gravitational models are preserved. (1) Balancing Newtonian gravitational and centrifugal forces at every point within the disk yields computed radial mass distributions that describe the measured rotation velocity profiles of mature spiral galaxies successfully. (2) There is no need for gravity deviations or ``massive peripheral spherical halos of mysterious Dark Matter''. (3) The calculated total galactic masses are in good agreement with star count data. (4) The addition of central bulges increases the calculated total galactic masses, possibly more consistent with the presence of galactic gases, dust, grains, lumps, planets and plasma in addition to stars. (5) Compared with the light distribution, our mass distributions within the disk are larger out toward the galactic periphery which is cooler with lower opactiy/emissivity (and thus darker). This is apparent from edge-on views of galaxies which display a dark disk-line against a much brighter galactic halo.

For mature spiral galaxies, the rotation velocities quickly increase from the galactic center and achieve a constant velocity from the core to the periphery. This dynamic behavior is described by models balancing Newtonian gravitational and centrifugal forces in rotating thin axisymmetric disks. Freeman's disk assumes a mass density decreasing exponentially with radius which correctly produces rotational velocities which increase from the galactic center to a maximum near the outer core, but then decreases out to the periphery contrary to measurements. Mestel's disk assumes a mass distribution decreasing more slowly (inversely with radius) that yields a constant rotational velocity across the entire disk, but has an unrealistic central mass singularity and does not describe the core rotation properly. Thus combine the Freeman and Mestel disks to utilize their strengths and eliminate their deficiencies. Utilize the Freeman formula for the central core, and the Mestel formula beyond the core to the galactic rim. This combined model produces rotation curves comparable to the measurements of mature spiral galaxies. For a more general thin-disk model, we develop an alternative computational method to solve for the mass density distributions for various measured rotation curves We compute radial mass densities that balance the Newtonian gravitational and centrifugal forces at every point. The computational solutions show mass densities which decrease approximately exponentially in the central core (similar to Freeman), and transition to a slower inverse radial decrease (similar to Mestel) to the periphery. Thus these diverse approaches yield similar results and are mutaully self-consistent with Newtonian gravity/dynamics.

The measured rotation velocity profiles of mature spiral galaxies are successfully described with a gravitational model consisting of a thin axisymmetric disk of finte radius. The disk is assumed uniformly thin but with variable radial mass density. The governing integral equation is based on mechanical balance between Newtonian gravitational and centrifugal forces (due to galaxy rotation) at each and every point in a finite set of concentric rings. The nondimensionalized mathematical system contains a dimensionless parameter we call ``galactic rotation parameter'' which concisely crystallizes perspective. Computational solutions are obtained for the radial mass distributions that satisfy the measured rotational velocity profiles. Together with a constraint equation for mass conservation, the galactic rotational parameter is also determined from which the total galactic mass is calculated from measured galactic radii and maximum rotation velocities. These calculated total galactic masses are in good agreement with data. Our deduced exponentially decreasing mass distributions in the central galactic core are in agreement with almost all others. However our mass distributions differ toward the galactic periphery with more ordinary baryonic mass in these outer disk regions which are cooler with lower opactiy/emissivity (and thus darker).