"I get [x=- 6.8493749048472302280238361002495b0, x=- 1.5748919531068678397256811243939b-1, x=1.1962965565382104296086949728645b-1] for the 3 roots." You have to set logexpand to false to get three roots BTW. log(x^2) is not equivalent to 2*log(x) for negative x. Regards. Rich