A standard way to get E(X) on a continuous random variable is to take the integral over the full range of the product of x and the function. The same would go for E(Y). In this case, it would be a double integral. Both x and y run from 0 to 1, so both will be evaluated for the extreme values 0 and 1. I can't put those limits on the integral symbols on this editor, so I'm explaining it here.
∫∫ x · 4xy dy dx = ∫∫ 4x2y dy dx
We do the integral with respect to y first, working from the inside out, and we get
∫ 4x2y2/2 [from 0 to 1] dx
When we evaluate the interval for y = 1 and y = 0 and subtract, we get
∫ 2x2 dx [from 0 to 1]
This integral is 2x3/3, evaluated from 0 to 1, which is 2/3.
I could do the one for E(Y), which uses ∫∫ 4xy2 dy dx, but you may begin to notice that we're doing exactly the same operations with y as with x, and using the same numbers. So guess what E(Y) will be? That's right! 2/3 also.
Now the one for P(Y ≥ 1 − X) will involve an integral of the given pdf, but the part that is above the y = 1 - x line. But perhaps an easier way is to evaluate the integral below the 1 - x line and then subtract from 1.
So here's roughly what the double integral will look like:
∫ [0 to 1] ∫ [0 to (1 - x)] 4xy dy dx
∫ 4xy2/2 [0 to (1 - x) dx
∫ 2x (1 - x)2 dx
∫ (2x3 - 4x2 + 2x) dx [0 to 1]
(2x4/4 - 4x3/3 + x2) [0 to 1] = 1/6.
This is P(Y < 1 - X), so P(Y > 1 - X) is 1 - 1/6 = 5/6.