> < ^ From:

> < ^ Subject:

On Mon, Jul 28, 2003 at 01:32:25PM +0100, Derek Holt wrote: > Dear GAP Forum, > > (I sent this message early in July, but I did not see it appear in the > GAP Forum, so apologies if it appears twice!) > > I have an expression, z, say, involving the n-th roots of unity E(n), which > evaluates to a real number (i.e. ComplexCojugate(z) = z.) For example > z := E(5) + E(5)^-1 + E(3) + E(3)^-1; > > Is there any way to test within GAP whether z is greater than 0, numerically?

I believe there is two answers:

1) You can't, because z is defined up to conjugacy and some conjugates

can be positive and other negative. For example z:=E(8)+E(-8) can be both

sqrt(2) and -sqrt(2).

2) You can use Sturm algorithm with the minimal polynomial of z.

At least that will tell you how many positive conjugates there are.

I don't know if this algorithm is available in GAP, but it is not

difficult to implement.

In your example, Sturm algorithm show there are two negative roots.

Cheers,

Bill.

> < [top]