Precision matrix expansion – efficient use of numerical simulations in estimating errors on cosmological parameters