velocity addition code snippets in matlab