Dear friends,

I ran into another problem while using Maple to do residue calculus. The following call

residue(1/x^(1/3)/(x^2+2*x*cos(phi)+1), x=-exp(I*phi));

returns zero -- it does not recognize the residue. On the other hand, if I do a subsitution like this

subs(x=-exp(I*phi), 1/x^(1/3)/(x^2+2*x*cos(phi)+1));

followed by

convert(%, trig); simplify(%);

then I get a divide by zero error, which shows that the pole is there. Even when we leave out the power of

x,

which might cause problems owing to the fact that we have a branch cut connecting zero to infinity, the residue is still not found. The call

residue(1/(x^2+2*x*cos(phi)+1), x=-exp(I*phi));

is zero again but the substitution

subs(x=-exp(I*phi), 1/(x^2+2*x*cos(phi)+1));

followed by

convert(%, trig); simplify(%);

produces the divide by zero. Curiously enough, the call to

singular(1/(x^2+2*x*cos(phi)+1),x);

does find the right singularities, but does not express them in their simplest form. If this is what's used internally that could explain why it does not recognize that the two expressions of the pole represent the same complex number.

Marko Riedel

Read more about it: here