Numerical solution of the soil freezing problem