quadrature weights code snippets in matlab