A Nonstandard Dynamically Consistent Numerical Scheme Applied to Obesity Dynamics

Figure 1

Profiles of subpopulations 𝑁(𝑡),𝐿(𝑡),𝑆(𝑡),𝑂(𝑡),𝐷_𝑆(𝑡), and 𝐷_𝑂(𝑡). Solutions are obtained by the proposed NSFD scheme with 𝜓(ℎ)=ℎ, where ℎ=0.2 and initial conditions are 𝑁(0)=0.462, 𝐿(0)=0.194, 𝑆(0)=0.2176, 𝑂(0)=0.09, 𝐷_𝑆(0)=0.0249, and 𝐷_𝑂(0)=0.0115.

(a)

(b)

(c)

(d)

(e)

(f)