gamma function code snippets in matlab