By Green's theorem,
[tex]\displaystyle\int_C\cos y\,\mathrm dx+x^2\sin y\,\mathrm dy=\iint_D\left(\frac{\partial(x^2\sin y)}{\partial x}-\frac{\partial(\cos y)}{\partial y}\right)\,\mathrm dx\,\mathrm dy[/tex]
where [tex]D[/tex] is the region with boundary [tex]C[/tex], so we have
[tex]\displaystyle\iint_D(2x+1)\sin y\,\mathrm dx\,\mathrm dy=\int_0^5\int_0^4(2x+1)\sin y\,\mathrm dy\,\mathrm dx=\boxed{60\sin^22}[/tex]