@Kitonum 

Thanks a lot. It works well, however, when I use Greek letters, for example, varepsilon instead of x, Maple does not display the symbol for varepsilon, it displays varepsilon instead of the symbol. Also if I use subscripts, like x__i, it displays x__i in the output.

Do you know how to fix this issue? Thanks

@Kitonum Thanks. How can I get this to work if the coefficients of the polynomial are variables (with assumptions on them) rather than numeric values?

@Carl Love Thanks. I was thinking there should be at least a warning which tells that the assumptions are ignored because I sometimes see those kinds of warning. Do you have any idea why there is no warning?