\r ellpn.gp
allocatemem();
allocatemem();
allocatemem();
allocatemem();
allocatemem();
fp2coefs([1,x,x,3,-2],x^2-x-10,2,20000);