Hi, I found that the code ex e=indexed(a+b*(c+d),idx(i,3)); e=e.expand(expand_options::expand_indexed); cout << e << endl; cout << e.expand(expand_options::expand_indexed) << endl; gives the output a.i+(d*b+c*b).i (c*b).i+a.i+(d*b).i I would have expected to see the second result twice instead. A patch is attached. Bye, Chris