In $[62, \S 2.4]$ the measure $\mathrm{d} \lambda(t)=\exp \left(-t^2-\sin t\right) \mathrm{d} t$ on $\mathbb{R}$ is considered as an "exotic choice" for illustrating the Riemann-Hilbert approach to Gaussian quadrature. The specific example given is $I=\int_{\mathbb{R}} f(t) \mathrm{d} \lambda(t)$ with $f(t)=\sin ^2 10 t$.
(a) Use the SOPQ routines sr_hermite.m and sgauss.m to obtain I accurate to 24 decimal digits.
(b) Use the QPQ routine modis.m to genenate the first $N$ recurrence coefficients $\alpha_k, \beta_k$, $k=0,1,2, \ldots, N-1$, for the measure $\mathrm{d} \lambda$.
(c) Reproduce Fig. 5(b) of [62] showing the absolute error of n-point Gaussian quadrature, $n=1,2, \ldots, 120$, for the integral $I$.