Change behaviour of integrate_times to prefer observer y.

The only reason to call this (over integrate_adaptive) is to get the
matrix of y values out, so this should just be produced as the main
output.