Jeroen Demeyer on Wed, 23 Mar 2016 12:04:14 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: Proposal to extend primes()


On 2016-03-16 10:32, Karim Belabas wrote:
Currently a single argument describes the interval (as in e.g. polrealroots),

I guess you mean polsturm(). But that actually does support polsturm(pol, a, b) as well as polsturm(pol, [a,b]). So the following are equivalent:

gp> polsturm(x^3+1,[-oo,1])
%1 = 1
gp> polsturm(x^3+1,,1)
%2 = 1