A simple test on Dahlquist:.
model = otp.linear.presets.Dahlquist(-1);
sol = model.solve('Reltol',1e-10, 'Abstol', 1e-10, 'Method', @ode15s);
sol.y(end)-exp(-1)
4.162909795304426e-04
if we switch to ode45
sol = model.solve('Reltol',1e-10, 'Abstol', 1e-10, 'Method', @ode45);
sol.y(end)-exp(-1)
1.209031486038015e-09