Expansion  of Horvay-Cahn Solution
Procedure:
Set b = P e
Re-express Horvay-Cahn solution in parabolic coordinates
Expand in powers of e for fixed value of P
Find the thermal field T(x,h,f,e), interface shape x = f(h,f,e), and
Stefan number S(e) as functions of e through 2nd order