\r ellpn.gp
allocatemem();
allocatemem();
allocatemem();
allocatemem();
allocatemem();
fp2coefs([0,2,1,-19-8*x,28+11*x],x^2-x-9,2,20000);