gaussian quadrature code snippets in matlab