Double Integral over function that returns a vector

I am trying to calculate a double integral using the integral2 command. The function I am trying to integrate over takes in two scalars, x and y, and returns a vector of size 5000. Currently, I am using a for loop which calculates the double integral for each element of the function output. However, I can't figure out how to skip the for loop and calculate the integral for each of the 5000 values all at the same time.