On 04/02/13 10:15, Stavros Macrakis wrote: I changed to use imagpart(rhs(eq3[i])) = 0 Thanks. Tom Dean