gauss-legendre code snippets in matlab