Has anyone worked out a way to allow numerical 2D integration, with vector inputs for the upper/lower bounds in the two dimensions? The mathworks FAQ page describes an m-file called "dblquad", which calls "quad.m" twice in a nested fashion to perform the integration over the two dimensions, but only returns a single value, i.e.:
x = dblquad('function_xy',xlow,xhigh,ylow,yhigh)
where x, xlow, xhigh, ylow, yhigh are scaler values (1x1).
I'm modeling a focal plane array in Matlab as an MxN array of pixels, with a continuous point-spread function. I want to evaluate the total energy over each pixel by integrating the 2D point-spread function over each pixel, like in the following pseudo-code...