Quantum mechanics and gravitation theory are unified here, when the unnormalized Schroedinger wave function of an isolated typical elementary particle, is given dimensions of potential energy, and is assumed to represent part of the particle's gravitational self potential energy, as well as representing the particle itself. This assumption is shown to be consistent with the normalization of the particle's wave function and its probabilistic interpretation. It leads directly to the derivation of a set of covariant partial differential equations which couple quantum mechanics, general relativity and electroweak and strong physics, and explains how particle rest mass arises quantum mechanically. The point of view taken here is that a necessary condition that a theory of elementary particle physics be fundamental, is that its defining equations be differential equations, exclusively. (For instance the electroweak field equations are the Bianchi identities satisfied by the gravitational field equations and lack parity invariant solutions). Gravitation theory is thereby inserted into quantum mechanics, and gravitation theory(including general relativity) and quantum mechanics are shown to be compatible theories. Interacting particle systems of arbitrary complexity are represented here as the states occupied by a single gravitationally self interacting particle. If this single gravitationally self interacting particle is composed of matter and not antimatter, then a consequence of the theory is matter/antimatter asymmetry. Since the theory involves the gravitational self interaction of a single particle, any multiparticle system, can be represented in a background space which is ordinary four dimensional space time. This contrasts with conventional Schroedinger theory where the background space has 3N spatial dimensions and time, where N is the number of particles in the system. This simplification may lead to the numerical solution of complex molecular interactions useful in drug, materials and energy research. Solutions of the aforementioned covariant equations are shown to represent the following five categories of elementary particle states: 1. 3 spin 1/2 lepton states, 2. 3 spin 1/2 lepton states representing the antiparticles of the particles in category 1, 3.6 spin 1/2 quarks, 4. 3 spin 1 bosons, and 5. 1 spin 0 boson. The strength of the gravitational self interaction at short range for leptons is shown to be proportional to 1/(r squared) and to be 1/(r to the fourth) for quarks, where r is the distance from the particle. The rotational stability of the galaxy based on elementary particle gravitational self interaction is shown.