In the same way, let:
F(x,y) = [-x²y,xy²]
Then:
div(F(x,y)) = x²+y²
The integral (moment of inertia) results:
I = int( dot( F(g(t)) , perp(g'(t) ) ) dt =
= int( -gx(t)²gy(t)d(gy(t))/dt + gx(t)gy(t)²d(gx(t))dt ) dt
Where:
g(t) = [gx(t),gy(t)]
Did I do any mistake ?