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);
