Abstract
This paper deals with the numerical modeling for transmission dynamics of Leptospirosis disease in juvenile and adult humans. A novel unconditionally stable Non-Standard Finite Difference (NSFD) numerical scheme is proposed and analyzed to investigate the behavior of the disease in a population. Numerical investigations are provided and results are compared with standard finite difference schemes being already used to handle such problems. Standard finite difference schemes give conditional convergence and do not behave well in certain scenarios. The proposed numerical scheme is dynamically consistent with the biological nature of the continuous model and preserves all of its essential properties.