'Efficiently compute Taylor shift ``f(x + a)``.'

def shift(f, a):

    if (not f.lev):
        return f.per(dup_shift(f.rep, f.dom.convert(a), f.dom))
    else:
        raise ValueError('univariate    polynomial    expected')
