The pertinent coupled equations describing the melting of granular polymeric solids on a hot moving interface were solved approximately with collocation method. A non-Newtonian boundary layer thickness function was defined, and the melting rate was evaluated in terms of this non-Newtonian boundary layer and the averaged velocity component in the moving direction. The results of this analysis were comparable to those by numerical simulation, successive iteration, and experimental data.