point spread function code snippets in matlab