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

gap> Loge10:=0.434294481903;; #Log_10(e)
gap> T:=Product([ 2, 2, 2, 2, 3, 4, 5, 16, 31, 223, 857, 13913, 106103, 
>                 1267133821, 141880858099, 1290077458421, 92119226692240157, 
>                 124336989193985597 ]);;
gap> 1.0*Log(Product(T,10)/(Loge*Norm(I));
0.0147988
