integrate w*(arcsin(min(1,w)/w)-arcsin(sqrt(max(0,w^2-1))/w)) dw from 0 to sqrt(2)