<!DOCTYPE html>
<html lang="en">
<head>
<script>
  MathJax = { 
    tex: {
		    inlineMath: [['\\(','\\)']]
	} }
</script>
<script type="text/javascript" src="https://cdn.jsdelivr.net/npm/mathjax@3/es5/tex-chtml.js">
</script>
<meta name="generator" content="plasTeX" />
<meta charset="utf-8" />
<meta name="viewport" content="width=device-width, initial-scale=1" />
<title>Weighted quadrature formulas <br />for semi-infinite range integrals\(^\ast \): Weighted quadrature formulas <br />for semi-infinite range integrals\(^\ast \)</title>
<link rel="stylesheet" href="/var/www/clients/client1/web1/web/files/jnaat-files/journals/1/articles/styles/theme-white.css" />
</head>

<body>

<div class="wrapper">

<div class="content">
<div class="content-wrapper">


<div class="main-text">


<div class="titlepage">
<h1>Weighted quadrature formulas <br />for semi-infinite range integrals\(^\ast \)</h1>
<p class="authors">
<span class="author">Gradimir V. Milovanović\(^\S \)</span>
</p>
<p class="date">November 6, 2015.</p>
</div>
<p>\(^\ast \) This paper was supported by the Serbian Ministry of Education, Science and Technological Development (No. #OI 174015). </p>
<p>\(^\S \) Serbian Academy of Sciences and Arts, Beograd, Serbia &amp; State University of Novi Pazar, Serbia, e-mail: <span class="tt">gvm@mi.sanu.ac.rs</span>. </p>
<p>Dedicated to prof. Ion Păvăloiu on the occasion of his 75th anniversary </p>

<div class="abstract"><p> Weighted quadrature formulas on the half line \((a,+\infty )\), \(a{\gt}0\), for non-exponentially decreasing integrands are developed. Such \(n\)-point quadrature rules are exact for all functions of the form \(x\mapsto x^{-2}P(x^{-1})\), where \(P\) is an arbitrary algebraic polynomial of degree at most \(2n-1\). In particular, quadrature formulas with respect to the weight function \(x\mapsto w(x)=x^\beta \log ^m x\) (\(0\le \beta {\lt}1\), \(m\in \mathbb {N}_0\)) are considered and several numerical examples are included. </p>
<p><b class="bf">MSC.</b> 65D30, 65D32 </p>
<p><b class="bf">Keywords.</b> Gaussian quadrature rules, nodes, Christoffel numbers, non-exponentially decreasing integrands. </p>
</div>
<h1 id="a0000000002">1 Introduction</h1>
<p>In this paper we consider weighted quadrature formulae on the half line \((a,+\infty )\), </p>
<div class="equation" id="GQCBS">
<p>
  <div class="equation_content">
    \begin{equation} \label{GQCBS} \int _a^{+\infty }w(x) f(x){\, \mathrm{d}}x=\sum _{k=1}^n A_k f(x_k)+R_n(f), \end{equation}
  </div>
  <span class="equation_label">1.1</span>
</p>
</div>
<p> where \(a\) is a finite real number and \(x\mapsto w(x)\) is a given weight function. Such a quadrature formula for \(a=0\) and \(w(x)=x^\alpha {\mathrm{e}}^{-x}\), \(\alpha {\gt}-1\), is the well known <i class="it">generalized Gauss-Laguerre</i> quadrature rule (cf. <span class="cite">
	[
	<a href="#gm-gvm-2008" >10</a>
	, 
	p.
	
	325
	]
</span>), which is exact for all algebraic polynomials of degree at most \(2n-1\), i.e., when \(f\in {\EuScript P}_{2n-1}\). </p>
<p>Error analysis and convergence of such Gaussian formulas on unbounded intervals (with the classical measures of Laguerre and Hermite) was given in 1928 by Uspensky <span class="cite">
	[
	<a href="#Usp" >17</a>
	]
</span>. Otherwise, the corresponding problems for quadrature rules on finite intervals was studied much earlier by <span class="cite">
	[
	<a href="#Posse" >15</a>
	]
</span>, Markov <span class="cite">
	[
	<a href="#Markov1885" >8</a>
	]
</span>, Stieltjes <span class="cite">
	[
	<a href="#Stieltjes1884" >16</a>
	]
</span>, etc. On some new results in this directions see books <span class="cite">
	[
	<a href="#GauBook04" >5</a>
	]
</span> and <span class="cite">
	[
	<a href="#gm-gvm-2008" >10</a>
	]
</span>, including the so-called truncated quadrature rules obtained by ignoring the last part of its nodes (see Mastroianni and Monegato <span class="cite">
	[
	<a href="#gm-gm-2003" >9</a>
	]
</span>). </p>
<p>Very recently Gautschi <span class="cite">
	[
	<a href="#Gautschi2" >6</a>
	]
</span> has constructed a special logarithmically weighted quadrature formula on \((0,+\infty )\), when \(x\mapsto (x-1-\log x)\mathrm{e}^{-x}\). Also, Xu and Milovanović <span class="cite">
	[
	<a href="#Xu_GVM" >18</a>
	]
</span> have developed generalized Gaussian quadrature rules of the form \((\ref{GQCBS})\), with \(x\mapsto w(x)=\mathrm{e}^{-x}\) on \((0,+\infty )\), which are exact on the set of basis functions \(\{ 1,\log x,x,x\log x,\ldots ,x^{n-1},x^{n-1}\log x\} \). In the other words, these rules are exact for each \(f(x)=p(x)+q(x)\log x\), where \(p,q\in {\EuScript P}_{n-1}\), so that they can calculate integrals with a sufficient accuracy, regardless of whether their integrands contain a logarithmic singularity, or they do not. For a similar approach for integrals on the finite intervals see <span class="cite">
	[
	<a href="#FILOMAT15" >12</a>
	]
</span> and <span class="cite">
	[
	<a href="#JCAM15" >14</a>
	]
</span>. </p>
<p>On the other side, a large number of integrals of the form \(\int _a^{+\infty }F(x){\, \mathrm{d}}x\) which appear in applications do not have exponentially decreasing integrands \(F(x)\), and in such cases Gauss-Laguerre quadrature rules are notoriously poor (see Evans <span class="cite">
	[
	<a href="#Evans2005" >2</a>
	]
</span>). As a starting simple example, Evans <span class="cite">
	[
	<a href="#Evans2005" >2</a>
	]
</span> has considered \(F(x)=1/(x^2+0.25)\), where the convergence of the corresponding integral depends on the \(1/x^2\) term for large \(x\). He has proposed a quadrature method based on the set of basis functions \(\{ 1/x^k\} \) and demonstrated its effectiveness on a series of numerical examples. </p>
<p>In this paper we develop a general approach for constructing a class of \(n\)-point generalized quadrature rules \((\ref{GQCBS})\) of Gaussian type on \((a,+\infty )\), \(a{\gt}0\), which are exact for all functions of the form \(x\mapsto x^{-2}P(x^{-1})\), where \(P\) is an arbitrary algebraic polynomial of degree at most \(2n-1\). In particular, we consider quadrature formulas with respect to the weight function \(x\mapsto w(x)=x^\beta \log ^m x\) (\(0\le \beta {\lt}1\), \(m\in \mathbb {N}_0\)), which reduces to the constant weight for \(\beta =0\) and \(m=0\). In order to show the efficiency of the obtained quadrature rules we present a few numerical examples. </p>
<h1 id="GWGR">2 Generalized Weighted Gaussian Rules</h1>

<p>Suppose \(a{\gt}0\), as well as that the weight function \(x\mapsto w(x)\) on \((a,+\infty )\) is such that </p>
<div class="equation" id="uslovi">
<p>
  <div class="equation_content">
    \begin{equation} \label{uslovi} 0<\int _a^{+\infty }\frac{w(x)}{x^2}{\, \mathrm{d}}x<+\infty . \end{equation}
  </div>
  <span class="equation_label">2.1</span>
</p>
</div>
<p> Following Evans <span class="cite">
	[
	<a href="#Evans2005" >2</a>
	]
</span>, we develop a general approach for constructing generalized Gaussian quadrature formulas of the form \((\ref{GQCBS})\). In the cases of integrals on \((\alpha , +\infty )\), when \(\alpha {\lt} a\), we simply take </p>
<div class="displaymath" id="a0000000003">
  \[ \int _\alpha ^{+\infty }w(x)f(x){\, \mathrm{d}}x=\int _\alpha ^{a}w(x)f(x){\, \mathrm{d}}x+ \int _a^{+\infty }w(x)f(x){\, \mathrm{d}}x \]
</div>
<p> and apply to first integral on the right hand side some of rules for calculating integrals on the finite intervals. Also, we mention here that a faster convergence of the corresponding quadrature process can be achieved by taking a greater value of \(a\). </p>
<p>Thus, the basic idea is to construct a quadrature formula of the form \((\ref{GQCBS})\), which is exact for all functions of the form </p>
<div class="displaymath" id="a0000000004">
  \[ x\mapsto \frac1{x^2}P_m\Big(\frac1{x}\Big),\quad m=0,1,\ldots ,2n-1, \]
</div>
<p> where \(P_m(t)\) are arbitrary selected algebraic polynomials in \(t\) of degree \(m\), i.e., </p>
<div class="equation" id="SystemEq">
<p>
  <div class="equation_content">
    \begin{equation} \label{SystemEq} \int _a^{+\infty }w(x) \frac1{x^2}P_m\Big(\frac1{x}\Big){\, \mathrm{d}}x=\sum _{k=1}^n \frac{A_k}{x_k^2}P_m\Big(\frac1{x_k}\Big), \quad m=0,1,\ldots ,2n-1. \end{equation}
  </div>
  <span class="equation_label">2.2</span>
</p>
</div>
<p><div class="unremark_thmwrapper " id="a0000000005">
  <div class="unremark_thmheading">
    <span class="unremark_thmcaption">
    Remark
    </span>
  </div>
  <div class="unremark_thmcontent">
  <p>Because of linearity, it is easy to see that this system of \(2n\) nonlinear equations in \(x_k\) and \(A_k\), \(k=1,\ldots ,n\), is equivalent to the corresponding system with monomials, i.e., when \(P_m(x)=x^m\), \(m=0,1,\ldots ,2n-1\).<span class="qed">â–¡</span></p>

  </div>
</div> </p>
<p>On the other side we consider the Gauss-Christoffel quadrature formula with respect to the weight function \(t\mapsto w(1/t)\) on \((0,1/a)\), i.e., </p>
<div class="equation" id="IntegFin">
<p>
  <div class="equation_content">
    \begin{equation} \label{IntegFin} \int _0^{1/a}w\Big(\frac1{t}\Big)g(t){\, \mathrm{d}}t= \sum _{k=1}^n B_k g(\tau _k)+R_n^G(g), \end{equation}
  </div>
  <span class="equation_label">2.3</span>
</p>
</div>
<p> where \(\tau _k\) and \(B_k\) are its nodes and Christoffel numbers, respectively, and \(R_n^G(g)\) is the corresponding remainder term. According to \((\ref{uslovi})\), such quadrature formulas exist uniquely, because the all moments \(\mu _k=\int _0^{1/a} w(1/t)t^k{\, \mathrm{d}}t\), \(k\ge 0\), exist and \(\mu _0{\gt}0\). </p>
<p>It is known that the nodes \(\tau _k\) in \((\ref{IntegFin})\) are eigenvalues of the following symmetric tridiagonal Jacobi matrix (cf. <span class="cite">
	[
	<a href="#gm-gvm-2008" >10</a>
	, 
	pp.
	
	325–328
	]
</span>) </p>
<div class="equation" id="Jac-matr">
<p>
  <div class="equation_content">
    \begin{equation} \label{Jac-matr} J_n=\left[\begin{array}{ccccc} \alpha _0& \sqrt{\beta _1}& & & {\mathbf{O}}\\[2mm] \sqrt{\beta _1}& \alpha _1& \sqrt{\beta _2}\\[2mm]& \sqrt{\beta _2}& \alpha _2& \ddots \\[2mm]& & \ddots & \ddots & \sqrt{\beta _{n-1}}\\[2mm] {\mathbf{O}}& & & \sqrt{\beta _{n-1}}& \alpha _{n-1} \end{array}\right], \end{equation}
  </div>
  <span class="equation_label">2.4</span>
</p>
</div>
<p> where \(\alpha _k\) and \(\beta _k\) are coefficients in the three-term recurrence relation </p>
<div class="displaymath" id="TTRR">
  \begin{eqnarray} \label{TTRR} \pi _{k+1}(t)& \! \! \! =\! \! \! & (t-\alpha _k)\pi _k(t)-\beta _k\pi _{k-1}(t),\quad k=0,1,\ldots ,\\[2mm] \pi _0(t)& \! \! \! =\! \! \! & 1,\  \  \pi _{-1}(t)=0,\nonumber \end{eqnarray}
</div>
<p> for the (monic) polynomials \(\{ \pi _k\} _{k\in \mathbb {N}_0}\) orthogonal with respect to the inner product </p>
<div class="equation" id="inner">
<p>
  <div class="equation_content">
    \begin{equation} \label{inner} (p,q)=\int _0^{1/a}w\Big(\frac1{t}\Big)p(t)q(t){\, \mathrm{d}}t. \end{equation}
  </div>
  <span class="equation_label">2.6</span>
</p>
</div>
<p> In fact, \(\pi _n(t)=(t-\tau _1)\cdots (t-\tau _n)\). </p>
<p>The weight coefficients \(B_k\) in (<a href="#IntegFin">2.3</a>) are given by </p>
<div class="displaymath" id="a0000000006">
  \[ B_k=\beta _0 v_{k,1}^2,\quad k=1,\ldots ,n, \]
</div>
<p> where \(v_{k,1}\) is the first component of the eigenvector \({\mathbf{v}}_k\  (=[v_{k,1}\  \ldots \  v_{k,n}]^{\mathrm T})\) corresponding to the eigenvalue \(\tau _k\) and normalized such that \({\mathbf{v}}_k^{\mathrm T}{\mathbf{v}}_k=1\), and \(\beta _0=\mu _0=\int _0^{1/a}w(1/t){\, \mathrm{d}}t\). </p>
<p>The most popular method for solving this eigenvalue problem is the Golub-Welsch procedure, obtained by a simplification of QR algorithm <span class="cite">
	[
	<a href="#GoWe" >7</a>
	]
</span>. This procedure is implemented in several packages including the most known <span class="tt">ORTPOL</span> given by Gautschi <span class="cite">
	[
	<a href="#gautschi3" >4</a>
	]
</span>. </p>
<p>As we can see from \((\ref{Jac-matr})\), for constructing Gauss–Christoffel quadratures \((\ref{IntegFin})\) for any number of nodes less than or equal to \(n\), we need the first \(n\) recursion coefficients \(\alpha _k\) and \(\beta _k\), \(k=0,1,\ldots ,n-1\), in \((\ref{TTRR})\). </p>
<p>In general, the recursion coefficients are known explicitly only for some narrow classes of orthogonal polynomialsc(e.g. for the classical orthogonal polynomials). In the case of the so-called <i class="it">strongly non-classical polynomials</i>, these recursion coefficients must be constructed numerically (cf. <span class="cite">
	[
	<a href="#GAU82" >3</a>
	]
</span>, <span class="cite">
	[
	<a href="#GauBook04" >5</a>
	]
</span>, <span class="cite">
	[
	<a href="#gm-gvm-2008" >10</a>
	, 
	pp.
	
	159–166
	]
</span>). However, recent progress in symbolic computation and variable-precision arithmetic today makes it possible to generate the recursive coefficients in \((\ref{TTRR})\) directly by using the original Chebyshev method of moments. Respectively symbolic/variable-precision software for orthogonal polynomials and Gaussian (and similar) quadratures is available. Our <i class="sc">Mathematica</i> package <span class="tt">OrthogonalPolynomials</span> (see <span class="cite">
	[
	<a href="#Aca" >1</a>
	]
</span> and <span class="cite">
	[
	<a href="#MilCvet-MB12" >13</a>
	]
</span>), is downloadable from the web site <span class="tt">http://www.mi.sanu.ac.rs/~gvm/</span>. Also, there is Gautschi’s software in <i class="sc">Matlab</i> (packages <span class="tt">OPQ</span> and <span class="tt">SOPQ</span>). </p>
<p>Now, we can give our main result: </p>
<p><div class="theorem_thmwrapper " id="a0000000007">
  <div class="theorem_thmheading">
    <span class="theorem_thmcaption">
    Theorem
    </span>
    <span class="theorem_thmlabel">2.1</span>
  </div>
  <div class="theorem_thmcontent">
  <p>Let \(x\mapsto w(x)\) be a weight function on \((a,+\infty )\), \(a{\gt}0\), such that the condition \((\ref{uslovi})\) holds. Assume also that \(\tau _k\) and \(B_k\), \(k=1,\ldots ,n\), are nodes and Christoffel numbers of the Gaussian quadrature formula \((\ref{IntegFin})\), respectively. Then there exists the generalized Gaussian quadrature formula </p>
<div class="equation" id="novaQ">
<p>
  <div class="equation_content">
    \begin{equation} \label{novaQ} \int _a^{+\infty }w(x) f(x){\, \mathrm{d}}x=\sum _{k=1}^n A_k f(x_k)+R_n(f), \end{equation}
  </div>
  <span class="equation_label">2.7</span>
</p>
</div>
<p> with </p>
<div class="equation" id="PARAM">
<p>
  <div class="equation_content">
    \begin{equation} \label{PARAM} x_k=\frac{1}{\tau _k}, \quad A_k=\frac{B_k}{\tau _k^2}>0,\quad k=1,\ldots ,n, \end{equation}
  </div>
  <span class="equation_label">2.8</span>
</p>
</div>
<p> which is exact for all functions of the form \(f(x)=x^{-2}P(x^{-1})\), where \(P\in {\EuScript P}_{2n-1}\). </p>
<p>The remainder term in this quadrature rule can be expressed in the following form \(R_n(f)=R_n^G(g)\), where \(g(t)=t^{-2}f(t^{-1})\). </p>

  </div>
</div> </p>
<p><div class="proof_wrapper" id="a0000000008">
  <div class="proof_heading">
    <span class="proof_caption">
    Proof
    </span>
    <span class="expand-proof">â–¼</span>
  </div>
  <div class="proof_content">
  
  </div>
</div> We start with the system of \(2n\) nonlinear equations \((\ref{SystemEq})\), whose solution determines the parameters of the quadrature formula \((\ref{novaQ})\). Our aim is to prove that this solution uniquely exists. </p>
<p>First, we take the sequence of orthogonal polynomials \(\{ \pi _m\} _{m=0}^{2n-1}\) in the system \((\ref{SystemEq})\) and then by a simple change of variables \(x=1/t\) in the integral on the left hand side we obtain </p>
<div class="displaymath" id="a0000000009">
  \begin{eqnarray*}  \int _a^{+\infty }w(x) \frac1{x^2}\pi _m\Big(\frac1{x}\Big){\, \mathrm{d}}x& =&  \int _0^{1/a}w\Big(\frac1{t}\Big)\pi _m(t){\, \mathrm{d}}t\\[1mm]& =& (\pi _0,\pi _m)\\[1mm]& =& \mu _0\delta _{0,m}, \end{eqnarray*}
</div>
<p> where the inner product is defined by \((\ref{inner})\) and \(\delta _{k,m}\) is Kronecker’s delta. </p>
<p>Evidently, this leads to the system of equations </p>
<div class="equation" id="sys1">
<p>
  <div class="equation_content">
    \begin{equation} \label{sys1} \sum _{k=1}^n A_k \frac{1}{x_k^2}\pi _m\Big(\frac1{x_k}\Big)= \mu _0\delta _{0,m},\quad m=0,1,\ldots ,2n-1, \end{equation}
  </div>
  <span class="equation_label">2.9</span>
</p>
</div>
<p> but, by an application of the Gaussian rule \((\ref{IntegFin})\), it gives also another system of equations </p>
<div class="equation" id="sys2">
<p>
  <div class="equation_content">
    \begin{equation} \label{sys2} \sum _{k=1}^n B_k \pi _m(\tau _k)=\mu _0\delta _{0,m},\quad m=0,1,\ldots ,2n-1, \end{equation}
  </div>
  <span class="equation_label">2.10</span>
</p>
</div>
<p> because \(R_n^G(g)=0\) for each \(g\in {\EuScript P}_{2n-1}\). The last system has the unique solution, and it represents the parameters \(\tau _k\) and \(B_k\), \(k=1,\ldots ,n\), of the Gaussian quadrature \((\ref{IntegFin})\). </p>
<p>Since the systems of equations \((\ref{sys1})\) and \((\ref{sys2})\) are equivalent, the statement of this theorem follows directly. <div class="proof_wrapper" id="a0000000010">
  <div class="proof_heading">
    <span class="proof_caption">
    Proof
    </span>
    <span class="expand-proof">â–¼</span>
  </div>
  <div class="proof_content">
  
  </div>
</div> </p>
<h1 id="SpecialC">3  Special cases and numerical examples</h1>

<p>In this section we consider special cases of quadrature formulas with respect to the weight function \(x\mapsto w(x)=x^\beta \log ^m x\), where \(0\le \beta {\lt}1\) and \(m\in \mathbb {N}_0\). For \(\beta =0\) and \(m=0\), it reduces to the constant weight \(w(x)=1\). In order to show the efficiency of the obtained quadrature formulas we present a few numerical examples. </p>
<p>We start this section with the weight function \(x\mapsto w(x)=x^\beta \), \(0\le \beta {\lt}1\). </p>
<p>The condition \((\ref{uslovi})\) is satisfied, because </p>
<div class="displaymath" id="a0000000011">
  \[ \int _a^{+\infty } \frac{w(x)}{x^2}{\, \mathrm{d}}x=\frac{a^{\beta -1}}{1-\beta }. \]
</div>
<p>Here we consider only the case \(\beta =0\), i.e., when \(w(x)=1\). Since </p>
<div class="displaymath" id="a0000000012">
  \[ \int _a^{+\infty }f(x){\, \mathrm{d}}x=a\int _1^{+\infty }f(ax){\, \mathrm{d}}x, \]
</div>
<p> we see that for this important case the following statement holds. </p>
<p><div class="proposition_thmwrapper " id="propKOEF">
  <div class="proposition_thmheading">
    <span class="proposition_thmcaption">
    Proposition
    </span>
    <span class="proposition_thmlabel">3.1</span>
  </div>
  <div class="proposition_thmcontent">
  <p> Let </p>
<div class="equation" id="const">
<p>
  <div class="equation_content">
    \begin{equation} \label{const} \int _a^{+\infty }f(x){\, \mathrm{d}}x=\sum _{k=1}^n A_k(a)f(x_k(a))+R_n(f;a),\quad a>0, \end{equation}
  </div>
  <span class="equation_label">3.1</span>
</p>
</div>
<p> be a generalized Gaussian quadrature \((\ref{novaQ})\) \((\)with the constant weight function \(w(x)=1)\). Then </p>
<div class="displaymath" id="a0000000013">
  \[ A_k(a)=aA_k(1)\quad \mbox{and}\quad x_k(a)=ax_k(1),\  \  k=1,\ldots ,n. \]
</div>

  </div>
</div> </p>
<p>This means that it is enough to know only quadrature parameters for \(a=1\). These parameters can be obtained directly using \((\ref{PARAM})\) and Gauss-Legendre parameters \(\tau _k\) and \(B_k\) for transformed interval \((0,1)\). </p>
<p>Recursive coefficients in \((\ref{TTRR})\), in this case for translated monic Legendre polynomials, are </p>
<div class="displaymath" id="a0000000014">
  \[ \alpha _k=\frac12,\  \  k\ge 0,\quad \beta _0=1,\  \  \beta _k=\frac{k^2}{4(4k^2-1)},\  \  k\ge 1. \]
</div>
<p> Otherwise, it can be obtained using our <i class="sc">Mathematica</i> Package <span class="tt">OrthogonalPolynomials</span> in symbolic form (see <span class="cite">
	[
	<a href="#Aca" >1</a>
	]
</span> and <span class="cite">
	[
	<a href="#MilCvet-MB12" >13</a>
	]
</span>). For example, if we need the first forty recurrence coefficients, then we start with the first eighty moments \(\mu _k=1/(k+1)\), \(k=0,1,\ldots ,79\), and then we use the standard Chebyshev algorithm (cf. <span class="cite">
	[
	<a href="#gm-gvm-2008" >10</a>
	, 
	160–162
	]
</span>: <span class="tt"><pre class="verbatim">
  &lt;&lt; orthogonalPolynomials`
  mom=Table[1/(k+1), {k,0,79}];
  {al,be} = aChebyshevAlgorithm[mom, Algorithm -&gt; Symbolic]
</pre></span> These recursive coefficients enable us to construct quadrature formulas \((\ref{IntegFin})\) for any number of nodes up to \(40\). </p>
<p>However, in this Legendre case (translated to \((0,1)\)) we can directly use <span class="tt">aGaussianNodesWeights</span> routine to construct nodes and weights in the Gaussian quadrature formula \((\ref{IntegFin})\), as well as ones in the quadrature formula \((\ref{novaQ})\): <span class="tt"><pre class="verbatim">
  &lt;&lt; orthogonalPolynomials`
  transLeg[n_] := (aGaussianNodesWeights[n, {aLegendre},
     WorkingPrecision -&gt; 70, Precision -&gt; 65] + {1,0})/2;
     parQF = Table[transLeg[n], {n,2,40,2}];
     For[m = 1, m &lt; 21, m++,
        parQF[[m]][[2]] = parQF[[m]][[2]]/parQF[[m]][[1]]^2;
        parQF[[m]][[1]] = 1/parQF[[m]][[1]];]
</pre></span> Thus, in this way for \(a=1\), we obtain quadrature parameters \(x_k\) and \(A_k\) for each \(n=2(2)40\). </p>
<p><div class="example_thmwrapper " id="a0000000015">
  <div class="example_thmheading">
    <span class="example_thmcaption">
    Example
    </span>
    <span class="example_thmlabel">3.2</span>
  </div>
  <div class="example_thmcontent">
  <p>In order to show the efficiency of our quadrature rule \((\ref{novaQ})\) we apply it to the integral </p>
<div class="displaymath" id="a0000000016">
  \[ J(a;c)=\int _a^{+\infty }\frac{1}{(x-2)^2+c^2}{\, \mathrm{d}}x =\frac{1}{2 c}\left[\pi -2 \arctan \big(\frac{a-2}{c}\big)\right], \]
</div>
<p> for different values of \(a{\gt}0\) and \(c{\gt}0\). In Figure <a href="#figures12">3.1</a> we present graphics of the function </p>
<div class="displaymath" id="a0000000017">
  \[  x\mapsto f(x;c)=\frac1{(x-2)^2+c^2}  \]
</div>
<p> for \(c=\frac18,\frac14,\frac12\), and \(1\), as well as the corresponding graphics of the exact values of this integral \(J(a;c)\) (right). </p>
<p>In order to test the quadrature formula \((\ref{const})\), we apply it to \(J(a;1)\) for \(a=\frac12\), \(1\), \(2\), \(3\), \(4\), and \(8\), when \(n=2(2)40\). </p>
<figure id="figures12">
  <div class="centered"> <div class="displaymath" id="a0000000018">
  \[  \includegraphics[width=0.50\textwidth ]{GraffunBOJA.pdf}\includegraphics[width=0.50\textwidth ]{TVRintBOJA.pdf} \]
</div> </div>
<figcaption>
  <span class="caption_title">Figure</span> 
  <span class="caption_ref">3.1</span> 
  <span class="caption_text">The function \(x\mapsto f(x;c)\) (left) and the integral \(a\mapsto J(a;c)\) (right) for \(c=\frac18\) (blue line), \(c=\frac14\) (black line), \(c=\frac12\) (brown line), and \(c=1\) (red line).</span> 
</figcaption>


</figure>
<p>Relative errors in the quadrature sums </p>
<div class="displaymath" id="a0000000019">
  \[ Q_n(f(\cdot ;c);a)=\sum _{k=1}^n A_k(a)f(x_k(a);c), \]
</div>
<p> defined by </p>
<div class="displaymath" id="a0000000020">
  \[ \operatorname {err}_n(f(\cdot ;c);a)=\Bigl|\frac{Q_n(f(\cdot ;c);a)-J(a;c)}{J(a;c)}\Bigr|,  \]
</div>
<p> are displayed in Figure&#160;<a href="#figure3">3.2</a> in a log-scale. </p>
<figure id="figure3">
  <div class="centered"> <div class="displaymath" id="a0000000021">
  \[  \includegraphics[width=0.80\textwidth ]{err_graf1.pdf} \]
</div> </div>
<figcaption>
  <span class="caption_title">Figure</span> 
  <span class="caption_ref">3.2</span> 
  <span class="caption_text">Relative errors \(\operatorname {err}_n(f(\cdot ;1);a)\) in quadrature sums \(Q_n(f(\cdot ;1);a)\).</span> 
</figcaption>


</figure>
<p> Numerical results show that the convergence is much faster if the parameter \(a\) is larger. For example, if \(a = 2\), then for \(n=10(10)40\), the relative errors are \(1.71\times 10^{-7}\), \(1.83\times 10^{-14}\), \(1.91\times 10^{-21}\), \(1.94\times 10^{-28}\), respectively, while the corresponding errors for \(a=4\) are \(5.52\times 10^{-15}\), \(1.21\times 10^{-29}\), \(1.40\times 10^{-44}\), \(1.44\times 10^{-59}\). </p>
<p>Otherwise, this integrand \(f(x;c)\) has poles at the points \(2\pm {\mathrm{i}}c\), which are approaching the real line when \(c\) tends to zero. In this case, for small values of \(a\) (near \(2\) or less than \(2\)), the convergence of the quadrature process slows down considerably, because of a strong influence of these singularities. This effect can be seen from Table&#160;<a href="#TableEx1">3.1</a>, where quadrature approximations and corresponding relative errors are presented for \((a,c)=(1,\frac14)\), \((\frac{21}{10},10^{-6})\), and \((4,10^{-6})\). In order to save space, in last case only relative errors are given. Digits in error are underlined, and numbers in parenthesis indicate the decimal exponents. </p>
<p>Notice that the integral \(J(a;0)\) for \(a\le 2\) does not exist. </p>
<p>Finally, the last column shows that for \((a,c)=(4,10^{-6})\), the convergence of the quadrature rule \((\ref{const})\) is very fast. </p>
<div class="table"  id="TableEx1">
   <div class="centered"><small class="small"><table class="tabular">
  <tr>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:right; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(n\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="2">
      <p> \((a,c)=(1,1/4)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="2">
      <p> \((a,c)=(21/10,10^{-6})\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \((a,c)=(4,10^{-6})\)</p>

    </td>
  </tr>
  <tr>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:right; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(2\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:right; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(\underline{2.83088}\  \) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(7.56(-1)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(\underline{4.21706255691703}\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(5.78(-1)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(5.92(-3)\  \, \)</p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:right; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(4\) </p>

    </td>
    <td  style="text-align:right; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(\underline{5.38719}\  \) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(5.35(-1)\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(\underline{8.01223217799471}\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.99(-1)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(9.70(-6)\  \, \)</p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:right; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(6\) </p>

    </td>
    <td  style="text-align:right; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(\underline{7.41379}\  \) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(3.60(-1)\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(9.\underline{47887835712778}\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(5.21(-2)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.24(-8)\  \, \)</p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:right; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(8\) </p>

    </td>
    <td  style="text-align:right; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(\underline{8.88711}\  \) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(2.33(-1)\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(9.\underline{88043864297441}\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.20(-2)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.42(-11)\)</p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:right; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(10\) </p>

    </td>
    <td  style="text-align:right; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(\underline{9.89102}\  \) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.46(-1)\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(9.9\underline{7447558340612}\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(2.55(-3)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.53(-14)\)</p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:right; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(20\) </p>

    </td>
    <td  style="text-align:right; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(11.\underline{45438}\  \) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.14(-2)\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(9.99999\underline{276505451}\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(7.23(-7)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.47(-29)\)</p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:right; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(30\) </p>

    </td>
    <td  style="text-align:right; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(11.5\underline{7808}\  \) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(7.23(-4)\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(9.99999999\underline{813998}\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.53(-10)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.08(-44)\) </p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:right; border-right:1px solid black; border-left:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p>\(40\) </p>

    </td>
    <td  style="text-align:right; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(11.586\underline{06}\  \) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(3.41(-5)\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(9.999999999666\underline{38}\) </p>

    </td>
    <td  style="text-align:left; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(2.86(-14)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(6.99(-60)\) </p>

    </td>
  </tr>
</table> </small><figcaption>
  <span class="caption_title">Table</span> 
  <span class="caption_ref">3.1</span> 
  <span class="caption_text">Quadrature sums \(Q_n(f(\cdot ;c);a)\) and their relative errors &#8195;&#8195;&#8195;&#8195;&#8195;&#8195;&#8195;&#8195;<br />\(err_n(f(\cdot ;c);a)\) for integrals \(J(a;c)\).</span> 
</figcaption> </div>
</div>

  </div>
</div> </p>
<p>In the sequel we consider quadrature rules with respect to the weight function \(x\mapsto w(x)=x^\beta \log x\), \(0\le \beta {\lt}1\). Here we suppose that \(a\ge 1\). The condition \((\ref{uslovi})\) is satisfied, because </p>
<div class="equation" id="U1">
<p>
  <div class="equation_content">
    \begin{equation} \label{U1} 0<\int _a^{+\infty }\frac{w(x)}{x^2}{\, \mathrm{d}}x= \frac{a^{\beta -1}}{(1-\beta )^2}\left[1+(1-\beta )\log a\right]. \end{equation}
  </div>
  <span class="equation_label">3.2</span>
</p>
</div>
<p> In this case, the moments </p>
<div class="displaymath" id="a0000000022">
  \[ \mu _k=\int _0^{1/a}w(1/t)t^k{\, \mathrm{d}}t=\int _0^{1/a}t^{k-\beta }\log \frac1t{\, \mathrm{d}}t \]
</div>
<p> can be expressed in the form </p>
<div class="equation" id="Mom1">
<p>
  <div class="equation_content">
    \begin{equation} \label{Mom1} \mu _k=\frac{a^{\beta -k-1} [(k+1-\beta )\log a+1]}{(k+1-\beta )^2},\quad k\ge 0. \end{equation}
  </div>
  <span class="equation_label">3.3</span>
</p>
</div>
<p>Taking the first one hundred moments (<span class="tt">mom</span>) (e.g. for \(a=1\) and \(\beta =1/4\)) and using <i class="sc">Mathematica</i> Package <span class="tt">OrthogonalPolynomials</span>, we can get the first fifty recurrence coefficients \(\alpha _k\) and \(\beta _k\) (denoted by <span class="tt">{al1,be1}</span>) in the three-term recurrence relation (<a href="#TTRR">2.5</a>) in a symbolic form <span class="tt"><pre class="verbatim">
  &lt;&lt; orthogonalPolynomials`
  mom=Table[(a^(-1+b-k)(1+(1-b+k)Log[a]))/(1-b+k)^2, {k,0,99}];
  mom1=mom/. {a-&gt;1, b-&gt;1/4}
  {al1,be1} = aChebyshevAlgorithm[mom, Algorithm -&gt; Symbolic]
</pre></span> For example, first four coefficients are <small class="small"><div class="displaymath" id="a0000000023">
  \[ \alpha _0=\frac{9}{49},\, \alpha _1=\frac{209897}{452025},\,  \alpha _2=\frac{6582284926939}{13538179995075},\,  \alpha _3=\frac{7618613698603068100869609}{15464687102113919816429449}  \]
</div></small> and <small class="small"><div class="displaymath" id="a0000000024">
  \[ \beta _0=\frac{16}{9},\   \, \beta _1=\frac{11808}{290521},\   \, \beta _2=\frac{213147564896}{3717280400625},\   \, \beta _3=\frac{421267942813254097088}{6997413354065613077481}. \]
</div></small> </p>
<p><div class="example_thmwrapper " id="a0000000025">
  <div class="example_thmheading">
    <span class="example_thmcaption">
    Example
    </span>
    <span class="example_thmlabel">3.3</span>
  </div>
  <div class="example_thmcontent">
  <p>As a test example we consider the function </p>
<div class="displaymath" id="a0000000026">
  \[  x\mapsto f(x)=\frac1{(x+1)^{2}}  \]
</div>
<p> and integral (see <span class="cite">
	[
	<a href="#Evans2005" >2</a>
	]
</span>) <small class="small"><div class="displaymath" id="a0000000027">
  \begin{align*}  I(f;a)=& \int _a^{+\infty }\frac{x^{1/4}\log x}{(x+1)^2}{\, \mathrm{d}}x\\ =& \frac{1}{36 a^{\frac74} (a+1)} \left\{ -9 (a+1)\Phi \Bigl(-\frac{1}{a},2,\frac{7}{4}\Bigr)\right.\qquad \quad \\[2mm]& \qquad \left.+4 a \left[3 (a+1)(\log a+4) \,  _2F_1\Bigl(\frac{3}{4},1;\frac{7}{4};-\frac{1}{a}\Bigr) +4 a+9 a \log a+4\right]\right\} , \end{align*}
</div></small> </p>
<p>where \(\Phi \) and \(_2F_1\) are the Lerch transcendent and Gauss hypergeometric function, defined by </p>
<div class="displaymath" id="a0000000028">
  \[ \Phi (z,s,a)=\sum _{k=0}^{+\infty } \frac{z^k}{(k+a)^s}\quad \mbox{and}\quad _2F_1(a,b;c;z)=\sum _{k=0}^{+\infty }\frac{(a)_k(b)_k}{(c)_k}\frac{z^k}{k!}, \]
</div>
<p> respectively, and \((a)_k=a(a+1)\ldots (a+k-1)\) is the Pochhammer symbol. </p>
<p>We consider this integral for two values of the lower bound: \(a=1\) and \(a={\mathrm{e}}\), i.e., <small class="small"><div class="displaymath" id="a0000000029">
  \[  I(f;1)={1.35974328097600895}\ldots \  \  \mbox{and}\   I(f;\mathrm{e})={1.22897618668037255}\ldots \  . \]
</div></small> Applying Gauss-Laguerre rule to \(I(f;1)\) (translated from \((1,+\infty )\) to \((1,+\infty )\)) gives poor results. Relative errors in the corresponding Gauss-Laguerre quadrature sums are presented in Table <a href="#Tab:one">3.2</a>. As we can see only two two decimal digits are true in quadrature sum with \(2048\) nodes! </p>
<div class="table"  id="Tab:one">
   <div class="centered"><small class="small"><table class="tabular">
  <tr>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(n=2\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=8\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=32\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=128\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=512\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=2048\)</p>

    </td>
  </tr>
  <tr>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-left:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p>\(6.72(-1)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(3.60(-1)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(1.64(-1)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(7.00(-2)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(2.90(-2)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(1.18(-2)\) </p>

    </td>
  </tr>
</table> </small><figcaption>
  <span class="caption_title">Table</span> 
  <span class="caption_ref">3.2</span> 
  <span class="caption_text">Relative errors in Gauss-Laguerre quadrature sums &#8195;&#8195;&#8195;&#8195;&#8195;&#8195;&#8195;&#8195;<br />with \(n=2,8,32,128,512\) and \(2048\) nodes.</span> 
</figcaption> </div>
</div>
<p>Now, we apply our quadrature formula \((\ref{novaQ})\), with parameters given by \((\ref{PARAM})\), to \(I(f;1)\) and \(I(f;\mathrm{e})\), with only \(n=2(2)12\) nodes. The relative errors in the quadrature sums \(Q_n(f;a)=\sum _{k=1}^n A_kf(x_k)\), </p>
<div class="displaymath" id="a0000000030">
  \[ \operatorname {err}_n(f;a)=\Bigl|\frac{Q_{n}(f;a)-I(f;a)}{I(f;a)}\Bigr|,  \]
</div>
<p> are presented in Table&#160;<a href="#Tab:two">3.3</a>. </p>
<div class="table"  id="Tab:two">
   <div class="centered"><small class="small"><table class="tabular">
  <tr>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(a\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=2\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=4\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=6\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=8\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=10\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=12\)</p>

    </td>
  </tr>
  <tr>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(1\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(2.94(-3)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(4.24(-6)\  \, \) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(5.15(-9)\  \) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(5.72(-12)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(4.74(-13)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(7.07(-13)\) </p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:center; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\({\mathrm{e}}\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(2.40(-4)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.64(-8)\  \, \) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(\  8.91(-13)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(8.83(-14)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(5.31(-14)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(3.80(-14)\) </p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:center; border-right:1px solid black; border-left:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p>\({\mathrm{e}}^2\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(7.18(-6)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(1.28(-11)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(\  3.10(-14)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">&nbsp;</td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">&nbsp;</td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">&nbsp;</td>
  </tr>
</table> </small><figcaption>
  <span class="caption_title">Table</span> 
  <span class="caption_ref">3.3</span> 
  <span class="caption_text">Relative errors \(\operatorname {err}_n(f;a)\) in quadrature sums \(Q_n(f;a)\) &#8195;&#8195;&#8195;&#8195;&#8195;&#8195;&#8195;&#8195;<br />for different number of nodes \(n\) and three values of \(a\) \((=1,{\mathrm{e}},\  \mbox{and}\  {\mathrm{e}}^2)\).</span> 
</figcaption>  </div>
</div>
<p>As we can see, the convergence is faster when \(a\) is bigger. In the third line of the same table we also present the corresponding relative errors when \(a={\mathrm{e}}^2\) and \(n=2,4\), and \(6\). </p>

  </div>
</div> </p>
<p>Finally, we mention that this approach can be applied also in the case of the weight functions </p>
<div class="displaymath" id="a0000000031">
  \[ w(x)=w_m(x)=x^\beta \log ^m x,\quad 0\le \beta {\lt}1,\  m=2,3,\ldots \  , \]
</div>
<p> on the interval \((a,+\infty )\), with \(a\ge 1\). </p>
<p>The condition \((\ref{uslovi})\) for \(U_m=\int _a^{+\infty }x^{-2}w_m(x){\, \mathrm{d}}x\) is also satisfied, because </p>
<div class="displaymath" id="a0000000032">
  \[ U_m=\frac1{1-\beta }\left[mU_{m-1}+a^{\beta -1}\log ^m a\right],\quad m=2,3,\ldots \, , \]
</div>
<p> where \(U_1\) is given in \((\ref{U1})\). The corresponding moments </p>
<div class="displaymath" id="a0000000033">
  \[ \mu _k^{[m]}=\int _0^{1/a}t^{k-\beta }\log ^m\frac1{t}{\, \mathrm{d}}t,\quad k=0,1,\ldots \, , \]
</div>
<p> can be expressed recursively in terms of the moments \(\mu _k^{[m-1]}\), </p>
<div class="displaymath" id="a0000000034">
  \[ \mu _k^{[m]}=\frac{1}{k+1-\beta }\left(m \mu _k^{[m-1]}+a^{\beta -k-1}\log ^m a\right), \quad m=2,3,\ldots \, , \]
</div>
<p> where the moments \(\mu _k^{[1]}\  (\equiv \mu _k)\) are given by \((\ref{Mom1})\). For example, for \(m=2\) we get </p>
<div class="displaymath" id="a0000000035">
  \[ \mu _k^{[2]}=\frac{a^{\beta -k-1} \left[(k+1-\beta )^2\log ^2a +2(k+1-\beta ) \log a+2\right]}{(k+1-\beta )^3},\quad k\ge 0. \]
</div>
<p>The corresponding recursive coefficients, for example for \(\beta =0\) and \(a=1\), are <small class="small"><div class="displaymath" id="a0000000036">
  \[ \alpha _0=\frac{1}{8},\, \alpha _1=\frac{115}{296},\,  \alpha _2=\frac{28200187}{62721512},\,  \alpha _3=\frac{28003451041760695}{59414538084233528}, \  \ldots  \]
</div></small> and <small class="small"><div class="displaymath" id="a0000000037">
  \[ \beta _0=2,\, \beta _1=\frac{37}{1728},\  \, \beta _2=\frac{211897}{4620375},\   \, \beta _3=\frac{945381680572419}{17600932734728000},\  \ldots  . \]
</div></small> </p>
<p><div class="example_thmwrapper " id="a0000000038">
  <div class="example_thmheading">
    <span class="example_thmcaption">
    Example
    </span>
    <span class="example_thmlabel">3.4</span>
  </div>
  <div class="example_thmcontent">
  <p>We consider the function \(x\mapsto 1/(1+x^2)\) and the corresponding weighted integral over \((a,+\infty )\) </p>
<div class="displaymath" id="a0000000039">
  \[ I(f;a)=\int _a^{+\infty }\frac{\log ^2x}{1+x^2}{\, \mathrm{d}}x, \]
</div>
<p> for two values of \(a\) \((=1\) and \(={\mathrm{e}})\), for which <small class="small"><div class="displaymath" id="a0000000040">
  \[ I(f;1)={1.93789229251873876096726969169}\ldots  \]
</div></small> and <small class="small"><div class="displaymath" id="a0000000041">
  \[ \  \  I(f;{\mathrm{e}})={1.80988687939786942602016447246}\ldots \  . \]
</div></small> </p>
<p>Applying our quadrature formula \((\ref{novaQ})\), with parameters given by \((\ref{PARAM})\), to \(I(f;1)\) and \(I(f;\mathrm{e})\), with \(n=2(2)12\) nodes, we get the corresponding quadrature approximations \(Q_n(f;a)\) with the relative errors \(\operatorname {err}_n(f;a)\) presented in Table&#160;<a href="#Tab:three">3.4</a>. </p>
<div class="table"  id="Tab:three">
   <div class="centered"><small class="small"><table class="tabular">
  <tr>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(a\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=2\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=4\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=6\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=8\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=10\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(n=12\)</p>

    </td>
  </tr>
  <tr>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black; border-left:1px solid black" 
        rowspan=""
        colspan="">
      <p>\(1\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.66(-4)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.31(-6)\  \  \) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(1.98(-10)\,  \) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(5.73(-12)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(2.08(-15)\) </p>

    </td>
    <td  style="border-top-style:solid; border-top-color:black; border-top-width:1px; text-align:center; border-right:1px solid black" 
        rowspan=""
        colspan="">
      <p> \(2.56(-17)\) </p>

    </td>
  </tr>
  <tr>
    <td  style="text-align:center; border-right:1px solid black; border-left:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p>\({\mathrm{e}}\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(5.33(-5)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(5.04(-10)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(1.86(-13)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(2.05(-17)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(1.22(-21)\) </p>

    </td>
    <td  style="text-align:center; border-right:1px solid black; border-bottom-style:solid; border-bottom-color:black; border-bottom-width:1px" 
        rowspan=""
        colspan="">
      <p> \(3.30(-26)\) </p>

    </td>
  </tr>
</table> </small><figcaption>
  <span class="caption_title">Table</span> 
  <span class="caption_ref">3.4</span> 
  <span class="caption_text">Relative errors \(\operatorname {err}_n(f;a)\) in quadrature sums \(Q_n(f;a)\) &#8195;&#8195;&#8195;&#8195;&#8195;&#8195; <br />for different number of nodes \(n\) and two values of \(a\) \((=1\) and \(={\mathrm{e}})\).</span> 
</figcaption> </div>
</div>
<p>Here also we can note a faster convergence when \(a\) has a larger value. </p>

  </div>
</div> </p>
<p><small class="footnotesize">  </small></p>
<div class="bibliography">
<h1>Bibliography</h1>
<dl class="bibliography">
  <dt><a name="Aca">1</a></dt>
  <dd><p><i class="sc">A.S. Cvetković</i> and <i class="sc">G.V. Milovanović</i>, <i class="it">The Mathematica Package “OrthogonalPolynomials”</i>, Facta Univ. Ser. Math. Inform., <b class="bf">19</b>, pp. 17–36, 2004. </p>
</dd>
  <dt><a name="Evans2005">2</a></dt>
  <dd><p><a href ="http://dx.doi.org/10.1080/00207160512331323399"> <i class="sc">G.A. Evans</i>, <i class="it">Some new thoughts on Gauss-Laguerre quadrature</i>, Int. J. Comput. Math. <b class="bf">82</b>, pp. 721–730, 2005. <img src="img-0001.png" alt="\includegraphics[scale=0.1]{ext-link.png}" style="width:12.0px; height:10.700000000000001px" />
</a> </p>
</dd>
  <dt><a name="GAU82">3</a></dt>
  <dd><p><a href ="http://dx.doi.org/10.1137/0903018"> <i class="sc">W. Gautschi</i>, <i class="it">On generating orthogonal polynomials</i>, SIAM J. Sci. Statist. Comput., <b class="bf">3</b>, pp. 289–317, 1982. <img src="img-0001.png" alt="\includegraphics[scale=0.1]{ext-link.png}" style="width:12.0px; height:10.700000000000001px" />
</a> </p>
</dd>
  <dt><a name="gautschi3">4</a></dt>
  <dd><p><a href ="http://dx.doi.org/10.1145/174603.174605"> <i class="sc">W. Gautschi</i>, <i class="it">Algorithm 726: ORTHPOL – A package of routines for generating orthogonal polynomials and Gauss-type quadrature rules</i>, ACM Trans. Math. Software, <b class="bf">20</b>, pp. 21–62, 1994. <img src="img-0001.png" alt="\includegraphics[scale=0.1]{ext-link.png}" style="width:12.0px; height:10.700000000000001px" />
</a> </p>
</dd>
  <dt><a name="GauBook04">5</a></dt>
  <dd><p><i class="sc">W. Gautschi</i>, <i class="it">Orthogonal Polynomials: Computation and Approximation</i>, Clarendon Press, Oxford, 2004. </p>
</dd>
  <dt><a name="Gautschi2">6</a></dt>
  <dd><p><a href ="http://dx.doi.org/10.1007/s11075-010-9366-0"> <i class="sc">W. Gautschi</i>, <i class="it">Gauss quadrature routines for two classes of logarithmic weight functions</i>, Numer. Algor. <b class="bf">55</b>, pp. 265–277, 2010. <img src="img-0001.png" alt="\includegraphics[scale=0.1]{ext-link.png}" style="width:12.0px; height:10.700000000000001px" />
</a> </p>
</dd>
  <dt><a name="GoWe">7</a></dt>
  <dd><p><a href ="http://dx.doi.org/10.2307/2004418"> <i class="sc">G. Golub</i> and <i class="sc">J.H. Welsch</i>, <i class="it">Calculation of Gauss quadrature rules</i>, Math. Comp. <b class="bf">23</b>, pp. 221–230, 1969. <img src="img-0001.png" alt="\includegraphics[scale=0.1]{ext-link.png}" style="width:12.0px; height:10.700000000000001px" />
</a> </p>
</dd>
  <dt><a name="Markov1885">8</a></dt>
  <dd><p><a href ="http://dx.doi.org/10.1007/BF01443287"> <i class="sc">A. Markov</i>, <i class="it">Sur la méthode de Gauss pour le calcul approché des intégrales</i>, Math. Ann., <b class="bf">25</b>, pp.&#160;427–432, 1885. <img src="img-0001.png" alt="\includegraphics[scale=0.1]{ext-link.png}" style="width:12.0px; height:10.700000000000001px" />
</a> </p>
</dd>
  <dt><a name="gm-gm-2003">9</a></dt>
  <dd><p><a href ="http://dx.doi.org/10.1137/S0036142901391475"> <i class="sc">G. Mastroianni</i> and <i class="sc">G. Monegato</i>, <i class="it">Truncated quadrature rules over \((0, \infty )\) and Nyström-type methods</i>, SIAM J. Numer. Anal. <b class="bf">41</b>, pp. 1870–1892, 2003. <img src="img-0001.png" alt="\includegraphics[scale=0.1]{ext-link.png}" style="width:12.0px; height:10.700000000000001px" />
</a> </p>
</dd>
  <dt><a name="gm-gvm-2008">10</a></dt>
  <dd><p><i class="sc">G. Mastroianni</i> and <i class="sc">G.V. Milovanović</i>, <i class="it">Interpolation Processes – Basic Theory and Applications</i>, Springer-Verlag, Berlin – Heidelberg – New York, 2008. </p>
</dd>
  <dt><a name="Stud15">11</a></dt>
  <dd><p><i class="sc">G.V. Milovanović</i>, <i class="it">Construction and applications of Gaussian quadratures with nonclassical and exotic weight function</i>, Stud. Univ. Babeş-Bolyai Math., <b class="bf">60</b>, pp. 211–233, 2015. </p>
</dd>
  <dt><a name="FILOMAT15">12</a></dt>
  <dd><p><i class="sc">G.V. Milovanović</i>, <i class="it">Generalized Gaussian quadratures for integrals with logarithmic singularity</i>, FILOMAT (to appear). </p>
</dd>
  <dt><a name="MilCvet-MB12">13</a></dt>
  <dd><p><i class="sc">G.V. Milovanović</i> and <i class="sc">A.S. Cvetković</i>, <i class="it">Special classes of orthogonal polynomials and corresponding quadratures of Gaussian type</i>, Math. Balkanica, <b class="bf">26</b>, pp. 169–184, 2012. </p>
</dd>
  <dt><a name="JCAM15">14</a></dt>
  <dd><p><a href ="http://dx.doi.org/10.1016/j.cam.2014.10.009"> <i class="sc">G.V. Milovanović, T.S. Igić</i> and <i class="sc">D. Turnić</i>, <i class="it">Generalized quadrature rules of Gaussian type for numerical evaluation of singular integrals</i>, J. Comput. Appl. Math. <b class="bf">278</b>, pp. 306–25, 2015. <img src="img-0001.png" alt="\includegraphics[scale=0.1]{ext-link.png}" style="width:12.0px; height:10.700000000000001px" />
</a> </p>
</dd>
  <dt><a name="Posse">15</a></dt>
  <dd><p><i class="sc">C. Posse</i>, <i class="it">Sur les quadratures</i>, Nouv. Ann. Math. (2) <b class="bf">14</b>, pp. 49–62, 1875. </p>
</dd>
  <dt><a name="Stieltjes1884">16</a></dt>
  <dd><p><i class="sc">T.J. Stieltjes</i>, <i class="it">Quelques recherches sur la théorie des quadratures dites mécaniques</i>, Ann. Sci. Éc. Norm. Paris, Sér. 2, <b class="bf">1</b>, pp. 409–426, 1884. </p>
</dd>
  <dt><a name="Usp">17</a></dt>
  <dd><p><a href ="http://dx.doi.org/10.1090/S0002-9947-1928-1501444-8"> <i class="sc">J.V. Uspensky</i>, <i class="it">On the convergence of quadrature formulas related to an infinite interval</i>, Trans. Amer. Math. Soc., <b class="bf">30</b>, pp. 542–559, 1928. <img src="img-0001.png" alt="\includegraphics[scale=0.1]{ext-link.png}" style="width:12.0px; height:10.700000000000001px" />
</a> </p>
</dd>
  <dt><a name="Xu_GVM">18</a></dt>
  <dd><p><i class="sc">Zhenhua Xu</i> and <i class="sc">G.V. Milovanović</i>, <i class="it">Efficient method for the computation of oscillatory Bessel transform and Bessel Hilbert transform</i> (submitted). </p>
</dd>
</dl>


</div>
</div> <!--main-text -->
</div> <!-- content-wrapper -->
</div> <!-- content -->
</div> <!-- wrapper -->

<nav class="prev_up_next">
</nav>

<script type="text/javascript" src="/var/www/clients/client1/web1/web/files/jnaat-files/journals/1/articles/js/jquery.min.js"></script>
<script type="text/javascript" src="/var/www/clients/client1/web1/web/files/jnaat-files/journals/1/articles/js/plastex.js"></script>
<script type="text/javascript" src="/var/www/clients/client1/web1/web/files/jnaat-files/journals/1/articles/js/svgxuse.js"></script>
</body>
</html>