  This paper is devoted to derivation of plasma equations, considering chemical
reactions and radiations. After deriving the governing equations of
multi-particle plasmas, the ordinary differential equations of plasma shock
layer equations are derived. These ODEs are in the form of a singularly
perturbed systems of ordinary differential equations. Many of physical
environments consist of two or more types of particles, each having its own
characteristics, such an environment may undergo some evolutions in the
presence of external forces, e.g., in the presence of electromagnetic or
gravitation fields, and may also be affected by a source of energy or by
interaction with its external environment. The situation becomes more complex
when there are some internal interactions due to collisions between particles,
chemical reactions and various dissipation mechanisms. From microscopic point
of view, the random motion of different particles and local forces, energy or
mass transaction and their recursive dependence with that of global makes the
qualitative behaviors of the environment even more complex. Heat convection in
a liquid mixture, a vibrating or rotating viscose fluid, chemically reacting
multiphase media and plasma in an electric and magnetic field are examples, to
say a few. Up to the author knowledge there is no general and precise
mathematical formulation efficiently describing such environments, although
statistical formulation in Kinetic theory combined with some physical
interpretations is a good approach but insufficient, say in describing
radiation. Our aim in this note is to give a set of partial differential
equation governing evaluation of a multispecies environment consisting of
charged particles in an electromagnetic field.
