Presents a survey of intrinsic techniques and practical implementation of higher-order finite element schemes. This book addresses automatic hp-adaptivity in one and two spatial dimensions. It includes a discussion of suitable data structures and efficient linear solvers. It is suitable for both applied mathematicians and engineers.