An alternative solution to the one I linked in the comments:
Let
z = f(x, y) = √(64 - x² - y²)
Then the integral of y² over the surface S is
We have
∂f/∂x = -x/√(64 - x² - y²)
∂f/∂y = -y/√(64 - x² - y²)
so that the integrand simplifies to
To compute the remaining integral, convert to polar coordinates using
x = r cos(θ)
y = r sin(θ)
dx dy = r dr dθ
The region x² + y² ≤ 16 is given by the set
D = {(r, θ) : 0 ≤ r ≤ 4 and 0 ≤ θ ≤ 2π}
Then the integral becomes
Now,
and by substituting u = 64 - r²,
So, the surface integral has a value of