finite difference method code snippets in matlab