integrate x*y/(1-x*y)*(1/log(x)+1/(1-x))*(1/log(y)+1/(1-y)-1)+eulergamma^2-eulergamma+Pi^2/24,x=0,1,y=0,1