radiation code snippets in matlab