We introduce a new approximation method for the distribution of functions of random variables that are real-valued. The approximation involves moment matching and exploits properties of the class of normal inverse Gaussian distributions. In the paper we examine the how well the different...