In a recent research, the authors established an approximation to the modified error function (MEF) 𝜑 𝛿 for a small positive thermal conductivity coefficient 𝛿 > 0, leaving the problem open for the general case −1 < 𝛿 < ∞. The approximation represents an approximate solution to the Stephan problem with linear thermal conductivity in a semi-infinite body, which converges to the classical MEF as 𝛿 → 0 + . In this paper, a new approximation to the MEF is found for −1 < 𝛿 < ∞ and converges to the classical MEF as 𝛿 → 0 ± . A comparative analysis shows that our proposed approximation gives a better estimate than the one in this recent reasearch by Ceratani et al.