We want to calculate the area of the shape $S$ bounded by the curve $\gamma$ defined by $\begin{cases}x = a\cos^3(\theta) \\ y = a\sin^3(\theta)\end{cases}$, where $a \in \mathbb R$ and $0 \leq \theta \leq 2\pi$.
I solved it. But I'm looking for a different approach, seeing as mine leads to a not very nice integral.
What I did:
$\gamma$ is a closed simple curve. Define the function $F(x,y) = (P(x,y), Q(x,y)) = (0, x)$. It's continuous everywhere with continuous derivatives and everything's nice,so we can use Green's theorem:
$\displaystyle \oint_{\gamma}\vec{F}d\vec{r} = \iint_S \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}dxdy = \iint_S1dxdy = Area(S)$
So we just need to calculate $\oint_{\gamma}\vec{F}d\vec{r}$ which by definition is $\int_{0}^{2\pi}(0,a\cos^3(\theta)) \cdot(-3a\cos^2(\theta)\sin(\theta), 3a\sin^2(\theta)\cos(\theta))d\theta = 3a^2\int_{0}^{2\pi}\sin^2(\theta)\cos^4(\theta)d\theta$.
This integral is doable, but hardly pleasant.
Is there a nicer way I'm not seeing?
Generally, when trigonometric functions are involved, you should try using $$A = \frac12\oint_C -y\,dx+x\,dy$$ rather than $\oint_C x\,dy$, as you did. This tends to make the calculations (far) simpler, as usually you can take advantage of the standard identities.
In this case, you will end up with $\displaystyle\int_0^{2\pi} 3\sin^2\theta\cos^2\theta\,d\theta$, which is easy with the double angle formula and the standard antiderivative $\int\sin^2u\,du$.