Package: pari

Version: 2, 6, 0

 

Here are some inputs and the resulting output:

 

(06:29) gp > nsb(n) = numerator(subst(bernpol(n), x, 3));
(06:30) gp > anb(nn) = {for(n=0, nn, print1(nsb(n), ", "););}
(06:30) gp > anb(20)
1,
  ***   at top-level: anb(20)
  ***                 ^-------
  ***   in function anb: for(n=0,nn,print1(nsb(n),", ");
  ***                               ^--------------------
  *** print1: bug in PARI/GP (Segmentation Fault), please report.
(06:30) gp >

But one by one, no problem:

(06:33) gp > nsb(0)
%4 = 1
(06:33) gp > nsb(1)
%5 = 5
(06:33) gp > nsb(2)
%6 = 37

 

Do you see a mistake I did ?

 

Thanks.

 

Michel