Discussion :

[najoua01]

bonjour,
je suis débutante j'ai une petite projet d'equation de la chaleur en C bon voici mon code

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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
#include<conio.h>
#include<stdio.h>
#include<math.h>
int main()
{
int i,j,n,m;
float dx,dt,ni=1,u[100][100],tmax,f[100][100],landa,;
float pi=3.141592;
FILE *g;
printf("donner le pas d'espace:\n");
scanf("%f",&dx);
printf("donner la largeur de l'intervalle du temps:\n");
scanf("%f",&tmax);
printf("donner le pas du temps:\n");
scanf("%f",&dt);
n=(int)(1/dx)+1;
m=(int)(tmax/dt);
//d‚finition de la fonction f
for(j=0;j<=m;j++)
for(i=0;i<=n;i++)
f[i][j]=2;
//conditions aux limites
for(j=0;j<=m;j++)
u[0][j]=u[n][j]=0;
//conditions initiales
for(i=1;i<n;i++)
u[i][0]=sin(pi*i*dx)+i*dx*(1-i*dx);
//sch‚ma d'Euler progressif
landa=dt*ni/(dx*dx);
for(j=0;j<m;j++)
for(i=1;i<n;i++)
u[i][j+1]=(1-2*landa)*u[i][j]+landa*(u[i-1][j]+u[i+1][j])+dt*f[i][j];
//affichage des r‚sultats
g=fopen("c:\\EUexplicite.txt","w");
for(j=0;j<=m;j++){
for(i=0;i<=n;i++)
fprintf(g,"%2.3f\t",u[i][j]);
fprintf(g,"\n\n");
}
fclose(g);
getch();
return 0;
}

mais le problème que après l'exécution il me donne pas des vrais valeur il m'affiche d'autre valeur
voici ce que je doit afficher

[ppphil]

Si c'est du C, t'auras plus de chances dans le forum C.... Non ?

[najoua01]

voici la photo
oui j suis sur que c'est C

[jbrasselet]

Tu devrais plutôt poster là alors :
http://www.developpez.net/forums/sea...archid=6703594

[najoua01]

quelqu'un peu m'aider pour términer ce mini projet
merci d'avance

[souviron34]

2 remarques vite faites :

• tu prends les boucles en sens inverse...
  
  tab[][j] s'explore avec
  
  

  Code :

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

  1
  2
  3
  4
  5
  6
  7
  for ( i = 0 ; i <  nli ; i ++ )
  {
     for ( j = 0 ; j < nco ; j++ )
    {
         tab[i][j] = ... 
     }
  }

• D'autre part tu n'as pas mis les parenthèses permttant de définir les blocs
  
  

  Code :

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

  1
  2
  3
  for(j=0;j<=m;j++)
  for(i=0;i<=n;i++)
  f[i][j]=2;

  revient à faire :
  
  

  Code :

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

  1
  2
  for ( i = 0 : i <= n ; i++ )
    f[i][n+1]=2;

  car la première boucle est effectuée sans rien, et e termine avec une valeur de j égale à n+1
  
  
  
• Dans toute initialisation d'un float ou double, il est préférable de le faire suivre d'un ".0" si c'est une valeur entière
  
  

  Code :

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

  float dx = 2 ;

  ne donnera pas forcément le bon résultat, car la seconde partie du nombre n'est pas définie..
  
  

  Code :

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

  float dx=2.0 ;

  donnera ce qu'on souhaite réellement.

[diogene]

D'autre part tu n'as pas mis les parenthèses permttant de définir les blocs

Code :

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

1
2
3
for(j=0;j<=m;j++)
for(i=0;i<=n;i++)
f[i][j]=2;

revient à faire :

Code :

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

1
2
for ( i = 0 : i <= n ; i++ )
  f[i][n+1]=2;

car la première boucle est effectuée sans rien, et e termine avec une valeur de j égale à n+1

Non, le code est équivalent à

Code :

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

1
2
3
4
5
6
7
for(j=0;j<=m;j++)
{
   for(i=0;i<=n;i++)
   {
      f[i][j]=2;
   }
}

Dans toute initialisation d'un float ou double, il est préférable de le faire suivre d'un ".0" si c'est une valeur entière

Code :

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

float dx = 2 ;

ne donnera pas forcément le bon résultat, car la seconde partie du nombre n'est pas définie..

Code :

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

float dx=2.0 ;

donnera ce qu'on souhaite réellement.

Non, le résultat est correct : le compilateur sait faire un conversion entier-> flottant exactement (si l'entier n'est pas très grand, mais avec 2 on est plus que tranquille).

[souviron34]

pour la première remarque, exact, j'avas lu trop vite J'avais cru voir un ";"...

Vraiment désolé..


Quant à la deuxième, je ne suis pas d'accord, car même si c'est vrai pour 2, il est de toutes façons plus prudent d'avoir une attitude consistante quelle que soit la valeur d'initialisation.... Et cela a l'avantage également de confirmer d'un coup d'oeil à la lecure qu'on a bien affaire à un réel...

En tous cas c'est ce que je prône comme "bonne pratique"...

[diogene]

Quant à la deuxième, je ne suis pas d'accord, car même si c'est vrai pour 2, il est de toutes façons plus prudent d'avoir une attitude consistante quelle que soit la valeur d'initialisation.... Et cela a l'avantage également de confirmer d'un coup d'oeil à la lecure qu'on a bien affaire à un réel...

En tous cas c'est ce que je prône comme "bonne pratique"...

je t'accorde que c'est une bonne pratique et dans ce cas il vaut mieux spécifier une constante en accord avec le type :

Code :

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

float dx = 2.0f;

Mais on ne pouvait dire

....ne donnera pas forcément le bon résultat, car la seconde partie du nombre n'est pas définie..

[souviron34]

c'est vrai..

Mais au vu du nombre de posts où la division entière est en jeu (à l'étonnement des posteurs) à cause de ce genre de choses, il serait à mon avis bon, surtout pour un débutant, d'insister sur ce point...

Mais par rapport au pbe du PO, je suis d'accord que ça n'est pas le cas..