polynomial-solving code snippets in matlab