EX3
EX3
EX3
EE : 2 HEURES
DOCUMENTS NON AUTORIS
ES
NBRE PAGES :2
EXAMEN
Exercice 1 : (6 pts)
Soit la fonction denie par : f(t) =
1
1 + t
2
. On admet que
_
+
0
f(t)dt =
2
.
Donner les instructions Maple permettant de :
1. Denir la fonction f.
2. Verier la partite de f, en achant un message texte adequat.
3. Calculer le developpement limite de f au voisinage de 0 dordre 10.
4. Denir la fonction denie par g(x) =
_
x
0
f(t)dt.
5. Calculer la limite en + de g.
6. Denir la suite denie par : S
n,a,b
=
ba
n
i=n1
i=0
f(a + i
ba
n
).
7. Pour n assez grand, on admet que S
n
converge vers
_
b
a
f(t)dt.
Ecrire une procedure en Maple Ap-
proxIntegrale prenant en entree a, b et epsilon et renvoyant une valeur approchee de
_
b
a
f(t)dt avec une
precision de epsilon pr`es. (|S
n+1
S
n
| epsilon).
8. On denit la suite U
k
= S
n,0,k
. On admet la suite U
k
converge vers
_
+
0
f(t)dt =
2
. Denir une
procedure en Maple ApproxPI dargument un reel epsilon et renvoyant une valeur approchee de avec
une precision de epsilon pr`es.
Exercice 2 : (4 pts)
Une matrice carree A est dite nilpotente sil existe un entier k> 1 tel que A
k
= 0. Soit la matrice suivante :
A =
_
_
0 2 1
0 0 1
0 m 0
_
_
Donner les instructions Maple permettant de :
1. Denir A.
2. Calculer le determinant DT , linverse IV, et le rang RG de A.
3. Recuperer dans P le polynome caracteristique de A en x.
4.
Evaluer P pour x=3.
5. Donner le degre de P.
6. Trouver le reel m tel que P = x
3
.
7. Pour la valeur trouvee de m, evaluer A.
8.
Ecrire une procedure en Maple IndiceNilpotence qui retourne lindice de niloptence de A( On suppose
A est nilpotente pour la valeur trouvee de m) .
1
Probl`eme : (10 pts)
On se propose de resoudre un syst`eme lin`eaire par la Methode de Jacobi . Soit le syst`eme Ax = b o` u
A =
_
_
a
11
. . a
1n
a
21
. . a
2n
. . . .
a
n1
. . a
nn
_
_
, b =
_
_
b
1
b
2
.
b
n
_
_
, x =
_
_
x
1
x
2
.
x
n
_
_
A partir de ce syst`eme , la matrice de Jacobi J et le vecteur constant c sont denis par les formules suivantes :
_
_
J
ii
= 0
J
ij
=
a
ij
a
ii
1 i n, i = j
c
i
=
b
i
a
ii
La methode de Jacobi permet de calculer iterativement les composantes des solutions directement par la
formule de recurrence suivante :
_
x
(k+1)
= Jx
(k)
+ c
x
0
= 0
x
(k)
represente le vecteur x `a literation k et 0 cest le vecteur nul c-a-d x
0
i
= 0, pour tout 1 i n. Pour
que la methode de Jacobi converge vers un vecteur solution S il faut que A est `a diagonale strictement
dominante.
A est `a diagonale strictement dominante si |a
ii
| >
j=n
j=1,j=i
|a
ij
| pour tout 1 i n.
On denit le produit dune matrice M par un vecteur v par le vecteur resultat x = Mv o` u x
i
=
j=n
j=1
M
ij
v
j
avec M
ij
le coecient de M dindices i,j et v
j
le coecient dindice j du vecteur v. On denit la somme de
deux vecteurs u et v par le vectuer resultat x o` u x
i
= u
i
+v
i
La dierence de u et v est deni par x = uv
o` u x o` u x
i
= u
i
v
i
. (1 i n). On denit aussi la norme dun vecteur v comme suit :
norme(v) =
_
i=n
i=1
v
2
i
Le nombre diterations necessaires pour obtenir une precision doit satisfaire norme(x
(k+1)
x
(k)
) .
Dans la suite on suppose avoir eectue les declarations suivantes :
CONSTANTE NMAX=100
TYPE MATRICE =TABLEAU[1..NMAX,1..NMAX] DE REEL
TYPE VECTEUR =TABLEAU[1..NMAX] DE REEL
1.
Ecrire une procedure algorithmique SAISIR SYSTEME permettant saisir un entier n entre 1 et
Nmax , remplir une matrice A carree n n dont les elements de la diagonale sont dierents de 0 et
saisir un vecteur b de taille n.
2.
Ecrire une procedure algorithmique MATRTCE JACOBI permettant construire la matrice J et le
vecteur c `a partir de A et b.
3.
Ecrire une procedure algorithmique PRODUIT permettant de calculer le produit (Mv) dune matrice
M par un vectuer v de meme taille n.
4.
Ecrire une procedure algorithmique SOMME permettant de calculer la somme (u+v) de deux vecteurs
u et v de meme taille n.
5.
Ecrire une procedure algorithmique DIFFERENCE permettant de calculer la dierence (u v) de
deux vectuers u et v de meme taille n.
6. Ecrire une fonction algorithmique NORME retournant la norme dun vecteur v de taille n.
7. Ecrire une fonction algorithmique DIAGONALE DOMINANTE qui retourne vrai si une matrice
A de taille n est `a diagonale strictement dominante et faux sinon .
8. Ecrire une procedure algorithmique JACOBI darguments une matrice A carree n n et un vecteur
b de taille n , et permettant de calculer dans un vecteur S, une solution approchee du syst`eme AX=b
de epsilon pr`es
2