This book provides practical demonstrations of how to carry out definite integrals with Monte Carlo methods using Mathematica. Random variates are sampled by the inverse transform method and the acceptance-rejection method using uniform, linear, Gaussian, and exponential probability distribution functions.