Methods are developed for dealing with the various dynamical problems that arise because of convective zones in stars. A system of equations for stellar convection is derived from the full equations of compressible fluid dynamics with the aid of two major approximations. The first of these is the anelastic approximation, which involves both the filtering out of acoustic waves and a suitable linearization of the fluctuating thermodynamic variables. The second one approximates the horizontal structure of convection by expanding the motion in a set of horizontal cellular platforms and severely truncating the expansion. The resulting system of partial differential equations, referred to as the anelastic modal equations, is outlined along with suggested boundary conditions and techniques for solving the equations. Ways of assessing the overall validity of the present treatment are discussed.