make simulation code pyhon of a system in contact with a thermal reservoir at temperature T, the system is composed of N particles that can be found in two states with zero energy or greater than zero energy. calculate using Monte Carlo method with metropole algorithms to calculate average energy and specific heat