7.This ϕ0ðtÞ replaces the essisrepeated to get new m_ 0fðtÞ. If jm_ 0f
ðtÞm_ f ðtÞjnΔtr5 105 kg go to step 8.Else, change ϕ0ðtÞ ¼ ϕðtÞ½ϕ0ðtÞϕðtÞn½m_ 0f
ðtÞm_ f ðtÞ and go to step 3. This step is different from Ge or geet al. [30] as they use O/F rati of or checking conver-gence. As m_ f ðtÞ is in the denominator of O/F ratio using that couldlead tonon-convergence .Also, Ge or geetal.
[30] assumes an arbitrary value of ϕðtÞ, which could again lead to non-conver genceand oscill atory solution.