Discussion :

[herosion]

Bonjour à tous,

Je souhaite programmer le modèle de Onsager. Pour cela on me donne ces infos:
voir fichier joint

voila ce que j'ai fais:

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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
#include<stdio.h>
#include<math.h>
 
#define pie 3.1419
 
 
float k1(float,float);
float integK1(float,float,float,float);
float e1(float,float);
float integE1(float,float,float,float);
 
 
int main()
{
   float km,q,K,K1q,E1q,CK,t;
 
   float inf = 0;
   float sup = pie/2;
   float N = 0.5;
   float pas = (sup-inf)/N;
 
 
  for(km=0.5; km<=5; km+=0.1)
  {
      K = 1./km; //km = 1/k
      q=(2*sinh(2*K))/(cosh(2*K)*cosh(2*K)); //q(K)
      K1q = integK1(pas,inf,sup,q); //E1(q)
      E1q = integE1(pas,inf,sup,q); //K1(q)
      //CK = ; c'est ce que je veux au finale
 
  }
printf("%f %f\n",K1q,E1q);
 
 
  return 0;
}
 
 
float k1(float q, float t)
{
 
      return 1/(sqrt(1-(q*q)*(sin(t)*sin(t))));//ma fonction de K1(q)
 
}
 
float integK1(float pas,float inf, float sup,float q)//integrale de k1
{
     int i=0;
     float K1=0;
 
     for(i=inf; i<sup; i+=pas)
     {
          K1 = K1 + (k1(q,i)*pas);
     }
 
     //printf("%f",K1);
     return K1;
 
}
 
float e1(float q, float t)
{
      return sqrt(1-(q*q)*(sin(t)*sin(t))); //fionction de E1(q)
}
 
float integE1(float pas,float inf, float sup, float q)//integrale de ma fonction
{
     int i=0;
     float E1=0;
 
     for(i=inf; i<sup; i+=pas)
     {
          E1 = E1 + (e1(q,i)*pas);
     }
     return E1;
}

les intégrales que j'ai calculés (K1q et E1q) ne marchent pas.
Pourriez-vous m'aider svp

[herosion]

personne pour m'aider?

J'aimerais juste savoir pourquoi mes intégrales E1 et K1 marchent pas.
en gros:
- j'ai écris ma fonction k1 (float k1(float q, float t) // t=phi
- puis je calcul son intégrale (integK1): j'utilise la méthode des rectangles pour cela.

Je fais la même chose pour E1.

Cordialement
Ben

[jeroman]

Le problème, c'est que tu nous donnes le scan d'une feuille raturée de partout. On ne sait pas ce que c'est, on ne sait pas non plus quelle partie est bonne ou mauvaise.

Donne-nous les équations en question, au propre.

[herosion]

Voila j'ai réécris au propre les équations, j’espère que ça vous aidera.

Cordialement
Ben

[diogene]

Code :

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

1
2
   float N = 0.5;
   float pas = (sup-inf)/N;

La valeur de N est bizarre !

[herosion]

effectivement , j'ai confondu N et la valeur du pas au final. lorsque j'affiche les valeurs de mon intégrale K1, elles restent les mêmes, de km=0.5 a km<5, c'est pas normale.