  A numerical method for free boundary problems for the equation \[
u_{xx}-q(x)u=u_t \] is proposed. The method is based on recent results from
transmutation operators theory allowing one to construct efficiently a complete
system of solutions for this equation generalizing the system of heat
polynomials. The corresponding implementation algorithm is presented.
