Fit a straight line y = a + b t to a series (with t = 0, 1, ..., n-1) by least squares and forecast h steps ahead.
Implement trend_forecast(y, h) returning (a, b, forecast) where forecast is an array of the next h values (t = n, ..., n + h - 1).