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.

Derivative Matrices 145<br />

The matrix is given by squaring the equivalent of the matrix in (7.4.11), for a general<br />

n ≥ 2. In fact, (7.4.23) can be decoupled, w<strong>it</strong>h the subst<strong>it</strong>ution ˆp := p ′′ , into two<br />

consecutive second order systems.<br />

We finally conclude w<strong>it</strong>h some examples relative to unbounded domains. Let us<br />

consider first the case I =]0,+∞[. Here, problems are formulated in the space of<br />

Laguerre functions (see sections 6.7 and 9.5). For Q ∈ Sn−2, n ≥ 2, σ ∈ R and<br />

µ > 0, we are concerned w<strong>it</strong>h finding the solution P ∈ Sn−1 of the problem<br />

(7.4.25)<br />

⎧<br />

⎨<br />

⎩<br />

−P ′′ (η (n)<br />

i ) + µP(η (n)<br />

i ) = Q(η (n)<br />

i ), 1 ≤ i ≤ n − 1,<br />

P(η (n)<br />

0 ) = σ,<br />

where the η (n)<br />

i ’s are the Laguerre Gauss-Radau nodes (see section 3.6). By setting<br />

q(x) := Q(x)e x , p(x) := P(x)e x , ∀x ∈ [0,+∞[, we get an equivalent formulation in<br />

the space of polynomials<br />

(7.4.26)<br />

⎧<br />

⎨<br />

⎩<br />

−p ′′ (η (n)<br />

i ) + 2p ′ (η (n)<br />

i ) + (µ − 1)p(η (n)<br />

i ) = q(η (n)<br />

i ), 1 ≤ i ≤ n − 1,<br />

p(η (n)<br />

0 ) = σ.<br />

This leads to a n×n linear system. The entries of the corresponding matrix are obtained<br />

in the usual way. Similar considerations hold for other types of boundary cond<strong>it</strong>ions. We<br />

remark that we only have one boundary point for a second-order equation. Nevertheless,<br />

problem (7.4.25) adm<strong>it</strong>s a unique solution. Actually, another cond<strong>it</strong>ion is implic<strong>it</strong>ly<br />

assumed by observing that P decays to zero at infin<strong>it</strong>y exponentially fast. More<br />

details are given in section 9.5.<br />

For the Herm<strong>it</strong>e case there are no boundary points. Let us discuss an example.<br />

The generating problem will be given in section 10.2. Let Q ∈ Sn−1, n ≥ 1, be a<br />

Herm<strong>it</strong>e function and µ > 0. We want to find P ∈ Sn−1 such that<br />

(7.4.27) −P ′′ (ξ (n)<br />

i ) − 2ξ (n)<br />

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

i ) + µP(ξ (n)<br />

i ) = Q(ξ (n)<br />

i ), 1 ≤ i ≤ n,

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

Saved successfully!

Ooh no, something went wrong!