You can use Matlab's 'trapz' successively for this purpose. The easiest case is integrating over a rectangular region in (x,y) - that is, where the tabular data, T(x,y), is defined over a rectangular mesh of (x,y) values. Let x and y be vectors of the (x,y) values along the two edges of the rectangular table.
[X,Y] = meshgrid(x,y);
T = (The array of tabular values corresponding to X and Y)
I = trapz(y,trapz(x,T,2));
You should realize that the accuracy of this method is limited by how closely-spaced the x and y values are in the rectangle. For higher order integration than the trapezoidal method there are some algorithms given in the MatLab File Exchange.
Best Answer