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.