“…Non-linear partial differential equations play a vital role in modeling real-life phenomena, especially in the field of applied sciences such as meteorology, aerospace engineering, oceanography, geology, astronomy, and many other disciplines. For the realistic treatment of these non-linear models, several numerical techniques such as the barycentric interpolation collocation method [1], Fourier spectral method [2], and reproducing kernel method [3] have been adopted for the numerical treatment of these non-linear partial differential equations, such as the Klein-Gordon-Zakharov equations, Schrödinger equation, Duffing systems, and many more. The Hunter-Saxton equation is one among those non-linear partial differential equations that describe waves of a nematic liquid crystal in an enormous director field and arise at the short-wave limit of the Camassa-Holm equationan integrable model of shallow-water waves propagating uni-directionally across a flat bottom [4].…”