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