expr, logexpand:all factor(%) substpart(logcontract(piece),%,2) scanmap(factor,%) On Thu, May 31, 2012 at 10:49 AM, David Ronis <David.Ronis at mcgill.ca> wrote: > log(((k[1]^2+2*k[1]*k[0]+k[0]^2)*D+z > +Gamma)^2/(((k[1]^2-2*k[1]*k[0]+k[0]^2)*D+z+Gamma)^2)); >