In this paper, an algorithm will be presented that enables solving the two-phase inverse Stefan problem, where the additional information consists of temperature measurements in selected points of the solid phase. The problem consists in the reconstruction of the function describing the heat transfer coefficient, so that the temperature in the given points of the solid phase would differ as little as possible from the predefined values. The featured examples of calculations show a very good approximation of the exact solution and stability of the algorithm.