Some new GP features
lfuneuler
lfuneuler(L,p) returns the Euler factor at p of the L-function
L assuming it is known. This is especially useful at bad places.
? lfuneuler(Mod(2,5),3)
%47 = 1/(I*x+1)
? lfuneuler(Mod(2,5),5)
%48 = 1
? L=lfungenus2([x^2 + x,x^3 + x^2 + 1]);
? lfuneuler(L,11)
%50 = 1/(121*x^4+11*x^2+1)
? lfuneuler(L,13)
%51 = 1/(13*x^2+5*x+1)