This paper presents a novel numerical method based on the extended moving least squares finite difference method(MLS FDM) for solving 1-D Stefan problem. The MLS FDM is employed for easy numerical modelling of the moving boundary and Taylor polynomial is extended using wedge function for accurate capturing of interfacial singularity. Difference equations for the governing equations are constructed by implicit method which makes the numerical method stable. Numerical experiments prove that the extended MLS FDM show high accuracy and efficiency in solving semi-infinite melting, cylindrical solidification problems with moving interfacial boundary.