Consider the following paper:

"A note on global stability for a tuberculosis model" by Gao and Huang: https://doi.org/10.1016/j.aml.2017.05.004

The methodology is understood in this paper apart from some technical things; How did the author arrive at constants $m_i, i=1,...,5$? They achieved this by killing the variable terms ($xu, xv, z, u, v, yu, yv$), however when I try solving the system, it outputs all $m_i$'s equalling zero. Can someone show a detailed answer how they arrived at these constants and how they got their final form of the Lyapunov function (page 168)?

Code:

    sol = Solve[{-m1 \[Beta] is + m3 l \[Beta] ss is/ls + 
         m4 (1 - l) \[Beta] ss == 
        0, -m1 \[Rho]1 \[Beta] ts + m3 \[Rho]1 l \[Beta] ss ts/ls + 
         m4 (1 - l) \[Beta] \[Rho]1 ss ts/is == 
        0, -(\[Mu] + \[Delta]) m3 + m4 \[Delta] ls/is == 0, 
       m1 \[Beta] is + m2 \[Rho]2 \[Beta] is - 
         m4 (\[Mu] + \[Alpha] + \[Gamma]) + m5 (\[Mu] + \[Rho]) == 0, 
       m1 \[Rho]1 \[Beta] ts + m2 \[Rho]1 \[Rho]2 \[Beta] ts + 
         m3 \[Rho] ts/ls - m5 (\[Mu] + \[Rho]) == 
        0, -m2 \[Rho]2 \[Beta] is + m3 \[Rho]2 \[Beta] vs is/ls == 
        0, -m2 \[Rho]1 \[Rho]2 \[Beta] ts + 
         m3 \[Rho]1 \[Rho]2 \[Beta] vs ts/ls == 0 }, {m1, m2, m3, m4, m5}]

EDIT:

Trying with the authors co-efficients, I managed to get:

    m1 = ss;
    m2 = (\[Delta] vs )/(((1 - l) (\[Mu] + \[Delta]) + \[Delta] l));
    m3 = (\[Delta] ls )/(((1 - l) (\[Mu] + \[Delta]) + \[Delta] l));
    m4 = ((\[Mu] + \[Delta]) is )/(((1 - 
             l) (\[Mu] + \[Delta]) + \[Delta] l));
    m5 = (1/(\[Mu] + \[Rho])) (\[Rho]1 \[Beta] ss ts + \[Rho]1 \[Rho]2 \
    \[Beta] ((\[Delta] vs ts)/((1 - 
                 l) (\[Mu] + \[Delta]) + \[Delta] l)) + ((\[Rho] \[Delta] \
    ts)/((1 - l) (\[Mu] + \[Delta]) + \[Delta] l)));
----------
$xu$ term:

    FullSimplify[-m1 \[Beta] is + m3 l \[Beta] ss is /ls + 
      m4 (1 - l) \[Beta] ss]

    (* 0 *)

$xv$ term:

    FullSimplify[-m1 \[Rho]1 \[Beta] ts + m3 \[Rho]1 l \[Beta] ss ts/ls + 
      m4 (1 - l) \[Beta] \[Rho]1 ss ts/is]

    (* 0 *)

$z$ term:

    FullSimplify[-(\[Mu] + \[Delta]) m3 + m4  \[Delta] ls/is]

    (* 0 *)

$u$ term:

        FullSimplify[
         m1 \[Beta] is + m2 \[Rho]2 \[Beta] is - 
          m4 (\[Mu] + \[Alpha] + \[Gamma]) + m5 (\[Mu] + \[Rho])]

    (* (-is \[Alpha] (\[Delta] + \[Mu]) - 
     is (\[Gamma] + \[Mu]) (\[Delta] + \[Mu]) + 
     is ss \[Beta] (\[Delta] + \[Mu] - l \[Mu]) - (-1 + 
        l) ss ts \[Beta] \[Mu] \[Rho]1 + is vs \[Beta] \[Delta] \[Rho]2 + 
     ts \[Delta] (\[Rho] + \[Beta] \[Rho]1 (ss + 
           vs \[Rho]2)))/(\[Delta] + \[Mu] - l \[Mu]) *)

$v$ term:

    FullSimplify[
     m1 \[Rho]1 \[Beta] ts + \[Rho]2 \[Beta] \[Rho]1 ts m2 + 
      m3 \[Rho] ts/ls - m5 (\[Mu] + \[Rho])]

    (* 0 *)

$yu$ term:

    FullSimplify[m3 \[Rho]2 \[Beta] vs is/ls - \[Rho]2 \[Beta] is m2]

     (* 0 *)

$yv$ term:

    FullSimplify[
     m3 \[Rho]2 \[Rho]1 \[Beta] vs ts/ls - \[Rho]2 \[Rho]1 \[Beta] ts m2]

     (* 0 *)

The $u$ term is causing some problems here, why?