Sample Fortran code for computing Ising RDMs based on the stochastic series expansion algorithm. Input the system size, transverse field strength, subregion size, thermalization steps, and Monte steps. Return the reduced density matrices of the subregion.