Implementation of an X-FEM solver for the classical two-phase Stefan problem