I’ve also used discretized solvers in my StanCon 2018 submission. In my case, I’ve found the accuracy of Euler to be OK, but systematically biased compared to RK45. However, the Trapezoidal rule achieved (in my case) results literally indistinguishable from RK45 for just a tiny increase in computation cost and code complexity (compared to Euler).