cubic code snippets in matlab