square-root code snippets in matlab