“…The nonclassical method, first introduced by Bluman and Cole [8] (see, for example, [2] or [3]), seeks invariance of a given partial differential equation (PDE) augmented with the invariant surface condition. As the determining equations for these nonclassical symmetries are nonlinear, there seemed to be little hope for this new method; however, with the development of computer algebra systems, the nineties saw a huge explosion of interest as several authors took interest in the nonclassical method and continues today to be an active area of interest (e.g., [9][10][11][12][13][14][15][16][17][18][19][20][21][22][23] and references within).…”