special relativity code snippets in matlab