The integrate_1d function in the Stan language accidentally requires that you specify the relative tolerance as the last argument, which is contrary to the documentation. The default value was supposed to be sqrt(epsilon()).
2 Likes