I think this sounds right.  I'll leave the current GK integrator in place for now as an option, but I think I'll make the default integrator the trapezoid rule evaluated on `bandpass.wave_list` (have to write this still), or evaluated evenly between `blue_limit` and `red_limit` with N=100 steps if `bandpass.wave_list` doesn't exist.