Cubic Hermite interpolation ensures continuity of derivatives between elements, so as well as providing an excellent representation of a smooth function, it also accurately models its gradients. Another attractive feature of this interpolation is that because each node is shared by neighbouring elements, a mesh of cubic Hermite elements has no more unknowns than a mesh of quadratic Lagrange ele...