Uncertainty quantification in complex systems using approximate solvers