integrate((1/2)*ln(y^2/(y^2+pi^2))*ln(1-e^(-2*y))dy from 0 to 10)