[ \operatornameInd p(\gamma) ;=; \frac12\pi\int S^1 \frac(x - p_x),dy - (y - p_y),dx(x-p_x)^2 + (y-p_y)^2 ]