Matlab - ayuda ha realizar un problema

Yo quiero ver si alguien puede ayudarme a resolver este problema lo mas rapido que se pueda pues es para un trabajo que estoy haciendo... El caso es el siguiente:

Elaborar un programa que calcule integrales dobles iteradas en un rectángulo de R2, mediante el método de Simpson con una tolerancia de error especificada por el usuario. El programa debe devolver el valor aproximado de la integral y el error de integración que esta aproximación posee.

Gracias