Implicit Runge–Kutta methods have a special role in the numerical solution of stiff problems, such as those found by applying the method of lines to the partial differential equations arising in physical modelling. Of particular interest in this paper are the high-order methods based on Gaussian quadrature and the efficiently implementable singly implicit methods.