number theory code snippets in matlab