  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.
