TD Enonces

L3-Mag1 Phys. fond.
16-17 Informatique 4 octobre 2016
L3 et Magistère 1ère année, Physique fondamentale, année 2016-2017
Travaux dirigés de langage C appliqué à la

physique
— Textes d’exercices de TD (n◦ 1 à 8)
1. Types entier et réel, priorité des opérations

2. Limitations des valeurs entières et réelles en grandeur et en précision
Diagnostics en cas d’opérations illégales
3. Tests et boucles
4. Fonctions
5. Pointeurs
6. Tableaux dynamiques
7. Générateur aléatoire
8. Équations différentielles
— L’éditeur emacs
— Commandes Linux de base
1
L3-Mag1 Phys. fond., TD info C 16-17 n◦ 1 4 octobre 2016
1. Types entier et réel, priorité des opérations
Rappels :
— Ouvrir le site du cours avec les énoncés des TD, les notes de cours, et plein d’autres informations utiles :
site
— Il sera utile de créer un sous-répertoire nommé TD1 (avec la commande mkdir TD1) et s’y placer (avec la commande
cd TD1) pour faire ce TD.
— Ouvrir l’éditeur emacs :
emacs prog1.cpp & ou
em prog1.cpp
— En C il y a distinction entre minuscules et majuscules.
— Pour afficher à l’écran de la façon la plus simple, placer en tête du fichier contenant le programme :
#include<iostream>
using namespace std;
et utiliser, par exemple dans le cas d’une variable k :
cout << k << endl;
— Pour compiler et exécuter en une seule commande :
ccc nom_de_fichier.cpp
(Alternativement, il faut d’abord compiler : g++ -lm nom_de_fichier.cpp puis exécuter : ./a.out)
Prévoir le résultat des quatre programmes C dont le texte est donné ci-dessous, puis écrire chacun de ces programmes
dans un fichier propre, le faire exécuter et comparer au résultat prévu. Seul le programme 1 doit être tapé intégralement.
A partir du programme 2, il suffit d’utiliser la copie de fichiers et de ne taper que ce qui doit être modifié.
Types entier et réel

Ici on montre les règles qu’il faut suivre pour manier les constantes ou les variables entières et réelles dans les expressions
arithmétiques et algébriques.
Programme 1
#include<iostream>
int main() {
cout << 2/3 << " " << 2./3 << " " << 2/3. << " " << 2./3. << endl;
return 0;
}
Programme 2
#include<iostream>
int main() {
int i = 2, j = 3, k ;
double x = 2., q ;
cout << i/j << " " << x/j << endl;
k = x/j ; q = x/j ;
cout << k << " " << q << endl;
return 0;
}
Comment obtenir 0.6666667 pour les quatre quantités calculées ?
2
Règles de priorité des opérations

Programme 3
#include<iostream>
int main() {
cout << (2.*3.)/4. << " " << 2.*(3./4.) << " " << 2.*3./4. << endl;
cout << (2./3.)*4. << " " << 2./(3.*4.) << " " << 2./3./4. << endl;
cout << (2.+3.)/4. << " " << 2.+3./4. << endl;
return 0;
}
Mélange de types
Programme 4
#include<iostream>
int main() {
cout << (2.*3)/4 << " " << 2.*(3/4) << " " << 2.*3/4 << endl;
return 0;
}
Usage des parenthèses et des variables intermédiaires

∗ ∗ ∗ Paquet d’onde ∗ ∗ ∗
La densité de probabilité d’un paquet d’ondes gaussien libre en mécanique quantique est :
r " 2 #
2 1 2 x − ~k
m t
0
ρ(x) = exp − 2 t2
a2 1 + 4~
q
π a 1 + 4~2 t2 m2 a4
m2 a4
Prendre telles quelles les valeurs numériques suivantes qui sont exprimées dans un système d’unités (fermi, Mev) approprié
au cas d’une particule telle que l’électron :
a = 105 fm (paramètre déterminant la largeur du paquet d’ondes à t = 0 avec 1 fm= 1 Fermi= 10−15 mètre)
~ = 197 MeV.fm.c−1 (constante de Planck avec MeV = Méga électron-volt)
t = 105 fm.c−1 (instant considéré, c = vitesse de la lumière)
m = 0.5 MeV.c−2 (masse de la particule)
x = 105 fm (position à laquelle on calcule la densité de probabilité)
k0 = 5.10−4 fm−1 (paramètre déterminant la vitesse ~k0 /m du centre du paquet d’ondes)
Calculer ρ(x) dans les unités indiquées. Ne pas oublier d’utiliser des variables intermédiaires pour simplifier l’expression
finale (comme par exemple u = 2~t/(ma2 )).
Remarques
Il peut être utile de laisser des parenthèses superflues si elles rendent l’expression plus lisible.
Quelle est l’unité de ρ(x) ?
3
2. Limitations des valeurs entières et réelles

en grandeur et en précision
Préambule :
Créer dans le répertoire initial un sous-répertoire nommé TD2 et s’y placer pour faire ce TD.
(Si vous ne l’avez pas fait la semaine dernière : créer d’abord un sous-répertoire nommé TD1 et copier dans
TD1 tous les fichiers déjà écrits pour faire le TD n◦ 1 (se reporter à la partie Linux du polycopié de TD).)
Remarques
Les limitations et diagnostics que l’on va observer ici se retrouvent sur tous les ordinateurs mais avec des
variantes selon le type de machine. Ce qui est étudié ici n’est donc valable, en détail, que pour les ordinateurs
que nous utilisons.
Pour écrire tous les programmes qui suivent utiliser au maximum la copie (de lignes dans un fichier et d’un
fichier dans un autre).
Résumer les résultats produits par les programmes proposés dans le tableau donné à la fin de ce texte de TD.
Limitations en grandeur et précision

Pour les entiers
Les entiers de type int sont écrits sur 4 octets (32 bits), ceux de type long int sur 8 octets (64 bits).
Un bit est réservé au signe, les autres aux chiffres.
Pour les int ne peuvent donc être écrits que les entiers compris entre -231 =-2 147 483 648 et 231 -1=2 147 483 647 inclus.
Pour les long int ne peuvent être écrits que les entiers compris entre -263 =-9 223 372 036 854 775 808 et 263 -1=9 223
372 036 854 775 807 inclus.
Écrire le programme suivant et remarquer le résultat produit.
#include<iostream>
int main() {
int m;
m = 2147483646;
cout << m << " " << m+1 << " " << m+2 << " " << m+3 << endl;
long int n;
n = 9223372036854775806;
cout << n << " " << n+1 << " " << n+2 << " " << n+3 << endl;
return 0;
}
À part cette limitation les opérations arithmétiques entre entiers sont exactes (compte tenu de ce que la division donne
la partie entière du résultat).
Pour les réels
Pour les réels il y a également des limites et de plus, un nombre n’est en général écrit que de façon approchée, avec un
nombre de chiffres significatifs exacts limité. Lors d’une opération algébrique entre réels le résultat sera lui-même arrondi
ce qui entraîne une nouvelle approximation. Dans certains cas les erreurs peuvent s’accumuler au cours des calculs.
Les réels de type double sont écrits sur 8 octets (= 64 bits).
Pour les réels il existe aussi le type float (sur 4 octets) et le type long double (sur 16 octets) mais nous ne considérerons
4
ici que le type double.

Le tableau suivant résume les limites et la précision des réels de type double pour les ordinateurs utilisés :
Type Nb de Nb de Nb de bits Plus petite Plus grande Nb minimum de

bits bits chiffres valeur absolue valeur absolue chiffres significatifs
signe exposant significatifs représentable représentable exacts en décimale
double 1 11 53 ' 2.225e-308 ' 1.798e308 15

(' 4.941e-324)
Rappel : 1.798e-308 est la notation informatique pour 1.798 · 10−308 . Pour les nombres sub-normaux (limites entre paren-
thèses) il n’y a pas de diagnostic, mais le nombre de chiffres significatifs exacts est moindre.
Pour les valeurs absolues inférieures à la plus petite valeur sub-normale ou supérieures à la plus grande valeur normale il
y a un diagnostic de dépassement.
Contrairement au cas des entiers il y a donc un diagnostic de dépassement pour les réels mais l’exécution se poursuit.
Ecrire et faire exécuter le programme suivant :
#include<iostream>
int main() {
double x = 1.e-308, y = 1.e-10;
cout << "1.e-308 = " << x << endl;
x = x*y;
cout << "1.e-318 = " << x << endl;
x = x*y;
cout << "1.e-328 = " << x << endl;
return 0;
}
puis :
#include<iostream>
int main() {
double x = 1.e308, y = 1.e10;
cout << "1.e308 = " << x << endl;
x = x*y ;
cout << "1.e318 = " << x << endl;
x = x*y ;
cout << "1.e328 = " << x << endl;
return 0;
}
et remplir les lignes correspondantes du tableau récapitulatif donné en fin d’énoncé.
∗ ∗ ∗ Étalement d’un paquet d’onde ∗ ∗ ∗

La largeur du paquet d’ondes étudié au TD n◦ 1 vaut, à l’instant t :
r
a 4~2 t2
∆x(t) = 1+ 2 4
2 m a
a) A t = 0 un électron libre a pour fonction d’onde un paquet d’ondes gaussien avec a = 1Å= 10−10 m.
Calculer ∆x une seconde plus tard. On peut soit chercher les valeurs de ~ et m en unités SI sur internet, soit se reporter
au dernier exercice du TD1 pour trouver les valeurs de ~, m et a en unités (fermi, MeV). Ne pas oublier de convertir la
valeur de t = 1 s dans ce dernier cas.
b) Même question pour un proton (chercher la masse du proton sur internet).
5
c) Même question pour une sphère de diamètre 1 mm, de densité 1, avec a = 1 µm.
Commentaire sur cet exercice, à lire attentivement après avoir fait les questions
2~t
Il faut calculer la quantité sans dimension u = , l’élever au carré et l’ajouter à 1. On peut utiliser indiffé-
ma2
remment les unités microscopiques introduites au TD n◦ 1 ou les unités SI. On obtient les ordres de grandeur
suivants :
u u2
Electron 2.3e+16 5.3e+32
Proton 1.2e+13 1.5e+26
Sphère 4.0e-16 1.6e-31
Il faut ensuite ajouter u2 à 1. Dans le cas de l’électron et du proton l’intervalle entre les deux nombres réels
représentables dans l’ordinateur qui encadrent la valeur de u2 est beaucoup plus grand que 1. Ajouter 1 à
u2 n’aura donc strictement aucun effet, le 1 ne sera pas pris en compte, mais cela n’entraîne qu’une erreur
négligeable. Dans le cas de la sphère, c’est l’inverse : l’intervalle entre 1 et le premier nombre représentable
qui le suit est beaucoup plus grand que u2 , qui n’est donc pas pris en compte. Le résultat du calcul indique
alors un élargissement nul, alors qu’en réalité il est extrêmement petit mais non nul. Pour le calculer il faut
faire un développement limité p
1 + u2 = 1 + u2 /2
justifié par le fait que, dans ce cas, l’erreur commise est inférieure ou égale à u4 /8, et qui indique que la
variation relative de la largeur du paquet d’onde est de l’ordre de 10−31 .

Écrire et faire exécuter le programme suivant :
#include<math.h>
#include<iostream>
int main() {
double x = 3., y = 0., z = -1.;
cout << endl << "Diviser par zéro : " << x/y << endl;
cout << "Prendre la racine carrée d’un réel négatif : " << sqrt(z) << endl << endl;
return 0;
}
et remplir le tableau :
Cas Noter ici le résultat obtenu
Entier trop grand ou trop petit
Réel trop petit en valeur absolue
Réel trop grand en valeur absolue
Diviser par 0
Prendre la racine carrée d’un réel négatif
6
3. Tests et boucles
Remarques
Certains programmes demandés ici pourront être utiles dans des TD ultérieurs en les transformant en fonc-
tions.
Un des buts des exercices suivants est de montrer que des formules équivalentes mathématiquement ne le sont
pas du point de vue de la programmation.
Pour savoir comment changer la notation des réels, obtenir le nombre de chiffres voulu lors de l’affichage d’un
réel avec la fonction cout, écrire et lire dans un fichier, consulter le chapitre du cours : Entrées, sorties,
fichiers, en C++ .
∗ ∗ ∗ Factorielle ∗ ∗ ∗
Calculer n! avec :
1. une boucle for
2. une boucle while
Jusqu’à quelle valeur de n le calcul est-il juste si on travaille en entiers de type int ? Pour répondre à cette question,
trouver une méthode pour tester si la valeur maximale d’une variable du type int est dépassée pendant les calculs.
∗ ∗ ∗ Combinaison ∗ ∗ ∗
Calculer le nombre de combinaisons de k éléments pris parmi n :
n!
Cnk =
k!(n − k)!
en utilisant la formule de récurrence :
n−i+1
Cni = Cni−1
i
Peut-on tirer profit de la propriété :
Cnk = Cnn−k ?
Ce programme permet-il de calculer le nombre de possibilités de choisir 10 élèves parmi 45 ? Et 44 parmi 45 ?
Pourquoi n’applique-t-on pas directement la formule en calculant n!, k!, (n − k)! puis n!/(k!(n − k)!) ?
∗ ∗ ∗ Pluie ∗ ∗ ∗
On considère un évènement qui a la probabilité p d’avoir lieu lors d’une épreuve. On note q = 1 − p la probabilité pour
qu’il n’ait pas lieu. La probabilité b(k) pour qu’il se produise k fois lors de n épreuves est donnée par la loi binomiale :
b(k) = Cnk pk (1 − p)n−k

√
(valeur moyenne de < k >= np et écart-type de k = npq)
Sachant qu’en moyenne il pleut 12 jours sur 30 au mois de novembre, quelle est la probabilité, sur une période de quatre
jours, d’avoir respectivement 0, 1, 2, 3, 4 jours de pluie ? Utiliser une boucle pour afficher ces cinq résultats ; à noter
qu’une seule boucle suffit, on n’a pas besoin d’une boucle séparée pour le calcul des Cnk .
On admet que le temps qu’il fait un jour est indépendant de celui des jours précédents.
∗ ∗ ∗ Suite ∗ ∗ ∗
La limite de la suite :
1 A
u0 = A un+1 = (un + )
2 un
√
est A.
7
a) Calculer ainsi une racine carrée avec une boucle for avec compteur, le nombre d’itérations étant fixé à l’avance à
nmax. On fera afficher les nombres avec le maximum de chiffres significatifs (15) en écrivant, en début de programme,
l’instruction :
cout << setprecision(15);
(setprecision nécessite la directive #include<iomanip>)
b) Ajouter une condition pour que le calcul s’arrête comme au a) ou si :
un+1 − un
| |<ε
un
ε étant fixé à l’avance.
Pour cela on pourra ajouter un break dans le for ou utiliser une boucle do ... while().
Pourrait-on utiliser une boucle while() au lieu de la boucle do ... while() ?
c) Avec les réels double combien faut-il d’itérations pour satisfaire la condition du b) lorsque A = 2 et ε=10−10 ?
d) Au vu des résultats affichés à l’écran, combien faut-il d’itérations, lorsque A = 2, pour que le résultat obtenu soit
identique avec celui donné par la fonction de la bibliothèque mathématique sqrt (inclure la directive #include<math.h>) ?
∗ ∗ ∗ Sinus ∗ ∗ ∗
1) Calculer et afficher sin x de 0◦ à 90◦ de 10◦ en 10◦ à l’aide du développement :
x3 x2n+1
sin x = x − + ... + (−1)n + ...
6 (2n + 1)!
en additionnant tous les termes de la série dont la valeur absolue est supérieure à 10−6 . Pour faire cela, il faut d’abord
déduire une formule de récurrence pour calculer le nouveau terme à partir du terme précédent, vu qu’on ne veut surtout pas
calculer x2n+1 et (2n + 1)! explicitement pour chaque terme à cause des problèmes élucidés dans les exercices précédents.
Comparer au résultat donné par la fonction sin de la bibliothèque mathématique.
2) Pour x = 20◦ par exemple, combien faut-il de termes de la série pour obtenir 6 chiffres significatifs identiques au
résultat de la fonction sin ? Même question pour 12 chiffres significatifs.
∗ ∗ ∗ Test ∗ ∗ ∗
Le test de l’égalité de deux réels peut ne pas être fiable, car la plupart des réels ne peuvent pas être écrits de manière
exacte dans l’ordinateur, le nombre de bits consacré à chacun d’eux étant fini (64 pour les doubles). Quel est le résultat
du programme suivant ?
#include<iostream>
int main() {
int i;
double x;
x = -1.;
for(i = 1; i <= 10; i++)
x = x + 0.1;
if(x != 0.)
cout << "x = " << x << " le test trouve x différent de 0" << endl;
else
cout << "x = " << x << " le test trouve x égal à 0" << endl;
return 0;
}
∗ ∗ ∗ Équation du second degré ∗ ∗ ∗
Calculer les racines réelles de l’équation du second degré à coefficients réels :
ax2 + bx + c = 0
On supposera que a, b, c peuvent avoir des valeurs réelles quelconques, y compris a = b = c = 0, car s’ils sont issus d’un
calcul précédent il faut prévoir tous les cas. Le programme doit afficher la valeur d’une variable entière qui vaut :
8
0 s’il n’y a pas de racine réelle

1 s’il y a une seule racine réelle, simple ou double
2 s’il y a deux racines réelles distinctes
3 dans le cas a = b = c = 0 (équation indéterminée, tout x est solution)
et les valeurs des racines quand il y en a.
∗ ∗ ∗ Températures à Orly en juillet 2015 ∗ ∗ ∗

Copier le fichier temperatures.dat par la commande :
cp /public/mphyo/temperatures.dat .
Ce fichier contient les températures relevées toutes les trois heures à Orly du 1er juillet 2015 à 00h au 31 à 21h. Il se
présente sous la forme :
23.8
19.7
23
...
Les 8 premières lignes sont donc les températures à 0h, 3h, 6h, 9h, 12h, 15h, 18h, 21h le 1er juillet, les 8 suivantes les
températures aux mêmes heures le 2 juillet, etc. Au total le fichier comporte donc 31 × 8 = 248 lignes.
Écrire un programme qui trace sur un même graphe :
- la courbe des minimas journaliers en fonction du jour du mois
- la courbe des maximas journaliers en fonction du jour du mois
Pour le tracé utiliser la notice Graphisme avec Python Matplotlib du site de l’enseignement d’informatique L3-Mag1.
∗ ∗ ∗ Conjecture de Syracuse ∗ ∗ ∗ (facultatif )

Soit u0 un entier positif et la suite :
un+1 = un /2 si un est pair
un+1 = 3un + 1 si un est impair
Calculer les termes de cette suite jusqu’à trouver 1.
(Il existe une conjecture selon laquelle on finit toujours par trouver 1. Pour plus de renseignements chercher « Conjecture
de Syracuse » sur google).
9
4. Fonctions
Remarque
Dans tous les cas on fera afficher les données de départ et les résultats par le programme principal.
∗ ∗ ∗ Force gravitationnelle ∗ ∗ ∗
Le module de la force gravitationnelle exercée par une sphère homogène de rayon R et de masse M sur une masse
ponctuelle m située à une distance d du centre de M , vaut :
Mm M md
G pour d ≥ R et G pour d ≤ R
d2 R3
G étant la constante de la gravitation universelle.
Écrire une fonction d’arguments M , R, m et d qui retourne le module de la force.
∗ ∗ ∗ Combinaisons ∗ ∗ ∗
Transformer en fonction le programme, écrit dans un TD précédent, qui calcule le nombre de combinaisons de k éléments
choisis parmi n, ces deux paramètres étant transmis en argument à la fonction.
∗ ∗ ∗ Fonction créneau ∗ ∗ ∗
Ecrire une fonction dont l’argument est un réel x et dont la valeur est celle de la fonction créneau :
O a x
Figure 1 –
a et b étant deux paramètres. Traiter successivement les deux cas :

a et b sont définis en variables locales dans la fonction
a et b sont définis en variables globales, connus par le programme principal et la fonction.
∗ ∗ ∗ Fonction de répartition ∗ ∗ ∗
1) Écrire une fonction d’arguments x et n, retournant la valeur de :
Z x
1 t2
√ exp(− ) dt
2π 0 2
10
calculée numériquement par la formule des rectangles :

Z b n
X 1
f (x) dx ' h f a + (i − )h
a i=1
2
où n est le nombre de rectangles et h = (b − a)/n le pas.

2) Ajouter une fonction qui effectue un ajustement automatique de la valeur de n de la façon suivante : on part d’une
valeur quelconque de n et on la double tant que la valeur absolue de la variation relative de l’intégrale quand on passe de
n à 2n est supérieure à ε.
3) Sachant que :
Z +∞
1 t2
√ exp(− ) dt = 1
2π −∞ 2
modifier la fonction précédente pour qu’elle retourne la valeur :
Z x
1 t2
√ exp(− ) dt
2π −∞ 2
qui est la fonction de répartition de la loi normale réduite.
∗ ∗ ∗ Lentilles simples ∗ ∗ ∗ (facultatif )

1) Sur un axe Ox on note :
l l’abscisse du centre d’une lentille simple de focale f
x l’abscisse d’un objet
x0 l’abscisse de l’image de l’objet donnée par la lentille
Écrire une fonction qui retourne la valeur x0 et dont les arguments sont l, f et x.
2) On place maintenant sur l’axe Ox :
un objet en O
une lentille de focale f > 0 à l’abscisse d > 0
une lentille de focale F à l’abscisse u avec 0 ≤ u ≤ d
On appelle y l’abscisse de l’image de O donnée par le système des deux lentilles quand on prend le point focal de la
seconde lentille comme origine (il faut donc effectuer une translation de l’origine à la fin du calcul). Écrire les instructions
qui calculent y en fonction de u, F , d et f .
3) On donne les valeurs numériques : f =20mm, d=500mm, F =50mm.
Tracer, avec Python Matplotlib 1 , le graphe de y(u), lorsque u varie entre 0 et d.
4) Au système précédent on ajoute un écran à l’abscisse d + f . Le système lentille de focale f plus écran schématise le
cristallin et la rétine d’un œil. On suppose que l’image paraît nette pour l’œil lorsqu’elle se forme dans une zone comprise
entre la rétine et le plan situé 2mm au-delà de la rétine (donc du côté opposé au cristallin).
Indiquer approximativement pour quelles valeurs de u cela se produit, en portant sur le graphe les deux droites horizon-
tales y = 0 et y=2mm.
5) Ajouter sur le graphe la valeur du grandissement du système en fonction de u.
1. Voir la notice Graphisme avec Python Matplotlib du site de l’enseignement d’informatique L3-Mag1.
11
5. Pointeurs
∗ ∗ ∗ Comparaison variable ordinaire, pointeur ∗ ∗ ∗

Prévoir et expliquer le résultat du programme suivant :
#include<iostream>
int main() {
int i, j;
i = 1; j = i;
cout << "i = " << i << " j = " << j << endl;
i = 2;
cout << "i = " << i << " j = " << j << endl;
int *x, *y;
x = &i;
*x = 1; y = x;
cout << "*x = " << *x << " *y = " << *y << endl;
*x = 2;
cout << "*x = " << *x << " *y = " << *y << endl;
return 0;
}
Prévoir ce qui se passerait si on ne mettait pas l’instruction x = &i;.
∗ ∗ ∗ Passage par adresse d’une variable ordinaire ∗ ∗ ∗

Exemple 1
Soit l’équation en x :
ax + b = 0
Suivant les valeurs de a et b, elle possède zéro, une, ou une infinité de racines. On veut écrire une fonction qui
retourne à la fois un entier et la racine.
Écrire :
a) une fonction de type int :
— valant 1 si l’équation a une seule racine x0 , et 0 sinon
— ayant pour arguments a, b et x0 , où x0 sera donc utilisé pour retourner la valeur de la racine,
s’il y en a.
b) un exemple d’appel de cette fonction.
Exemple 2
Soit la suite :
un+1 = aun + b avec u0 quelconque
Écrire :
a) une fonction de type void ayant trois arguments de type double :
deux en entrée pour transmettre a et b
un qui contient un en entrée et un+1 en sortie
b) un programme principal imprimant les N premiers termes de la suite.
12
∗ ∗ ∗ Dérivée ∗ ∗ ∗
Étant donné une fonction f (x), une valeur approchée de sa dérivée en x est donnée par la formule à cinq points suivante :
1
f 0 (x) = [f (x − 2∆x) − 8f (x − ∆x) + 8f (x + ∆x) − f (x + 2∆x)] (1)
12 ∆x
∆x étant une variation de x choisie par l’utilisateur.
1) Écrire une fonction à laquelle on fournit en argument :
une fonction f (x)
une valeur x0
une valeur de la variation ∆x
et dont la valeur est la valeur approchée de la dérivée en x0 calculée par la formule (1).
(Tester la précision de la formule avec la fonction arc tangente dont la dérivée est :
1
1 + x2
ou tout autre fonction au choix.)
2) Étudier la précision du résultat en fonction de la valeur de ∆x.
13
6. Tableaux dynamiques
∗ ∗ ∗ Matrice ∗ ∗ ∗
Soit la matrice suivante :
0.1 10.31 8.7 2
2841.573 54 6.57 185.34
5.03 0.73 1065 0.07
L’utiliser pour initialiser un tableau et la faire afficher à l’écran.
∗ ∗ ∗ Remplirtab ∗ ∗ ∗
Écrire un programme qui remplit un tableau d’entiers n(10,9) de la façon suivante :
000000000 ______ second indice
100000001
200090002
300090003
400090004
500090005
600090006
700090007
800090008
999999999
|
|
|
premier indice
Pour vérifier le résultat, faire afficher ce tableau à l’écran.
∗ ∗ ∗ Norme d’un vecteur ∗ ∗ ∗

Les composantes d’un vecteur ~v dans un repère orthonormé à n dimensions sont définies dans un programme principal
par un tableau v.
Ecrire une fonction de v dont la valeur est la norme de ~v .
Dans le programme principal faire deux appels successifs à la fonction avec deux valeurs de n différentes.
∗ ∗ ∗ Vecteur unitaire ∗ ∗ ∗
Les composantes d’un vecteur ~v dans un repère orthonormé à n dimensions sont définies dans un programme principal
par un tableau v à n éléments.
1) Écrire une fonction de type int qui :
si ~v n’est pas nul, vaut 1 et retourne en argument les composantes du vecteur unitaire ~u parallèle à ~v et de
même sens que lui dans un tableau u
si ~v est nul vaut 0 et retourne le vecteur nul au lieu d’un vecteur unitaire dans le tableau u.
Comment faire pour interdire à priori toute modification des valeurs de v dans la fonction ?
2) Dans le programme principal, faire deux appels successifs de la fonction avec deux vecteurs de dimensions différentes.
∗ ∗ ∗ Trace d’une matrice ∗ ∗ ∗

Écrire une fonction dont la valeur est la somme des éléments diagonaux d’une matrice carrée de dimension quelconque
14
fournie en argument.
∗ ∗ ∗ Produit de matrices ∗ ∗ ∗
Écrire une fonction à laquelle on passe en argument une matrice a à m lignes et n colonnes et une matrice b à n lignes et
p colonnes et qui retourne en argument la matrice produit.
∗ ∗ ∗ Résolution de l’équation de Laplace : lentille électrostatique à fente ∗ ∗ ∗

Le problème étudié est à deux dimensions. Dans une enceinte conductrice rectangulaire au potentiel 0 sont placés des
conducteurs de potentiels imposés.
0 30 50 55 75 105
y (indice j)
20
0 volts 1000 volts
40
0 volts 1000 volts Toute l’enceinte
est a 0 volts
60
x
Les cotes sont exprimees en prenant
(indice i) comme unite le pas du reseau
Figure 2 – Agrandir cette figure à l’écran pour bien distinguer les traits pleins des traits pointillés
Le but de l’exercice est de calculer le potentiel électrostatique V (x, y) à l’intérieur de l’enceinte, hormis sur les conducteurs.
Il s’agit de déterminer numériquement la solution unique de l’équation :

∆V = 0 entre les conducteurs
V a les valeurs imposées sur les conducteurs
On effectue le calcul de V uniquement aux points d’un réseau régulier 2 appliqué sur l’enceinte, c’est à dire en des points
tels que x = i et y = j, i et j étant des entiers nuls ou positifs. La quantité V (i, j) est notée Vi,j . Dans ces conditions la
solution est obtenue de la façon suivante :
1. pour tous les points situés sur un conducteur, y compris l’enceinte, on attribue à V la valeur du potentiel de ce
conducteur
2. pour tous les autres points de l’enceinte on attribue à V une valeur initiale nulle puis on calcule tous les Vi,j par
la formule :
Vi−1,j + Vi,j−1 + Vi+1,j + Vi,j+1
Vi,j = (moyenne arithmétique des quatre points adjacents)
4
3. on réitére l’étape 2. jusqu’à obtenir une convergence satisfaisante.
1) Écrire un programme qui calcule le potentiel pour un nombre d’itérations fixé par l’utilisateur et imprime les valeurs
finales dans un fichier sous la forme :
V0,0 V0,1 ... V0,j ... V0,105
V1,0 V1,1 ... V1,j ... V1,105
... ... ... ... ... ...
Vi,0 Vi,1 ... Vi,j ... Vi,105
... ... ... ... ... ...
V60,0 V60,1 ... V60,j ... V60,105
2. dont on choisit le pas comme unité de longueur
15
2) Compléter le programme précédent pour que le calcul s’arrête dès que l’une des deux conditions suivantes est réalisée :
le nombre d’itérations fixé par l’utilisateur est dépassé
|Vi,j (itération n+1) − Vi,j (itération n)| ≤ ε pour chacune des valeurs de i et j hors conducteurs, ε étant choisi par
l’utilisateur.
3) Tracer la surface V (x, y) avec Python Matplotlib 3 sous forme de nappe puis sous forme de courbes de niveau.
∗ ∗ ∗ Échantillonnage d’une fonction ∗ ∗ ∗

On appelle échantillonnage de la fonction f (x) n valeurs de cette fonction correspondant à n valeurs de x régulièrement
espacées entre les bornes xmin et xmax .
1) Écrire une fonction ech(f,xmin,xmax,n,t) qui calcule le tableau d’échantillonnage t, f étant la fonction à échan-
tillonner.
2) Donner ensuite un exemple d’utilisation de cette fonction en écrivant :
un programme principal complet qui l’appelle en lui fournissant les arguments nécessaires
une fonction f test
3) Ajouter dans le programme principal l’écriture du tableau d’échantillonnage dans un fichier.
∗ ∗ ∗ Interpolation linéaire ∗ ∗ ∗
On suppose qu’on ne connaît une fonction mathématique f (x) qu’en n points d’abscisses xi . Les valeurs des xi et les
valeurs correspondantes yi = f (xi ) de la fonction sont données dans un fichier organisé de la façon suivante :
sur la première ligne n
sur les n suivantes les couples xi yi
Les xi ne sont pas nécessairement régulièrement espacées mais sont rangées par ordre croissant. L’abscisse minimum est
donc x0 et l’abscisse maximum xn−1 .
Exemple d’un tel fichier :
8
0.2 2.5
0.35 0.198669
...
0.67 -8.472
1) Écrire un programme qui stocke les valeurs des xi dans un tableau à un indice et les valeurs des yi dans un autre
tableau à un indice, puis qui calcule la valeur interpolée linéairement de f en un point d’abscisse x quelconque comprise
entre x0 et xn−1 .
2) Transformer la partie du programme qui fait le calcul d’interpolation en une fonction dont les arguments sont n, les
xi , les yi , x et qui retourne la valeur interpolée de f (x).
∗ ∗ ∗ Trinum ∗ ∗ ∗
On donne une liste de n réels dans un tableau. Réécrire cette liste dans le même tableau par ordre croissant. On pensera
à ce que l’on fait quand on remet un jeu de cartes dans l’ordre.
La liste des réels sera fournie dans le programme par initialisation du tableau.
3. Voir la notice Graphisme avec Python Matplotlib du site de l’enseignement d’informatique L3-Mag1.
16
∗ ∗ ∗ Statistique des tailles d’une population d’individus ∗ ∗ ∗

On a mesuré la taille de chacun des individus d’une population donnée et on a écrit ces tailles, exprimées en mètres, dans
un fichier, à raison d’une taille par ligne. On dispose donc d’un fichier du type :
1.625
1.576
1.733
...
On ne sait pas à priori combien de tailles, donc de lignes, contient le fichier.

On demande d’écrire un programme qui :
1) donne le nombre
√ de tailles mesurées, la valeur minimale et la valeur maximale, la valeur moyenne < t >,
et l’écart type < t2 > − < t >2
2) constitue dans un tableau la distribution statistique des tailles par tranches de un centimètre.
∗ ∗ ∗ Divergence et rotationnel d’un champ de vecteurs ∗ ∗ ∗ (facultatif )

On considère un champ de vecteurs :
→
− →
V (−
r)
dans l’espace à trois dimensions et on veut en calculer numériquement la divergence et le rotationnel en un point quelconque
(→
−
r est le vecteur position).
Le calcul numérique des dérivées sera fait par la formule approchée :
f (x + h) − f (x − h)
2h
h étant le pas de dérivation numérique à choisir par l’utilisateur.
Ce champ de vecteurs est défini par une fonction du type :
void champ(double *x,double *v) {

v[0] = une fonction quelconque de x[0],x[1],x[2]
v[1] = une seconde fonction quelconque de x[0],x[1],x[2]
v[2] = une troisième fonction quelconque de x[0],x[1],x[2]
}
→
−
où v[0],v[1],v[2] désignent les trois composantes du vecteur V et x[0],x[1],x[2] celles du vecteur →
−
r.
Écrire une fonction qui a pour arguments d’entrée :
un champ de vecteurs, c’est à dire une fonction de même prototype que champ
la position à laquelle on veut calculer la divergence et le rotationnel, c’est à dire le tableau x
le pas de dérivation numérique h.
et pour arguments de sortie :
la divergence
le rotationnel.
Écrire également un exemple de programme principal appelant cette fonction.
17
7. Générateur aléatoire
∗ ∗ ∗ Lancer de pièces ∗ ∗ ∗
On lance n pièces identiques de rayon r dans le fond d’une boîte carrée de côté 1. Les pièces tombent à plat et peuvent
se chevaucher. Si une des pièces heurte le bord de la boîte on ne compte pas le lancer.
Figure 3 –
On veut calculer la proportion de la surface du fond de la boîte qui est recouverte en moyenne par les pièces. Pour cela :
- simuler un lancer de pièces en tirant au hasard les coordonnées xi et yi des centres des n pièces.
- calculer la proportion de la surface recouverte par un lancer donné en tirant au hasard un très grand nombre de
points dans la boîte et en comptant la proportion de ceux qui se trouvent sur une pièce.
- répéter l’expérience un grand nombre de fois pour trouver la valeur moyenne de cette proportion.
Indiquer quelques conditions dans lesquelles on peut tester la justesse du programme.
Utiliser la fonction drand48() (nécessite la directive #include<stdlib.h>) dont la valeur varie aléatoirement, à chaque
appel, dans l’intervalle [0, 1] avec une distribution uniforme. Pour initialiser la suite des valeurs aléatoires à une valeur
qui sera différente à chaque exécution du programme (à savoir le nombre de secondes depuis le 1er janvier 1970, c’est à
dire l’heure UNIX), mettre la commande srand48(time(NULL)); au début du programme (ce qui nécessite en plus la
directive #include<time.h>).
18
8. Équations différentielles
∗ ∗ ∗ Étude du pendule simple sans faire l’approximation des petits angles ∗ ∗ ∗
z
S
θ m
O y
x
Figure 4 –
Un pendule simple est constitué d’une tige rigide sans masse, libre de tourner autour du point fixe S de l’axe Oz et d’une
masse m située à l’autre extrémité de la tige.
−→ −→
La position et la vitesse initiale sont contenues dans le plan yOz. On pose θ = (SO, Sm). Tracer la courbe de θ en fonction
du temps, sans faire l’approximation des petits angles :
a) par la méthode d’Euler en programmant soi-même la formule d’Euler
b) par la méthode de Runge-Kutta d’ordre 4, en utilisant la fonction rk4 4 décrite dans le cours.
Tracer sur un même graphe les deux courbes obtenues aux a) et b) et celle résultant de l’approximation des petits angles.
4. Pour pouvoir utiliser cette fonction, il suffit d’ajouter la directive #include<bibli_fonctions.h> et de compiler et exécuter avec ccc.
19
L3-Mag1 Phys. fond. 16-17 TD Info : emacs 4 octobre 2016
L’éditeur emacs
Introduction
Un éditeur de texte permet d’écrire des caractères dans un fichier à partir du clavier et de la souris.
Toutes les commandes emacs peuvent être effectuées à partir du clavier, mais seules certaines d’entre elles peuvent l’être à l’aide de
la souris, en utilisant les menus de la barre supérieure. Les commandes au clavier commencent par une des deux touches particulières
Ctrl ou Alt, suivie d’un caractère. Par exemple :
Ctrl-d signifie : maintenir la touche Ctrl enfoncée et appuyer sur la touche d (ce qui a pour effet d’effacer le caractère
situé à l’emplacement du curseur).
Alt-< signifie : maintenir la touche Alt enfoncée et appuyer sur la touche < (ce qui a pour effet de placer le curseur au
début du fichier).
Il est à noter que la combinaison spéciale Ctrl-x est toujours suivie par un deuxième caractère (avec ou sans Ctrl), tandis que la
combinaison spéciale Alt-x permet de taper le nom explicite d’une commande dans la barre de dialogue.
Dans la suite on notera C pour la touche Entrée.
Ouvrir et fermer une session avec emacs

Ouvrir
Pour créer un nouveau fichier ou ouvrir un fichier déjà existant, de nom fich.cpp par exemple, écrire sur la ligne de commande
Linux, après l’invite symbolisée par $ :
$ emacs fich.cpp & C
qu’on peut, dans notre environnement, écrire en abrégé :
$ em fich.cpp C
Ceci a pour effet d’ouvrir une nouvelle fenêtre dans laquelle on va pouvoir écrire. Elle comporte :
dans sa partie supérieure une barre de menus
dans sa partie inférieure :
une barre affichant le nom du fichier, l’heure et le numéro de ligne de l’emplacement du curseur
une barre affichant le dialogue emacs -utilisateur.
Lorsque emacs pose une question dans la barre de dialogue il faut y répondre en terminant par Entrée.
Fermer
Sauver et fermer avec le même nom de fichier Ctrl-x Ctrl-c
Sauver sans fermer avec le même nom de fichier Ctrl-x Ctrl-s
Sauver sans fermer dans le fichier de nom lili (il y a demande du nom du fichier) Ctrl-x Ctrl-w lili C
Ce qu’on écrit à l’écran est stocké dans une zone de la mémoire appelée buffer et n’est pas enregistré dans un fichier tant que l’on
n’a pas sauvé d’une des trois façons précédentes. En cas d’arrêt anormal, une partie ou la totalité du contenu du buffer peut donc
être perdue. A chaque fois qu’on sauve, la version précédente est conservée dans un fichier portant le même nom suivi d’un tilde :
ici, par exemple, fich.cpp∼. Si on ferme sans sauver (dans ce cas une double confirmation est demandée) le contenu du buffer non
sauvé est conservé dans un fichier de même nom encadré de dièses : ici, par exemple, #fich.cpp# 5
Se déplacer dans le buffer

Le numéro de la ligne sur laquelle est situé le curseur est indiqué sur la barre inférieure de l’écran (L15 pour ligne 15).
Déplacer le curseur vers le haut, à gauche, vers le bas, à droite ↑←↓→

Avancer (reculer) d’un mot Alt-f (b)
Aller en début de ligne Ctrl-a
Aller en fin de ligne Ctrl-e
Avancer d’un écran Ctrl-v
Reculer d’un écran Alt-v
Aller au début (à la fin) du buffer Alt-< (>)
Placer le curseur à la 15ème ligne Alt-x goto-line C 15 C
Placer le curseur au milieu de l’écran Ctrl-l
5. Pour avoir accès à ce fichier dans une commande Linux il faut mettre des antislash avant les dièses : \#fich.cpp\#
20
Supprimer des caractères

Supprimer le caractère à l’emplacement du curseur Ctrl-d
Supprimer le caractère précédent Backspace
Supprimer depuis le curseur inclus jusqu’à la fin de la ligne Ctrl-k
Supprimer le mot précédent Alt-Backspace
Supprimer le mot suivant Alt-d
Définir une région, la supprimer ou la déplacer

Définir une région 1) placer le curseur au début de la région
et Ctrl-barre d’espace
2) placer le curseur immédiatement après
la fin de la région
Mettre la région en mémoire et la supprimer (couper) Ctrl-w
Mettre la région en mémoire sans la supprimer (copier) Alt-w
Recopier la région mise en mémoire à l’emplacement du curseur (coller) Ctrl-y
Chercher
Chercher la chaîne zzz vers l’avant (l’arrière) à partir du curseur Ctrl-s (r) zzz
(puis Ctrl-s (r) pour chercher l’endroit suivant)
Chercher et remplacer avec confirmation la chaîne zzz par la chaîne yyy Alt-% zzz C yyy C
et répondre y pour remplacer, n sinon
et q pour quitter
Chercher et remplacer sans confirmation la chaîne zzz par la chaîne yyy Alt-x replace-string C zzz C yyy C
Echanges avec d’autres fichiers

Insérer le contenu du fichier lili à l’endroit du curseur Ctrl-x i lili C
Copier une région dans le fichier lili Alt-x write-region C lili C
Détruire le fichier lili Alt-x delete-file C lili C
Renommer le fichier lili en lulu Alt-x rename file C lili C lulu C
Définir un rectangle, le supprimer ou le déplacer

Définir un rectangle comme pour une région
Mettre le rectangle en mémoire et le supprimer Ctrl-x r k
Mettre le rectangle en mémoire sans le supprimer Ctrl-x r Alt-w
Recopier le rectangle mis en mémoire à l’emplacement du curseur Ctrl-x r y
Supprimer le rectangle sans le mettre en mémoire Ctrl-x r d
Indenter
Indenter une région d’un fichier C Alt-x indent-region C
Idem Alt-Ctrl-\
Ajouter une chaîne en début ou en fin de ligne

Ajouter la chaîne zzz en début (fin) de ligne Alt-x query-replace-regexp C ˆ ($) C zzz C
Recherche de commande par complétion

Cette méthode est très utile si on ne connaît pas le nom exact de la commande ou s’il est long à écrire.
Avoir la liste des commandes dont le nom commence par la chaîne zzz (longueur qqc) Alt-x zzz Tab (touche Tabulation)
Avoir la liste de toutes les commandes Alt-x Tab
Dans la liste affichée il suffit ensuite de cliquer sur la commande désirée pour la faire s’exécuter.
Par exemple si on veut ajouter la chaîne zzz en fin de ligne, on tape :
21
Alt-x qu Tab
emacs affiche alors la liste de toutes les commandes commençant par qu et, dans cette liste, se trouve la commande query-replace-
regexp sur laquelle il suffit de cliquer pour que la commande s’exécute.
En cas d’erreur
Interrompre une commande en cours Ctrl-g
Annuler les derniers changements un par un Ctrl-x u ou Ctrl-_
Annuler tous les changements depuis la dernière sauvegarde Alt-x revert-buffer C
Revenir à la version initiale, avant ouverture Ctrl-x Ctrl-f fich.cpp∼ C
Ctrl-x Ctrl-w fich.cpp C
Accéder à l’aide en ligne Ctrl-h
Débloquer l’écran Ctrl-q
Bibliographie : Emacs Debra Cameron O’REILLY collection précis & concis
22
L3-Mag1 Phys. fond. 16-17 Commandes Linux 4 octobre 2016
Commandes Linux
Dans ces pages on donne une liste de commandes courantes avec le minimum d’explications dans le but d’une utilisation
simple et immédiate.
Linux distingue entre les minuscules et les majuscules.
Les fichiers sont rangés dans des répertoires répartis en arborescence.
Les commandes doivent être écrites dans la fenêtre Terminal, après l’invite, et validées par un Entrée.
Les commandes cp, rm, mv des tableaux suivants sont accompagnées d’une demande de confirmation si on leur ajoute
l’option -i. Exemple :
rm fich destruction sans confirmation
rm -i fich destruction avec confirmation.
Les commandes suivies du signe (+) ont été ajoutées, elles sont propres à l’environnement du L3/Magistère.
Ctrl-b signifie : maintenir la touche Ctrl enfoncée et appuyer sur la touche b (ne signifie pas qu’il faut taper de - entre
Ctrl et b).
23
Commandes les plus courantes

Faire la liste des fichiers et des répertoires du répertoire courant,
sauf ceux qui commencent par un point ls
Faire la liste de tous les fichiers et répertoires du répertoire courant ls -a
Idem mais les noms de répertoires sont suivis d’un /
et ceux des fichiers exécutables d’une ∗ ls -aF
Idem avec plus d’informations et par ordre chronologique inverse ls -altr
Passer du répertoire courant au répertoire fils rep cd rep
Ouvrir le site du cours dans un navigateur site (+)
Editer le fichier fich emacs fich & ou em fich (+)
(voir la notice de l’éditeur emacs )
Faire défiler le contenu du fichier fich d’un seul coup cat fich
Faire défiler le contenu du fichier fich page par page less fich
Passer à la page suivante Barre d’espace
Passer à la ligne suivante Entrée
Passer à la page précédente b
Quitter q
Copier le contenu de fich1 dans fich2 cp fich1 fich2
(fich2 est écrasé s’il existait déjà, fich1 est conservé)
Changer le nom du fichier fich1 en fich2 mv fich1 fich2
(fich2 est écrasé s’il existait déjà)
Déplacer le fichier fich dans le répertoire fils rep mv fich rep
( fich est écrasé s’il existait déjà dans rep)
Détruire le fichier fich rm fich ou del fich (+)
Créer le répertoire rep mkdir rep
Détruire le répertoire rep rmdir rep
(Il faut d’abord l’avoir vidé de tout fichier)
Compiler prog.cpp g++ -lm prog.cpp
Exécuter a.out
Compiler et exécuter prog.cpp ccc prog.cpp (+)
Interrompre et tuer un processus en cours Ctrl-c
Passer un processus en cours en arrière-plan Ctrl-z bg
Complétion et rappel de commandes

On peut ne taper que les premières lettres d’une commande ou d’un nom de fichier du répertoire courant, suivies
de Tab. S’il n’y a pas d’ambiguïté le système complète tout seul, sinon il propose une liste des possibilités.
On peut également faire un copier-coller sur les noms de fichiers et de répertoires apparaissant dans la fenêtre
du terminal dans lequel on écrit les commandes Linux. Il faut double-cliquer avec le bouton de gauche de la
souris sur la chaîne de caractères que l’on veut saisir, puis cliquer sur le bouton du milieu.
On peut rappeler les commandes écrites antérieurement à l’aide des flèches verticales du clavier.
Il faut utiliser ces trois possibilités systématiquement, ce qui fait gagner beaucoup de temps.
24
Commandes un peu moins courantes

Ajouter le contenu du fichier fich2 au fichier fich1 cat fich2 >> fich1
Rediriger dans le fichier fich le résultat de la commande cmd cmd > fich
Caractère générique ∗
Il remplace n’importe quelle chaîne de caractères dans un nom de fichier ou de répertoire.
Désigner tous les fichiers dont le nom commence par si et se termine par do si∗do
Ceux qui commencent par fa fa∗
Ceux qui se terminent par sol ∗sol
etc... ∗re∗
mi∗fa∗sol
Tous les fichiers ∗
Très commode mais à utiliser avec précaution pour des opérations irréversibles :
tester d’abord le résultat obtenu avec la commande echo (sans risque).
Très utile aussi pour désigner un fichier unique en abrégé.
Faire le ménage
c’est à dire détruire (avec confirmation) tous les fichiers commençant par #,
finissant par ∼, les fichiers core et a.out pu (+)
Créer le fichier de commandes fichcom emacs fichcom &
Créer le fichier de commandes fichcom dans le répertoire bin à partir
du répertoire d’accueil et y placer la liste des commandes que l’on veut faire exécuter emacs bin/fichcom &
Rendre fichcom exécutable chmod 755 fichcom
Faire exécuter la liste de commandes fichcom
Créer l’abréviation com pour la commande commande alias com=’commande’
Pour le rendre permanent il faut placer cet alias dans le fichier .bashrc
Faire apparaître les différences entre les fichiers fich1 et fich2 diff fich1 fich2
ou kompare fich1 fich2 &
Demander une description détaillée de la commande com man com
Faire la liste de toute l’arborescence à partir du répertoire initial du -a
Faire afficher toutes les lignes du fichier fich qui contiennent
la chaîne chaîne grep chaîne fich
Faire la liste de tous les utilisateurs connectés who -u
Rechercher si l’utilisateur util est connecté who -u | grep util
Rechercher où est placée la commande com which com
Rechercher tous les noms de fichiers dont le nom contient la chaîne toto locate toto
25
Comment corriger une ligne de commandes ?

Les commandes sont les mêmes que celles de emacs . Ce qui est utile ici :
Aller au début de la ligne de commande Ctrl-a
Aller à la fin de la ligne de commande Ctrl-e
Reculer d’un caractère Ctrl-b ou ←
Avancer d’un caractère Ctrl-f ou →
Détruire le caractère sous le curseur Ctrl-d ou Del/Suppr
Détruire le caractère précédent Ctrl-h ou Backspace
Effacer la fin de la ligne Ctrl-k
Effacer toute la ligne Ctrl-u
Echanger le caractère sous le curseur et celui qui le précède Ctrl-t
Comment rappeler une ligne de commandes frappée précédemment ?

Faire défiler les lignes précédentes en remontant le temps ↑
Faire défiler les lignes précédentes en suivant le temps ↓
Conventions pour les noms de fichier

Il est conseillé d’adopter les conventions suivantes pour les noms de fichier :
C ou C++ Résultats Données
fich.cpp fich.res fich.dat
26

TD Enonces

Transféré par

Droits d'auteur :

Formats disponibles

TD Enonces

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

TD Enonces

Transféré par

Droits d'auteur :

Formats disponibles

L3-Mag1 Phys. fond.

16-17 Informatique 4 octobre 2016

L3 et Magistère 1ère année, Physique fondamentale, année 2016-2017

Travaux dirigés de langage C appliqué à la

— Textes d’exercices de TD (n◦ 1 à 8)

1. Types entier et réel, priorité des opérations

— Commandes Linux de base

1. Types entier et réel, priorité des opérations

Types entier et réel

Règles de priorité des opérations

Usage des parenthèses et des variables intermédiaires

2. Limitations des valeurs entières et réelles

Diagnostics en cas d’opérations illégales

Limitations en grandeur et précision

Pour les réels

ici que le type double.

Type Nb de Nb de Nb de bits Plus petite Plus grande Nb minimum de

double 1 11 53 ' 2.225e-308 ' 1.798e308 15

∗ ∗ ∗ Étalement d’un paquet d’onde ∗ ∗ ∗

b) Même question pour un proton (chercher la masse du proton sur internet).

Diagnostics en cas d’opérations illégales

Cas Noter ici le résultat obtenu

Entier trop grand ou trop petit

Réel trop petit en valeur absolue

Réel trop grand en valeur absolue

Prendre la racine carrée d’un réel négatif

b(k) = Cnk pk (1 − p)n−k

0 s’il n’y a pas de racine réelle

∗ ∗ ∗ Températures à Orly en juillet 2015 ∗ ∗ ∗

∗ ∗ ∗ Conjecture de Syracuse ∗ ∗ ∗ (facultatif )

a et b étant deux paramètres. Traiter successivement les deux cas :

calculée numériquement par la formule des rectangles :

où n est le nombre de rectangles et h = (b − a)/n le pas.

qui est la fonction de répartition de la loi normale réduite.

∗ ∗ ∗ Lentilles simples ∗ ∗ ∗ (facultatif )

∗ ∗ ∗ Comparaison variable ordinaire, pointeur ∗ ∗ ∗

∗ ∗ ∗ Passage par adresse d’une variable ordinaire ∗ ∗ ∗

Pour vérifier le résultat, faire afficher ce tableau à l’écran.

∗ ∗ ∗ Norme d’un vecteur ∗ ∗ ∗

∗ ∗ ∗ Trace d’une matrice ∗ ∗ ∗

∗ ∗ ∗ Résolution de l’équation de Laplace : lentille électrostatique à fente ∗ ∗ ∗

∗ ∗ ∗ Échantillonnage d’une fonction ∗ ∗ ∗

∗ ∗ ∗ Statistique des tailles d’une population d’individus ∗ ∗ ∗

On ne sait pas à priori combien de tailles, donc de lignes, contient le fichier.

∗ ∗ ∗ Divergence et rotationnel d’un champ de vecteurs ∗ ∗ ∗ (facultatif )

void champ(double *x,double *v) {

∗ ∗ ∗ Étude du pendule simple sans faire l’approximation des petits angles ∗ ∗ ∗

Ouvrir et fermer une session avec emacs

Se déplacer dans le buffer

Déplacer le curseur vers le haut, à gauche, vers le bas, à droite ↑←↓→

Supprimer des caractères

Définir une région, la supprimer ou la déplacer

Echanges avec d’autres fichiers

Définir un rectangle, le supprimer ou le déplacer

Ajouter une chaîne en début ou en fin de ligne

Recherche de commande par complétion

Bibliographie : Emacs Debra Cameron O’REILLY collection précis & concis

Commandes les plus courantes

void champ(double x,double v) {