Surface temperature of an airless planet: day, night, and below the ground

Afterwards you can turn a heat equation with a radiating surface into ODEs for solve_ivp and compute an airless body's day, night, and subsurface temperatures.