midpoint method code snippets in matlab