integrate (sin(x-phi)-sin(x))^2 over x from 0 to 2pi