imptaylor:= (eq, y, b, x, a, n) -> taylor(b +(x-a) * convert(taylor(RootOf(subs(y=b+(x-a)*_Z, eq)),x=a,n-1), polynom) + (x-a)^n, x=a, n);