restart; macro(e = epsilon); Order := 2; alias(alpha = RootOf(z^5-1, z)); f := u -> u^5 - e*(u+1); z := convert( series( RootOf(f(u),u), e), polynom); Delta := series(f(z), e, Order+2): map(simplify, Delta);