“…In this work, we used dassl as our ode solver. This technique has been used to solve a variety of important problems in the past, including the solutions of adsorption problems involving steep gradients in bidisperse particle [6], hyperbolic models of transport in bidisperse solids [7], transport problems involving steep concentration gradients [8], and modelling saltwater intrusion into coastal aquifers [9,10].…”