用对面积的曲面积分喽
假设曲面的方程是x^2+y^2+z^2=R^2,由对称性,只考虑第一卦限部分的面积
第一卦限的球面的方程是z=√(R^2-x^2-y^2),αz/αx=-x/z,αz/αy=-y/z
dS=Rdxdy/√(R^2-x^2-y^2)
第一卦限的球面在xoy面的投影区域是D:x^2+y^2≤R^2,x≥0,y≥0
所以,球面面积S=8∫∫Rdxdy/√(R^2-x^2-y^2)=8∫(0~π/2)dθ∫(0~R) R/√(R^2-ρ^2)ρdρ=4πR∫(0~R) 1/√(R^2-ρ^2)ρdρ=4πR^2