Use the transformation x= v/u, y= uv to find ∫∫ᵣ xy³ dA over the region R in the first quadrant enclosed by y = x. y = 3x, xy = 1, xy = 4.