Calculate $\iint_S \vec F \cdot \vec n dS$.
While $S$ is between $z=x^2+y^2$, $z=\sqrt{x^2+y^2}$.
$\vec F=(3x,4y,-z)$. (The normal to the outside of S)
Using gauss law: $$\iint_S \vec F \cdot \vec n dS=\iint_D(\int^{\sqrt{x^2+y^2}}_{x^2+y^2}6)=6\int^{2\pi}_{0}(\int^1_0(r-r^2)r dr)d\theta)=12\pi\frac{1}{12}=\pi.$$
Normal way: On the paraboloid "$S_1$" : $\vec n=(2x,2y,-1)$.
$$\iint_D(\vec F)\cdot(\vec n)dxdy = \iint_D6x^2+8y^2+(x^2+y^2)dxdy=\int_0^{2\pi}[\int_0^1(7r^2+2r^2sin^2(\theta)dr]d\theta=2\pi\frac{7}{4}=\frac{7\pi}{2}.$$
On the cone $S_2$ : $\vec n=(\frac{-x}{\sqrt{x^2+y^2}},\frac{-y}{\sqrt{x^2+y^2}},1)$.
$$\iint_D(\vec F)\cdot(\vec n)dxdy=\iint_D-[\frac{3x^2+4y^2}{\sqrt{x^2+y^2}}+\sqrt{x^2+y^2}]=\int_0^{2\pi}-[\int_0^1(\frac{3r^2+r^2sin^2\theta}{r}+r)rdr]d\theta=-\int_0^{2\pi}[1+\frac{1}{3}+sin^2\theta] d\theta =- \frac{8\pi}{3}.$$
Now adding them together: $\frac{5\pi}{6}$. and according to my first attempt it needs to be $\pi$.
Where's the hidden mistake? (I'm more sure of gauss law attempt since it was easier, but I just want to make sure I don't have any hidden bad understanding of solving integrals).
Any feedback is really appreciated, thanks in advance!