Geometric properties of lattice quantum gravity in two dimensions are studied numerically via Monte Carlo on Euclidean Dynamical Triangulations. A new computational method is proposed to simulate gravity coupled with fermions, which allows the study of interacting theories on a lattice, such as non-Riemannian gravity models. This was tested on Majorana spinors, where we obtained a Hausdorff dimension d(W) = 4.22 +/- 0.03, consistent with the bounds from the literature 4.19 < d(H) < 4.21.