I ran your numbers with
$$H(s) = \frac{2.123\times 10^6 \times (0.15s +1)}{s^3+12566.37s^2}$$ and unity negative feedback.
A system is stable if its parts are stable and all of the poles in its transfer function have negative real parts.
If a system has simple, unity negative feedback and the loop transfer function is \$P(s)/Q(s)\$, then the closed loop transfer function is
$$F(s) = \frac{\frac{P(s)}{Q(s)}}{1+\frac{P(s)}{Q(s)}} = \frac{P(s)}{P(s)+Q(s)}$$
The poles of \$F(s)\$ are then just the zeros of \$P(s)+Q(s)\$.
Your loop transfer function, with \$K_{vp}=0.15\$, and \$K_{vi}=1\$
is
$$\frac{P(s)}{Q(s)}=\frac{2.123\times 10^6 \times (0.15s +1)}{s^3+12566.37s^2}$$
If I have done my math correctly, we thus want to find the roots of
$$P(s)+Q(s) = (2.123\times 10^6 \times (0.15s +1))+(s^3+12566.37s^2)$$ $$P(s)+Q(s) = s^3 + 12566.37s^2 + 318450s + 2123000$$
The zeros of the above expression turn out to be
- r = -12541
- r = -12.6896 - 2.87376 j
- r = -12.6896 + 2.87376 j
Since the all of the roots have real parts which are negative, the closed loop system is stable. This is the definitive test for stability when either the loop gain, or closed loop transfer function is known in the form of a rational function. The Nyquist method works, but to be accurate, one still needs to know the number of right-half plane zeros, and the number of encirclements of -1 in each direction. The Bode plot method works with simple systems, but the gain margin test requires looking at all frequencies where the phase is \$180^{\circ}\$ (or \$0^{\circ}\$, depending). So, to get a definitive answer, when you know the transfer function, calculate as I have done, and rely on that, rather than on gain margin.
The key word in the text is "may"
In summary, if a negative feedback circuit has a loop gain that satisfies two condictions $$|H(j\omega_0)| \ge 1$$ $$\angle H(j\omega_0) = 180^{\circ}$$ then the circuit may oscillate at \$\omega_0\$.
(emphasis added).
Note: Roots of a polynomial may be found using online calculators such as Wolfram Alpha. In the days of slide rules, finding roots of higher order polynomials was not trivial. This fact is a contributing factor to the spread of graphical methods for testing stability. Another obvious factor is that we often don't know an accurate representation of the transfer function as a rational function, but can obtain Bode plots through simulation and/or Vector Network Analyzers.