Write a MATLAB function program to do the center-point method for integration of a function $f(x, y)$ on a rectangle $a \le x \le b$, $c \le y \le d$ using $m$ subintervals in $x$ and $n$ subintervals in $y$. Let the inputs be $(f, a, b, c, d, m, n)$. Include comments.