We have shown, that our implementation of the twodimensional Ising Model without the presence of a magentic field behaves, as expected from theoretical predictions. Especially showing relative errors in the order of \qtyrange{e-3}{e-4}{\percent} in quantities that are linear in the spin quantum numbers. Furthermore the statistical phenomena of phase change at the critical temperature has been proven. In future work it may be possible to improve the accuracy of the phase change behavior given more computing resources. We additionaly propose a parallelization approach for calculating the system evolution for multiple systems concurrently. We obtain a parallel efficiency of $\mathcal{E} = \num{.61}$ with the potential of further speed up under ideal conditions.
