The integral method of incremental hole-drilling is used extensively to determine the residual stress distribution in isotropic materials. When used with Tikhonov regularization, the method is robust and produces accurate results with minimal uncertainty. Alternatively, an optimal hole depth distribution can be found using the method of Zuccarello to improve the conditioning of the calibration matrices. If substantial measurement noise or a steep variation in stress exists, however, considerable uncertainty in, or distortion of, the calculated residual stress distribution can occur. Series expansion offers an alternative solution, but it has been reported to become unstable before meaningful accuracy can be achieved. Investigate the use of series expansion to determine a rapidly changing throughthickness residual stress distribution in an aluminium alloy 7075 plate subjected to laser shock peening treatment. Power series expansion of eigenstrains is used in finite element modelling to calculate the calibration coefficients. Monte Carlo simulation is used to determine robust uncertainties in the residual stress distributions. This allows the series order with the lowest RMS uncertainty in stress to be selected from those series orders that have converged. The best estimate of the residual stress distribution is thereby obtained. Series expansion is shown to be stable up to 8th order and convergence to a stress solution can be found before instability dominates. The method is insensitive to measurement errors due to the least-squares approach employed by the inverse solution. The use of series expansion reduces the RMS uncertainty in stress when compared to the regularized integral and Zuccarello methods.