Abstract
A wavelet-Galerkin procedure is derived and implemented for the numerical solution of inhomogeneous diffusion equations subject to mass specification involving non-polynomial functions. It is demonstrated that the accuracy of approximation of these functions by polynomial interpolation at Chebyshev nodes is in good agreement with the exact solution for several numerical experiments.