aberration code snippets in matlab