def METHOD_NAME(ts, coors, region=None):
    mtx = nm.ones((coors.shape[0], 1, dim), dtype=nm.float64)
    rhs = nm.arange(coors.shape[0], dtype=nm.float64)[:, None]
    rhs *= 0.1 / (coors.shape[0] - 1)
    return mtx, rhs