On the use of the geometric approach to global stability for three dimensional ODE systems: a bilinear case