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? 

Similar questions and discussions