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