Summary: 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 pre-classical solutions to the Schwarzschild interior geometry and the Bianchi I LRS model with cosmological constant, which have not appeared elsewhere.