Zeta = lfuncreate(1)
lfunan(L,20)
lfun(Zeta,2)
lfun(Zeta,0,1)
lfun(Zeta,1) 
lfun(Zeta,1+x+O(x^10))
lfunzeros(Zeta,20)
lfunlambda(Zeta,2)

G=znstar(4,1); G.clgp
Dir=lfuncreate([G,[1]]); Dir[2..5]
lfunan(Dir,30)
lfun(Dir,2)
Catalan
znconreyexp(G,[1])
lfun(Mod(3,4),2)

Dedek = lfuncreate(x^2+1); Dedek[2..5]
lfun(Dedek,2)
zeta(2)*Catalan
L=lfunmul(Zeta,Mod(3,4));
lfun(L,2)
L2=lfundiv(Dedek,1);
lfun(L2,2)



E = ellinit([0,-1,1,-10,-20]);
L=lfuntwist(E,Mod(2,5));
lfunan(E,10)
lfunan(Mod(2,5),10)
lfunan(L,10)

nf=nfinit(polcyclo(5,’a));
E2=ellinit(E[1..5],nf);
localbitprec(64); lfun(E2,2)
L2=lfuntwist(E,Mod(4,5));
lfun(E,2)*lfun(L2,2)*norm(lfun(L,2))

L=lfungenus2([x^2+x,x^3+1]);
L[2..5]
lfun(L,1)
lfunan(L,5)

bnf = bnfinit(a^2+23);
bnr = bnrinit(bnf, 1);
bnr.clgp
Hecke = lfuncreate([bnr,[1]]);
Hecke[2..5]
z=lfun(Hecke,0,1)
algdep(exp(z),3)

N = nfinit(x^6+108);
G = galoisinit(N);
[T,o] = galoischartable(G);
T~
galoisconjclasses(G)

L = lfunartin(N,G,T[,3],o);
lfuncheckfeq(L)
L[2..5]
z = lfun(L,0,1)
algdep(exp(z),3)
bnr = bnrinit(bnfinit(a^2+a+1),6);
lfunan([bnr,[1]],100)==lfunan(L,100)