Yes, nasty result but it seems correct :) I don't exactly understand how the procedure t -> has(evalf(t),I) works to get a boolean-value in your remove-function. How does it know which values to remove? Many thanks!