Function with three parameters and I want to integrate in one dimension

I have this function:

fn_E(x,y,f) = (1 + 0.25*sin(2*pi((x*f)/(4*c))))*exp(-1i*pi/4*sin(2*pi((x*f)/(4*c))));

f is a 1x71 vector of frequencies 0.140GHz to 0.280GHz. C is the speed of light. I want to integrate the function above at x (or r), where x (or r) is between -12.192/2 and 12.192/2. I want to plug this in for x (or r):

r*sin(theta(j) + Pi/2) %where theta(j) is equal to linspace(0,360,513) in radians with the 513th term removed. So its a 1x512 vector.

Now I can integrate the function fn_E. In the integral, fn_E is also being multiplied by

EX = exp(-1i*2*k(m)*r*cos(theta(j)+Pi/2)) %where r is replacing x at this point and k is the wave number that is a 1x71 vector as well

The output of the integral of from -12.192/2 to 12.192/2 of fn_E times EX needs to be a 512x71 matrix.

So the integral from -12.192/2 to 12.192/2 that needs to occur is

@(r) (1 + 0.25*sin(2*pi((r*f)/(4*c))))*exp(-1i*pi/4*sin(2*pi((r*f)/(4*c))))*exp(-1i*2*k(m)*r*cos(theta(j)+Pi/2))

%where r in @(r) is (r*sin(theta(j) +Pi/2))

THANK YOU.