Numerical simulation of random Dirac operators