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

Afterwards you can turn the heat equation of an airless body into ODEs for DifferentialEquations.jl and compute its day, night, and subsurface temperatures.