I am currently working on a problem involving Young´s double slit interference experiment. What I want to do is to evaluate a integral over a surface that looks something like this:

where the integral is to be taken over the red surface only, and the integrand depends on the radius of the circle. What would be the best way to do this numerically in matlab?