On désigne par {a,b:[0,1]\to\mathbb R} deux applications continues, avec {a(t)\geq0}. Soient {\lambda,\mu\in\mathbb R}. On considère le problème aux limites {f''(t)=a(t)f(t)+b(t),\quad f(0)=\lambda,\quad f(1)=\mu}On suppose qu’il admet une solution {f}, de classe {\mathcal C^4} sur {[0,1]}. On note {M_2(f)=\sup\limits_{[0,1]}|f''|} et {M_4(f)=\sup\limits_{[0,1]}|f^{(4)}|}.
L’objectif est d’étudier une méthode d’approximation de {f}. Les parties I et II sont indépendantes.
Première partie
Soit {g\in\mathcal C^2([\alpha,\beta],\mathbb R)}, avec {\alpha\lt\beta}, {g(\alpha)=g(\beta)=0}. On pose {M_2(g)=\sup\limits_{[\alpha,\beta]}|g''|}.
I.1. (a) Montrer que {|g(t)|\leq\dfrac{(t-\alpha)(\beta-t)}2M_2(g)} pour tout {t\in[\alpha,\beta]}.
On fixe {n\in\mathbb N}, on pose {h=1/(n+1)} et {t_k=kh} pour {0\leq k\leq n+1}. On dispose d’approximations {u_k} de {f(t_k)} et on note {\delta_n=\max_{0\leq k\leq n+1}|u_k-f(t_k)|}On définit {\varphi_n} par {\varphi_n(t_k)=u_k}, avec {\varphi_n} affine sur chaque {[t_k,t_{k+1}]}.
I.3. Montrer que {\sup\limits_{[0,1]}|f-\varphi_n|\leq M_2(f)/(8(n+1)^2)+\delta_n}.
On fixe {n\geq2}, des réels {\alpha_1,\ldots,\alpha_n\geq2} et {\beta_1,\ldots,\beta_n}. On considère le système tridiagonal {(S)} : {\begin{cases}\alpha_1u_1-u_2=\beta_1\\[4pt]-u_{k-1}+\alpha_ku_k-u_{k+1}=\beta_k\quad(2\leq k\leq n-1)\\[4pt]-u_{n-1}+\alpha_nu_n=\beta_n\end{cases}}
II.2. On définit {v_1=\beta_1\gamma_1} et {v_{k+1}=(\beta_{k+1}+v_k)\gamma_{k+1}}. Montrer que {(S)} admet une solution unique, donnée par {u_n=v_n,\quad u_k=v_k+\gamma_ku_{k+1}\quad(1\leq k\lt n)\quad(\Sigma)}
III.2. Avec les notations de I.3, montrer qu’il existe des réels {\varepsilon_k}, {1\leq k\leq n}, tels que {\begin{array}{l}-f(t_{k-1})+(2+h^2a(t_k))f(t_k)-f(t_{k+1})\\[4pt]=-h^2b(t_k)+\varepsilon_k,\quad |\varepsilon_k|\leq h^4M_4(f)/12\end{array}}
On connaît {u_0=\lambda=f(t_0)} et {u_{n+1}=\mu=f(t_{n+1})}. Pour approcher les valeurs intérieures, on résout {-u_{k-1}+(2+h^2a(t_k))u_k-u_{k+1}=-h^2b(t_k)}pour {1\leq k\leq n}.