Add variadic version of integrate_1d

[WardBrian]

The integrate_1d function already supports variadic arguments in it's C++ implementation, but not in the function exposed in the language.

Similar to #1268

We will likely want to give it a new name, possibly with the algorithm mentioned to support https://github.com/stan-dev/math/issues/3000. Maybe integrate_1d_double_exp?