On considère le problème aux limites suivant :
$
\begin{cases}
-u''(x)+c(x)u(x)=f(x), & 0<x<1,\
u(0)=u(1)=0.
\end{cases}
Les fonctions $f$ et $c$ sont continues et données sur $[0,1]$. On suppose $c(x)\geq0$ pour $x\in[0,1]$.
On subdivise lΓÇÖintervalle $[0,1]$ en $N+1$ sous-intervalles $[x_i,x_{i+1}]$ avec $x_i=ih$, pour $i=0,\ldots,N+1$, o├╣ $h=\frac{1}{N+1}$. On note $u_i$ lΓÇÖapproximation de la solution en $x_i$, $f(x_i)=f_i$ et $c(x_i)=c_i$.
Le schéma d’approximation vérifie
$
\begin{cases}
\frac{-u_{i+1}+2u_i-u_{i-1}}{h^2}+c_iu_i=f_i, & i=1,\ldots,N,\\
u_0=0,\quad u_{N+1}=0.
\end{cases}
- Écrire le problème approché sous la forme d’un système linéaire Au=b, où
$
u=(u_1,u_2,\ldots,u_N)^T,\qquad b=h^2(f_1,f_2,\ldots,f_N)^T.
2. Montrer que
$
v^TAv=v_1^2+v_N^2+h^2\sum_{i=1}^{N}c_iv_i^2+\sum_{i=1}^{N-1}(v_{i+1}-v_i)^2,
pour tout v=(v1,v2,…,vN)∈RN.
-
En déduire que le système linéaire admet une solution unique.
-
Pour i=1,2,…,N, on pose
$
R_i=\frac{1}{h^2}\left(2u(x_i)-u(x_{i-1})-u(x_{i+1})\right)+c(x_i)u(x_i)-f(x_i).
QuΓÇÖappelle-t-on $R_i$ ? Montrer que
$
\lvert R_i\rvert\leq\frac{h^2}{12}\sup_{[0,1]}\lvert u^{(4)}\rvert.
``$