Implements ultraspherical spectral methods for solving linear ODEs