programme d'equation de la chaleur

Version imprimable

1 pièce(s) jointe(s)

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

Code:

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

28/01/2011, 12h54
ppphil

Si c'est du C, t'auras plus de chances dans le forum C.... Non ?
28/01/2011, 13h07
najoua01

1 pièce(s) jointe(s)

voici la photo
oui j suis sur que c'est C
28/01/2011, 13h08
jbrasselet

Tu devrais plutôt poster là alors :
http://www.developpez.net/forums/sea...archid=6703594
28/01/2011, 15h23
najoua01

quelqu'un peu m'aider pour términer ce mini projet
merci d'avance
29/01/2011, 14h05
souviron34
2 remarques vite faites :
- tu prends les boucles en sens inverse...
  
  tab[][j] s'explore avec
  Code:
  
  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:
  
  1 2 3 for(j=0;j<=m;j++) for(i=0;i<=n;i++) f[i][j]=2;
  revient à faire :
  Code:
  
  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:
  
  float dx = 2 ;
  
  ne donnera pas forcément le bon résultat, car la seconde partie du nombre n'est pas définie..
  
  Code:
  
  float dx=2.0 ;
  
  donnera ce qu'on souhaite réellement.
29/01/2011, 19h48
diogene
Citation:
D'autre part tu n'as pas mis les parenthèses permttant de définir les blocs

Code:

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

revient à faire :

Code:

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:

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

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:

float dx = 2 ;

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

Code:

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).
30/01/2011, 01h56
souviron34

pour la première remarque, exact, j'avas lu trop vite :oops::oops::oops: 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"...
30/01/2011, 10h00
diogene

Citation:

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:

float dx = 2.0f;

Mais on ne pouvait dire

Citation:

....ne donnera pas forcément le bon résultat, car la seconde partie du nombre n'est pas définie..
30/01/2011, 12h10
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..