matrices code snippets in r