Suppose I'm getting a set of three solutions from a pde simulation(assume that each solution is a matrix of dimension 100x500). I run this until I get the steady state. Then at the final time where I have the solution, I need to numerically integrate one pde over the space and divide by the length of the domain(simply get the average of the solution). How can I do this with MATLAB?