求教关于 Simpson' s rule的实现

求教各位大神，由于课题需要，我需要实现连续层面图像上一个不规则肿块的体积，怎么通过划取每一层的肿块边界，最后积分出体积啊？Matlab怎么实现呢？我看文献上说用（modified） Simpson’s rule。我是小白，请各位大神赐教。论坛币不多，聊表心意。