Calculate surface integral: $$I = \iint \limits_{S} \frac{xy}{x^2 + y^2} d{S}$$ where S is surface determined by sides of pyramid $x + y + z = 1 \phantom\ (x, y, z \ge 0 ),$ inside sphere $(x - \frac{1}{2})^2 + (y - \frac{1}{3})^2 + z^2 \le \frac{1}{4}.$
How should I parameterize the surface $S$? I am also having a question regarding the case of the flux integral: How should I find the unit normal vector of the surface $S$?