Calculate $\int_CPdx+Qdy$ where $P(x,y)=xe^{-y^2}$ and $Q(x,y)=-x^2ye^{-y^2}+\frac{1}{x^2+y^2+1}$, $C$ is the boundary of the square determined by the inequalities $-a\leq x\leq a$, $-a\leq y\leq a$, oriented positively.
I have thought to do the following:
Using Green's theorem we get to that
$\int_CPdx+Qdy=\int\int_D(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y})dA= \int\int_D(-2xye^{-y^2}-\frac{2x}{(x^2+y^2+1)^2}+2xye^{-y^2})dA=\int\int_D-\frac{2x}{(x^2+y^2+1)^2}dA=\int_{-a}^{a}\int_{-a}^{a}-\frac{2x}{(x^2+y^2+1)^2}dxdy=-\int_{-a}^{a}\int_{a^2+y^2+1}^{a^2+y^2+1}u^{-2}dudy=-\int_{-a}^{a}0dy=0$
Is this fine? Thank you.
Since your region $D$ is symmetric about the $x$-axis, any integrable function $f$ satisfying $f(x,y) = -f(-x,y)$ will satisfy $$\int_D f(x,y) \, dA = 0$$ by symmetry.