code: def build_delta_t_table(delta_t_recent):
      ancient = load_bundled_npy('morrison_stephenson_deltat.npy')
      historic = load_bundled_npy('historic_deltat.npy')
      historic_start_time = historic[0,0]
      i = searchsorted(ancient[0], historic_start_time)
      bundled = concatenate([ancient[:,:i], historic], axis=1)
      recent_start_time = delta_t_recent[0,0]
      i = searchsorted(bundled[0], recent_start_time)
      row = ((0,),(0,))
      table = concatenate([row, bundled[:,:i], delta_t_recent, row], axis=1)
      century = 36524.0
      start = table[0,1] - century
      table[:,0] = start, delta_t_formula_morrison_and_stephenson_2004(start)
      end = table[0,-2] + century
      table[:,-1] = end, delta_t_formula_morrison_and_stephenson_2004(end)
      return table