• This paper analyzes a seven-dimensional system of ODEs as model malaria propagation and seeks qualitative properties the on numerical levels so-called dynamical consistency. The is solved by nonstandard θ -method including nonlocal scheme standard to right side step-size function approximate time derivatives with first-order consistency valid for any step size. addresses reliable mathematical...