11.08.2013 Views

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Numerical Integration 49<br />

(3.4.2) w (n)<br />

j<br />

(2n + α + β) Γ(n + α) Γ(n + β)<br />

= 2α+β<br />

n! Γ(n + α + β + 1)<br />

×<br />

<br />

P (α,β)<br />

n−1 (ξ(n) j ) d<br />

dx<br />

(α,β)<br />

P n (ξ (n)<br />

−1 j ) , 1 ≤ j ≤ n.<br />

For the proof we argue as follows. By equating the leading coefficients, (1.3.5) yields<br />

(3.4.3)<br />

un(x)<br />

x − ξ (n)<br />

j<br />

= (2n + α + β)(2n + α + β − 1)<br />

2n (n + α + β)<br />

un−1(x) + qj(x), 1 ≤ j ≤ n,<br />

where qj ∈ Pn−2, 1 ≤ j ≤ n, and un = P (α,β)<br />

n . Therefore, recalling (3.4.1) and (3.2.4),<br />

(3.4.4) w (n)<br />

j =<br />

=<br />

=<br />

1<br />

un−1(ξ (n)<br />

j )<br />

1<br />

1<br />

−1<br />

(un−1u ′ n)(ξ (n)<br />

j )<br />

1<br />

un−1(ξ (n)<br />

j )<br />

n<br />

i=1<br />

l (n)<br />

j un−1w dx =<br />

(l (n) (n)<br />

j un−1)(ξ i ) w (n)<br />

i<br />

1<br />

(un−1u ′ n)(ξ (n)<br />

j )<br />

(2n + α + β)(2n + α + β − 1)<br />

2n (n + α + β)<br />

Finally, (3.4.2) is obtained w<strong>it</strong>h the help of (2.2.10).<br />

1<br />

−1<br />

un<br />

x − ξ (n)<br />

j<br />

un−1w dx<br />

un−1 2 w, 1 ≤ j ≤ n.<br />

Legendre case - To obtain an expression for the weights, we simply set α = β = 0 in<br />

(3.4.2). This gives<br />

(3.4.5) w (n)<br />

j<br />

2<br />

=<br />

n<br />

<br />

Pn−1(ξ (n)<br />

j )P ′ n(ξ (n)<br />

−1 j ) , 1 ≤ j ≤ n.<br />

Chebyshev case - This is the most remarkable case. We recall (1.5.1), (3.1.16) and the<br />

relation Tn−1(ξ (n)<br />

j ) = (−1) j+n<br />

(3.4.2), one obtains<br />

(3.4.6) w (n)<br />

j<br />

<br />

1 − [ξ (n)<br />

j ] 2 , 1 ≤ j ≤ n. By taking α = β = − 1<br />

2 in<br />

π<br />

= , 1 ≤ j ≤ n.<br />

n

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!