simson's rule code snippets in matlab