“…The most attractive feature of this methodology is that all manipulations only involve the solutions of linear algebraic equations and the expansions can be carried out to any order in h on a symbolic manipulation platform. Also, our treatment of both the coating and half-space can be easily generalised to deal with anisotropic and/or inhomogeneous materials (Spencer, 1984(Spencer, , 2005, nonlinear effects (Erbay, 1997), liquid crystal elastomers (Goriely & Mihai, 2021;Liu et al, 2021), and the case when U and V depend on both of the in-plane coordinates.…”