Approximate pre-classical solutions in loop quantum cosmology

In this paper we introduce a numerical approximation technique to obtain pre-classical solutions to models of loop quantum gravity. In particular, we apply the technique to vacuum Bianchi I cosmological models and recover known solutions. We also present a pre-classical solution to the Bianchi I LRS model with cosmological constant, which has not appeared elsewhere.