If we take our pool equation and solve for y we get:

This is the length of the base of the square as the cross sections stack up along the x-axis.

But the cross-sections are squares:

We squared our function and multiplied by 4 because the integral gives the area in one quadrant of the ellipse.

What we end up with is a very easy integral once it is set up.