解(1)、曲线绕$y$轴旋转后的曲面方程为:\[x^2+3y^2+z^2=1,(x^2=x^2+z^2),\]过$P(x,y,z)$的法向量为
$\mathfrak{n}=\{x,3y,z\}.$而\[\rho (x,y,z)=\frac{1}{\sqrt{x^2+(3y)^2+z^2}}.\]于是\[\begin{align*}I
&=\iint\limits_{S}\frac{z}{\rho (x,y,z)} dS\\
&=\iint\limits_{S}z\sqrt{x^2+(3y)^2+z^2} dS \\
&=\iint\limits_{S}\frac{z(x^2+(3y)^2+z^2)}{\sqrt{x^2+(3y)^2+z^2}} dS \\
&=\iint\limits_{S}[\frac{zx\cdotp x}{\sqrt{x^2+(3y)^2+z^2}}+\frac{3zy\cdotp 3y}{\sqrt{x^2+(3y)^2+z^2}}+\frac{zz\cdotp z}{\sqrt{x^2+(3y)^2+z^2}} ]dS \\
&=\iint\limits_{S}[zx\cos \alpha +3zy\cos \beta +z^2\cos \gamma ]dS \\
&=\iiint\limits_{\Omega }(z+3z+2z)dV,(高斯公式) \\
&=6\iiint\limits_{\Omega }zdV \\
&=6\iint\limits_{\Omega }dxdy\int_{0}^{\sqrt{1-x^2-3y^2}}zdz \\
&=3\iint\limits_{\Omega }(1-x^2-3y^2)dxdy \\
&=3\int_{0}^{1}r(1-r^2)dr\int_{0}^{2\pi } d\theta \\
&=6\pi (\frac{1}{2}r^2-\frac{1}{4}r^4)|_0^1\\
&=\frac{3}{2}\pi .
\end{align*}\]
(2)、与(1)结果相同。