Discussion :

[Jerome Briot]

Il suffit de mettre bout à bout tous les codes.

Si tu exécutes le code ci-dessous, tu verras qu'il prend énormément de temps pour afficher toutes les lignes dans la fenêtre de commande d'Octave.
Il faudra appuyer sur Ctrl+C pour interrompre le code.

Je ne vois pas bien l'intérêt d'afficher autant de résultat avec fprintf…

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
clc
clear all;
%%Définition des constantes :
%Une année en secondes
An = 365.25*24*3600;
 
%Paramètres du maillage
Nz = 81; %de 0 à 80km PROFONDEUR
Nt = 3001; %de 0 à 45Ma TEMPS
 
%Cette boucle permets de définir une température de 825°C entre 0 et 11km à t=0
for i=1:12;
  T(i,1)=825;
end
%Cette boucle permets de définir une température de 100°C entre 13 et 80km à t=0
for i=13:Nz-1;
  T(i,1)=100;
end
%Cette boucle permets de définir une température de 100°C à 80km en tout temps
for j=1:Nt;
  T(Nz,j)=100;
end
 
for i = 1:Nz
  for j = 1:Nt
    fprintf("Temps = %d ; profondeur = %d ; température = %.3f\n", i, j, T(i,j)); 
  endfor  
endfor

[Thibautswiitter]

Super ! C'est exactement ce que je voulais ! Merci !
Du coup, j'ai ajouté mes conditions et j'aimerais savoir si c'est possible maintenant de convertir les valeurs dans la boucle en valeur réelle ?
Je m'explique : dans la boucle, j'ai utilisé les valeurs du pas en temps et du pas en profondeur

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
%Pas en profondeur 
Nz = 81 %nombres de points pour aller de 0 à 80 km
%Pas en temps
Nt = 3001 %nombres de points pour aller de 0 à 45 Ma

Et j'ai utilisé le code suivante :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
clc
 
for i = 1:Nz
 if ((i==7) || (i==9) || (i==11) || (i==14) || (i==19) || (i==45));
  for j = 1:Nt
   if ((j==1) || (j==15) || (j==30) || (j==75) || (j==150) || (j==1500) || (j==3001));
    fprintf("profondeur = %d; temps = %d; température = %.3f°C\n", i, j, T(i,j)); 
   endif
  endfor 
 endif 
endfor

Est-ce que c'est possible de convertir les valeurs dans la boucle ? Ou plutôt à la sortie pour passer de j==3001 à j==45 Millions d'années
Encore merci pour vos réponses ^^

[Jerome Briot]

Ou plutôt à la sortie pour passer de j==3001 à j==45 Millions d'années

Tu dois trouver comment passer d'un intervalle [1 3001] à un intervalle [0 45] et appliquer la formule à chaque valeur de j affichée.

Tu trouveras un début de réponse ici : Comment normaliser les valeurs d'une matrice dans l'intervalle [0 1] ?