First, recall that Gaussian quadrature is based around integrating a function over the interval [-1,1], so transform the function argument accordingly to change the integral over [1,5] to an equivalent one over [-1,1].



So,

Let

. With

, we're looking for coefficients

and nodes

, with

, such that

You can either try solving for each with the help of a calculator, or look up the values of the weights and nodes (they're extensively tabulated, and I'll include a link to one such reference).
Using the quadrature, we then have
