A jump process with finite variation and infinite activity is approximated by replacing the small jumps by their mean and a suitably scaled Brownian motion. Integration by parts for the approximated process with respect to the Brownian motion is investigated by numerical experiments. In particular, a Monte Carlo method, involving integration by parts for computation of a sensitivity measure Delta of a European put option in models with Normal Inverse Gaussian log returns, is applied.