use of prodeulerrat in a script

[Alessandro Languasco <languasc@math.unipd.it>]

Dear all,

I am experimenting with prodeulerrat and I am trying to generate

(randomly ) a rational function g(x)= a(x)/b(x) and then compute

prodeulerrat(g(1/p))

g is such that deg(b)>deg(a)+2, so the product converges.

Now, I’d like to do that in a script; so 

- in input I have a parameter y (an integer),

- I generate first g(x,y) as a polynomial in x

- then I would like to compute prodeulerrat(g(1/p,y))

If I do so in the shell, it works nicely. 

But I would like to write a script that do so and I am stuck on this

because I do not know how to “pass” the information about 

perform the infinite product on p.

What am I missing here ?

Do you have any suggestion ? 

Thanks in advance,

Alessandro