integrate((ax+r)^2/(1+e^x))