Sumas dobles de riemann

Realizar un programa en matlab que calcule el volumen de un sólido mediante una función arbitraria(dada por el usuario) usando el punto medio de los subintervalos. Un ejemplo en realidad es el uso de la regla del punto medio, solo que no se debe de realizar con integrales si no solo con la suma de riemann, preferentemente si es con ciclo for, se le debe de pedir al usuario los limites tanto inferiores como superiores de x(a,b) y y (c,d), al igual que n y m; n representando en cuántas partes estará dividida X y m representando en cuántas partes estará dividida Y. Por último, calcular deltaX con (b-a)/n y deltaY con (d-c)/m para así poder calcular con mayor facilidad la suma