jjWj.sin.j2

# arXiv:0708.4233v2
#
# Section 3.1.2
#
# Table 1, Row 2

include("{{ SIN_LOCAL }}/ggW.sin")

process j1 = j, j => e1, N1, j

# calls(adapt) = 5e7, calls(integ) = 5e7
integrate (j1) {iterations = 10:500000:"gw", 5:1000000}

printf "RES %8e %8e" (integral(j1), error(j1))