CCINP Maths 2 PSI 2001

SESSION 2001 \ PSIOOG

A

CONCOURS (OMMUNS POlYTECHNIOUES

ÉPREUVE SPÉCIFIQUE - FILIÈRE PSI

MATHÉMATIQUES 2

DURÉE : 4 heures

Les calculatrices programmables et alphanumériques sont autorisées, sous 
réserve des conditions
définies dans la circulaire n° 99-186 du 16.11.99 .. BOEN n°42 du 25.11.99.

L'objet de ce problème est de définir un algorithme de calcul approché d'une 
intégrale, utilisant la
méthode des trapèzes.

Dans la première partie, on étudie le procédé d'extrapolation de Richardson.
Dans une deuxième partie, on établit la formule d'Euler -- Mac Laurin.

La troisième partie utilise les deux premières parties pour définir la méthode 
de Romberg, qui est

une méthode d'intégration basée sur l'accélération de la convergence à partir 
de la méthode des
trapèzes.

La deuxième partie est indépendante de la première partie.

De nombreuses questions de ce problème sont simples ; le candidat s'attachera à 
les résoudre avec
soin et complètement.

On note R l'ensemble des nombres réels, N l'ensemble des entiers naturels et N* 
l'ensemble des
entiers naturels non nuls.

Etant donné un intervalle I de R, on note C °°(I, R) l'ensemble des fonctions 
définies sur I à valeurs
dans R, indéfiniment dérivables.

Etant donné un entier se N et une fonction  O,

@@

tS

qui signifie que le quotient est borné lorsque t ---> 0, t # 0.
t,, étant le terme général d'une suite qui ne s'annule pas et qui tend vers 0 
10rsque n -----> +oo, on note
u,.=O(tn) lorsque n --> +oo, le terme général d'une suite telle que le quotient 
y--£'-- est borné lorsque

ln
ll "9 +00.

Tournez la page S.V.P.

PREMIÈRE PARTIE

Procédé d'extrapolation de Richardson

On désigne par A une fonction définie sur R à valeurs dans R, et on suppose que 
A admet un
développement limité à tout ordre au voisinage de 0.

On note A( t ) : ao + alt + -- --- + aktk + O(tk+l) son développement limité à 
l'ordre k au voisinage de
0, les coefficients a,, étant des réels.

1.1.1 Etant donné un réel p non nul et un entier se N , on suppose que (p est 
une fonction qui
vérifie (p(t) : O((pt)') lorsque t ----> 0 .Montrer que (p(t) : O(t') lorsque t 
----> 0.

1.1.2 Pour k EUR N * , on suppose que ç0(t) : O(t") lorsque t ---> 0. 
Déterminer la limite lorsque 1 ---9 0,

t # O , du quotient (;)/((i) .

1.2.1. Montrer que A( t) admet une limite lorsque t ----> 0 et déterminer cette 
limite.

Soit r un réel vérifiant r > 1. On définit la suite des fonctions A,, par :

n
, , r A __ t --- A __ rt
pour t réel, A0(Ï) =A(t) plus, pour t reel et n E N*, A,, 
(t)=---------'--'--l%)----Ë--LÇ--).

r ---1

1.2.2. Montrer qu'il existe un réel a... , que l'on déterminera, tel que le 
développement limité de A,
à l'ordre k au voisinage de 0 soit A1 (t): a() + al'2t2 + - . - + O(t""1 ).
1.2.3. En déduire qu'il existe un réel a,,_n+1_ que l'on ne demande pas de 
déterminer, tel que le
développement limité de A,, à l'ordre k au voisinage de 0 soit A,, (t') = a() + 
a,m+lt"+1 + - - - + 0(t""1 ).

1.2.4. Soit to un réel non nul fixé. Montrer que la suite de terme général 
A(r--"' to) converge vers ao
lorsque m -----> + oo.

Dans la suite de la première partie, on suppose que pour tout to # 0 fixé et r 
>1 fixé, on sait calculer
les premiers termes A(to), A(r "'to), ..., A(r--'" to) de la suite.

Le procédé de Richardson consiste à extrapoler ces valeurs pour obtenir, grâce 
à un procédé
d'accélération de convergence, la valeur de ao_

Pour p E N, on note Ap,o : A0 (r"p to) puis, pour q entier vérifiant 1 5 q S p, 
on note Ap_q : Aq(r'p to).

1.3.1. Justifier l'égalité A,... : ao + O(r"p ) lorsque p ----> +00 .

1.3.2. Déterminer un entier naturel a(p,q)>0, que l'on expliciterà, tel que 
AP,q : aO+O(r" °"""")
lorsque p ---> +oo.

1.3.3. Pour p 2 l, justifier l'égalité : Ap,l 1
r ..

1.3.4. Pour 1 5 q 5 p, justifier l'égalité :
r"A _] --- A

prq "'1, --1
qu : ' p q :AP (14 + (AP (H _A
' rq ----1 r" ----1

Dans la pratique on range les valeurs AP,q pour 0 5 q 5 p .<... m dans un tableau triangulaire : Ao,o A1,0 Al,l A2,0 A2,1 A2,2 A A m,2 m ,m Am,0 A m,l 1.4. Déterminer la plus petite valeur et la plus grande valeur de oc (p,q) pour 0 5 q S p 5 m . Lorsque m --> +00 , de laquelle des valeurs Ap,q du tableau peut--on attendre 
la meilleure

approximation de ao (on pourra utiliser 1.1.2 pour justifier la réponse) ?

On écrira cette valeur sous la forme ao + O(r" 0(m)) lorsque m ---> +00 et on 
précisera la valeur de
l'entier naturel a(m) > O .

On considère une fonction g E C °°(R, R) et on note g(oc+h) : co + c;h + ...+ 
cu h2k + O(h2k+l) son

développement limité à l'ordre 2k au voisinage de oc.

1.5.1. Exprimer les coefficients c[) pour 0 S p _<_ 2k, en fonction de g et de ses dérivées successives. - Pour h : 0, on note G(h) = ËË--Ülz'lË--(Y--ZËÎ. 1.5.2. Montrer que la fonction G est paire. Montrer que G se prolonge par continuité en 0 par une valeur que l'on déterminera. On note G la fonction G prolongée en 0 par cette valeur. N 1.5.3. Exprimer à l'aide des coefficients cp le développement limité de G à l'ordre 2k ----1 au voisinage de 0. Pour l réel positif, on note A(t)=G(JÏ). h 1.6.1. On choisit h > 0 et on considère la suite de valeurs G(h), G(--2--), 
..., G(--£-fi--].
Déterminer un réel to > 0 et un réel r >1 tels que cette suite de valeurs soit 
A(to), A(r"1 to),
A(fm IQ).

...,

Tournez la page S.V.P.

On utilise les notations des questions précédentes avec A0( t) : A(t), A,,,g : 
Ao(r' %) puis A..., ,
pour les valeurs r et to déterminées dans 1.6.1.

1.6.2. Quelle est la limite EUR de A,... lorsque p----> +oo '? On exprimera EUR 
à l'aide de la fonction g
et de oc.

Dans ce qui suit, on prend g(x) : Ën x, ou = 3 et h = 0,8.

I.7.1. Calculer les valeurs A...) pour 0 .<. p S 3. Donner le tableau des valeurs A..., pour 0 5 q 5 p 5 3. I.7.2. Quelle est la valeur exacte de EUR? Parmi les valeurs Am trouvées, quelle est la meilleure approximation de [? DEUXIÈME PARTIE Formule d'Euler --Mac Laurin Pour p EUR N, on définit la suite B,) de polynômes par : (i) Pour tout t E R, Bo (t) = 1 (ii) Pour pE N* et te R, B;,(r)= po_l(t) et Ipr(t)dt=0, () et on note b,, : B,,(O). II.1.1. Déterminer les polynômes 31, 32, B3. 11.1.2. Pour 0 .<. p S 3, calculer bp et comparer b,, à B,) (1). II.1.3. Montrer que pour p .>. 2, on a b,) : Bp (1).

11.2.1. Pour p e N et : E R, on définit Ëp(r)= (-1)P BP (1 --r).
Montrer que la suite de polynômes ËP vérifie les relations (i) et (ii). En 
déduire que ËP : B p .
Il.2.2. Montrer que pour p EUR N*, on a 172,0": 0.

Soit fe C °°([0,1], R) ; on note f"" la dérivée d'ordre p de la fonction f.
11.3.1. Montrerl'égalité [ 1f(t)dr= [ 1B.,(r) f(t)dt=â--(f(0)+ f(1))--j IBI(t) 
f'(r)æ.
0 0 ' 0

II.3.2. Pour n 2 2, montrer l'égalité

â--(f(0)+ f(1))= J;f(t) .. + Ë<-- 1)? " _L 2 P-' (f """(l)-- f'""'(0))+ (-1)flflj(') En?) f...(t) dt. ." II.3.3. En déduire que pour n-- -- 2kbzp on a l'égalité ... à<:@»:@> =)j):@ .:.. @P(f<@--1111f@--...) J@@fa( p=1 1(2p)! 0 (2k)! Pour t E R, on note E(t) la partie entière de t. Pour p E N, on définit la fonction Dp par: pour t E R, D,,(t) : B,,(t-E(t)). II.4.1. Montrer que D,) est une fonction périodique de période 1. Montrer que D[) est une fonction de classe C °°par morceaux sur R. Dans la suite la fonction f appartient à C °°([O,N], R) où N E Navec NZ 2. Pour q entier vérifiant 1 5 q S N, on définit les fonctions fq de [0,1] dans R par fq (Ï)=f(ï+CI--1) II.4.2. Montrer que les fonctions f() appartiennent à C oe([O,l], R) et qu'elles vérifient les égalités : pour m E N et q entier tel que 2 .<.. q 5 N, f:""'(0)=f'm)(0), :°"'(0 0=) f.'£'9(l) :S"')(l)=f(m)(N). 11.4.3. En appliquant (1) aux fonctions fq, en déduire la formule d'Euler --Mac Laurin sur [O,N] : (2) â-f(O)+Êf(q)+ --â--f(N)= JZf(t)d:+Ê2 b2p !----)((f(2p 1 (N)--fl @ 1)(0))_ JNP_2ÀQ'f(uW p=1(2) P) ("')- TROISIÈME PARTIE Méthode de Romberg Dans cette partie on note [a, b] un intervalle de R et f une fonction de C °°([a, b], R). _ --1 Etant donné N E N* et h=b a, on note Tf(h)=1{--â--f@)+î f(a+qh)%f(b)} si N 2 , q=1 N Tf (h) : hË-- f (a)+ % f (b)] si N = 1 , la valeur approchée de l'intégrale Jb f (t) dt obtenue par la méthode des trapèzes pour le pas h. III. 1. On suppose N> 2. En appliquant la formule (2) à la fonction g(t)= f(a + 
ht) définie sur
[0 ,N], montrer la formule

bD2k("Î) ...
(%)! f(( (tdt)

(3) T(h )=î(:)jf d(::+pË_h2p(2; )!--(f(2p 1((b)--f(2p'a )) th)

Tournez la page S.V.P.

III.2. Montrer que la formule (3) peut s'écrire : (4) Tf(h) =)bf(t) 
dt+Îdph2"+0(hu) où les dp
a p=1

désignent des nombres réels.

Pour : > 0, on définit A(t) = T, (J? )

111.3.1. Déterminer lim A(t).
t-->O

III.3.2. On prend N = 1 et donc h = b -- a , et on calcule la suite de valeurs

T,(h), T,%), Tf{--2b--).

Déterminer un réel to >O et un réel r > 1 tels que cette suite de valeurs soit

A(t ), A(r"to),..., A(r--mm).

A(t), Ap,0 =Ao(r"Ptg) puis, AM, zo
b--a
2P '

On utilise les notations de la première partie avec Ao(t)

et r étant les valeurs trouvées en III.3.Z. On note h p :

III.4.l. Exprimer Ap,o et Ap-l,0 à l'aide de Tf et de hp.

III.4.2. Pour p _>_ 1, on définit A;... =hp2 f (a+(2q+l) hp), la somme étant 
étendue aux entiers q

telsque a