“…The values of M 3 and P 3 are determined analytically. Whereas for the HT model (9) , the Euler corrector algorithm contains four equations to determine x 3 , y 3 , T 3 and θ 3 . In this case, M 3 , ρ 3 and P 3 are determined analytically, and for the PG model (3) , the Euler corrector algorithm always contains four equations to determine, in this case, x 3 , y 3 , M 3 and θ 3 .…”