Université Sultan Moulay Slimane FST Béni Mellal Département de Physique
Année universitaire : 2023-2024
Filière d'Ingénieur d'État 1ère Année
Génie Physique - Matériaux et Énergie
(Option : Energie)
COMPTE RENDU :
TRAVAUX PRATIQUES
Module : Informatique et Analyse
Numérique
Préparé par : Sous la direction de :
AIT IBOUH Rida MR. Amin BOUKDIR
LAHYOULI Bilal
LADIB Youssef
TP N￮1
Exercice 1:
clc;
clear all;
U=[1 -4 -1 1;1 7 1 -2;1 4 -1 2;3 -10 -2 5] % pour la matrice U
B1=[U(2,2:3);U(3,1:2)] % pour déterminer la matrice B1
B2=U(2:4,3:4) % pour déterminer la matrice B1
U1=triu(U) % pour déterminer la matrice B1
I=eye(4)
A=I-U*U % pour déterminer la matrice A
inv(A) % pour calculer l'inverse de la matrice A
det(A) % pour calculer la déterminant de la matrice A
rank(A) % pour calculer le rang de la matrice A
eig(A) % pour calculer les valeurs propres de la matrice A
eig(A) % pour calculer les vecteurs propres de la matrice A
Résultat 1 :
1
Exercice 2 :
clc;
clear all;
x=[0:0.1:1]
Résultat2 :
F1=sin(x).*sqrt(x+1);
subplot(3,1,1)
plot(x,F1,'r')
xlabel('x')
ylabel('F1')
title('f1(x)')
F2=exp(x)+x.^2;
subplot(3,1,2)
plot(x,F2,'g')
xlabel('x')
ylabel('F2')
title('f2(x)')
hold on
F3=1./(1+abs(x));
subplot(3,1,3)
plot(x,F3,'b')
xlabel('x')
ylabel('F3')
title('f3(x)')
2
Exercice 3 :
a=input('entrer la valeur de a=');
b=input('entrer la valeur de b=');
c=input('entrer la valeur de c=');
Résultat3 :
if a>=b && a>=c
max_val = a;
elseif b>=a && b>=c
max_val = b;
else
max_val = c;
end
disp(max_val);
Exercice 4 :
a.
clc;
clear all;
S=0;
i=1;
while i<=100
S=S+2*sin(i);
i=i+1;
end;
disp(S);
3
b.
clc;
clear all;
n=59
while n<=99
n=n+1
if n/2==round (n/2)
'n est paire'
else
'n est impaire'
end
end
4
TP N￮2
Méthode de Horner:
1)Évaluation à l'aide de la forme initiale de Pn:
2)Évaluation par la forme de Horner:
5
Matrice de Vandermonde:
6
TP N￮3
1)Rectangles simple :
7
2)Rectangles composée :
8
Méthode des trapèzes :
1)Trapèze simple :
2)Trapèze composée :
9
Méthode de Simpson :
1)Simpson simple :
2)Simpson composée :
10
