Sorry, this appears to be a bug. A workaround: substpart(ratsimp(piece),expr,2,1,2); Check ? substpart ? piece for more info. -s On Sun, Jan 24, 2010 at 9:12 AM, gelu precup <geluprecup at yahoo.com> wrote: > expr:x^g(x^(b+1/b)); >