Vectorised generation of multiple time-series

Allow for generation of arrays of time series thanks to contributions from
Alex Spaeth. The second positional argument "samples" was renamed to "size" to
reflect this. This change is backwards compatible unless "samples" was used as
a keyword argument.