Two-step Runge-Kutta methods with quadratic stability functions