“…It is an initial value method even if we solve a different model for each iteration when the governing differential equation is not invariant under every scaling group of point transformation. Its versatility has been shown by solving several problems of interest: free boundary problems [33,23,28,29], a hyperbolic moving boundary problem [25], the Homann and the Hiemenz flows governed by the Falkner-Skan equation in [26], one-dimensional parabolic moving boundary problems [30], two variants of the Blasius problem [32], namely: a boundary layer problem over moving surfaces, studied first by Klemp and Acrivos [38], and a boundary layer problem with slip boundary condition, that has found application to the study of gas and liquid flows at the micro-scale regime [18,40], parabolic problems on unbounded domains [34] and, recently, see the preprints: [22] parabolic moving boundary problems, and [21] an interesting problem in boundary layer theory: the so-called Sakiadis problem [45,46].…”