Two-step runge-kutta methods with quadratic stability functions