quadrature points code snippets in matlab