We present an algorithm for numerically computing an absolutely continuous
invariant measure associated with a piecewise C^2 expanding mapping S:A -> A
on a bounded region A in R^n. The method is based on the Galerkin projection
principle for solving an operator equation in a Banach space. With the help of
the modern notion of functions of bounded variation in multidimension, we
prove the convergence of the algorithm