Write and call a function in Matlab to implement Gauss quadrature. At a minimum, the function

should require as input: $x,$ function to be integrated. At a minimum, the function should output:

value of the integral. The program should be capable of $2-4$ point Gauss quadrature. Test it by

integrating

${\displaystyle \underset{0}{\overset{1}{}}}{x}^{0\mathrm{.}1}(1\mathrm{.}2-x)(1-e^{20(x-1)})dx$

You should output the calculated value of the integral and the absolute relative true error ${}_t$ for

the $2,3$ and $4$ point Gauss quadrature $($use the true value of $0\mathrm{.}602298).$