“…So far in the literature, UA methods have been constructed by different techniques by several authors [1,3,4,6,9,10,11,12,13,15,16,18]. It is out of the scope of this paper to describe all of them and for an extensive comparison, we refer the interested reader to [25]. However, all these methods have in common that they necessitate not only an analytic expression of f θ in (1.1), but also higher-order derivatives ∂ p u f θ in explicit form: instead of solving directly (1.1) the system is first transformed into a nonstiff one which is more amenable to numerical integration.…”