Ralf Stephan on Tue, 29 Jun 2004 17:27:19 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
list / pointers lost |
Reading this into 2.2.8: /*************************** start */ f(pol)= { local(i,g,pr,p); i = 1; g = gcd(pol, subst(pol, x, x+1)); pr = pol / subst(g, x, x-1); p = listcreate(1000); listput (p, gcd (pol/g, pr)); while (g<>1, i = i + 1; pol = g; pr = pr / p[i-1]; g = subst (pol/pr, x, x+1); listput (p, gcd (pol/g, pr)) ); p } f((x-3)*(x-2)^2*(x+1)^4*(x+2)^2*(x+3)*(x+5)^3*(x+6)) /*************************** end */ gives *** significant pointers are lost in gerepile !!! (please report) ralf