Exact solution of inverse one phase Stefan problem is represented in the form of linear combination of integral error functions. Heat flux function is reconstructed and coefficients of solution function are found exactly. Test problem was considered for engineering purposes and it was shown that by collocation method the error for three points does not exceed 0.01%. Error estimate was calculated by maximum principle