torus code snippets in matlab