Revision history [back]

It does seem that a sort command or option is missing. I can get closer to what you want with:

var('c0 c1 c2 x0 x1 x')
p=(c0+(x-x0)*(c1+(x-x1)*c2))
p.expand().collect(x)

This collects terms better, but they are not in order by degree.