You can get symmetric interpolants as a special case of tree interpolants. Here's an example:
(declare-const a Int)
(declare-const b Int)
(declare-const c Int)
(declare-const d Int)
(compute-interpolant
(and
(interp (<= 0 a))
(interp (and (= a b) (= b c)))
(interp (and (= c d) (<= d -1)))))
Here is z3's output:
unsat
(>= a 0)
(= a c)
(<= c (- 1))
Have a look at the tutorial at http://rise4fun.com/iZ3/tutorial/guide under "Tree interpolants".