gap> Q:=QuadraticNumberField(-1);;
gap> Lfunction(Q,2)/(18*3.142);
0.0161957

gap> 1.0*Log(Product(AbelianInvariants(F)),10)/Norm(I);
0.00913432
