Code attachment to article Semi-analytic benchmark for the Stefan problem - assessing accuracy of enthalpy-based methods.

Run with python3 stefan-benchmark dim, where dim stands for 1d, 2d, or 3d.