11.08.2013 Views

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Numerical Integration 63<br />

(3.10.5) ˜ω (n)<br />

j<br />

= Γ(n + α + 1)<br />

4n 2 n!<br />

4n + η (n)<br />

j<br />

ˆL (α)<br />

n (η (n)<br />

j )<br />

where we noted that, for n ≥ 2<br />

(3.10.6)<br />

Again we have<br />

(3.10.7)<br />

where ˆp := pS (α)<br />

n .<br />

d<br />

<br />

S<br />

dx<br />

(α)<br />

n−1 (x)<br />

−1 <br />

:= ˜w (n)<br />

j<br />

<br />

S (α)<br />

n (η (n)<br />

−2 j )<br />

d<br />

dx ˆ L (α)<br />

n−1 (η(n) j ) + ˆ L (α)<br />

n−1 (η(n)<br />

n−1 <br />

j )<br />

m=1<br />

=<br />

+∞<br />

p<br />

0<br />

2 (x) x α e −x dx =<br />

<br />

S (α)<br />

n−1 (x)<br />

−1 n−1 <br />

n−1 <br />

j=0<br />

m=1<br />

1<br />

4m + η (n)<br />

j<br />

1<br />

, x ∈ [0,+∞[.<br />

4m + x<br />

ˆp 2 (η (n)<br />

j ) ˜ω (n)<br />

j , ∀p ∈ Pn−1,<br />

The same procedure applied to the Herm<strong>it</strong>e weights in (3.4.9), yields<br />

(3.10.8) ω (n)<br />

j<br />

=<br />

√ π n!<br />

2 n−1 [(n/2)!] 2<br />

(3.10.9) ω (n)<br />

j<br />

=<br />

:= w (n)<br />

j<br />

<br />

S (−1/2)<br />

n/2 ([ξ (n)<br />

j ] 2 −2 )<br />

<br />

ˆH ′<br />

n(ξ (n)<br />

−2 j ) , 1 ≤ j ≤ n, if n is even,<br />

:= w (n)<br />

j<br />

√ π n!<br />

2 n−1 [((n − 1)/2)!] 2<br />

<br />

S (1/2)<br />

(n−1)/2 ([ξ(n) j ] 2 −2 )<br />

<br />

ˆH ′<br />

n(ξ (n)<br />

−2 j ) , 1 ≤ j ≤ n, if n is odd.<br />

−1<br />

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

Here the ξ (n)<br />

j ’s are the Herm<strong>it</strong>e zeroes. Relations (1.7.9), (1.7.10), (1.7.11) and (1.7.12)<br />

have been taken into account.<br />

,

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

Saved successfully!

Ooh no, something went wrong!