Abstract
The stellar hydrodynamical equations give the components of stress required to support a stellar system against gravitational collapse. The paper derives a numerical algorithm to solve the stellar hydrodynamical equations in spheroidal coordinates, widely used for galactic modelling.