1

This is probably a pretty stupid question--but I can't seem to find a relevant explanation in the pari-gp literature. I'm looking for a function in the library that takes a polynomial, or a series, or a constant, and spits out the constant portion (in the constant case, it'll just spit out the input).

The tool I've been using is polcoef(Expression,0) but this doesn't work for multiple variables--we have to call the function again.

For example,

polcoef(1+z+z^2,0)
%2 = 1

polcoef(1+y+z,0)
%3 = z + 1

polcoef(polcoef(1+y+z,0),0)
%4 = 1

Is there a built in function that'll do this for an arbitrary number of variables? So that, I don't have to iterate the polcoef function. This would clean up my code considerably; and I know I could probably build a recursive function myself, but I'm wondering if there's something built-in that does this.

Any help, comments, requests to edit and clarify are greatly appreciated.

Regards.

Edit:

Forgot to add, my desired output is,

Const(1+x+y+z)
1

Const(Pi+x^2+p)
Pi
Richard Diagram
  • 209
  • 1
  • 7

1 Answers1

5

You can use the built-in substvec function to replace several variables at once. For example:

substvec(1+y+z, ['y,'z], [0,0])
%1 = 1

You even can easily implement a short-cut function not to create the list of zeros of proper length:

Const(poly: t_POL) = {
    my(vars = variables(poly));
    substvec(poly, vars, vector(#vars))
};
Const(1+x+y+z)
%1 = 1

Const(Pi+x^2+p)
%2 = Pi
Piotr Semenov
  • 1,761
  • 16
  • 24
  • 1
    Thank you so much! That looks much better than the method I was thinking which was to iterate ```polcoef(Expression,0)``` until its constant, lol. – Richard Diagram Aug 23 '21 at 18:18