heun's method code snippets in matlab