Discussion :

[Isamaz]

Bonjour,

Je souhaite créer un programme pour résoudre des intégrales afin de tracer au final une courbe. Pour l'instant je me concentre sur la premiere integral qui se résume à l'estimation de cette somme:


Mon soucis est que ma donné psi évolue selon n mais aussi en fonction des valeurs qu'elle a pris pour n-1. Je suis un peu perdu en C pour savoir comment lui dire d'utilisé les valeurs calculées avant pour faire la valeur suivante. Dois-je passé par un tableau? peut-il me le remplir au fur et a mesure des données?
J'ai tenté d’écrire un programme le plus simple possible (en estimant toutes les constantes), mais pour l'instant je ne peux pas lancer l'application, car c'est certainement beaucoup trop simple....:

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
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define pi 3.1416
#define R 8.144621
#define T 298.15
#define F 96485
#define alpha 1
 
float k1,to,Eimp,dj;
int j, n;
 
int main()
{
    k1=10.0;to=5.5;Eimp=0.6;dj=0.001;n=500;
    psi=0.0;
    dt=1/n;
    for(i=1;i<=n;i=i++); 
    { 
        for(j=i-1;i<n;j+=dj);
        {
            psi+=(1+1.33*1/(sqrt(pi))*sqrt(dt)*(psi(j)*pow(n-j+1,(3/2))+2*pow(n-j,(3/2))+pow(n-j-1,(3/2)))*dj)/(exp((-alpha*F/(R*T)))+1.33*(1/sqrt(pi))*sqrt(dt));
    }
 
    printf("l'integral de psi =%3.3f",psi);
    return 0;
    }
 
 
}

Auriez-vous des conseils? J'ai bouquiné tout le tuto C en 20h, quelques idées me sont venu mais ça reste un peut flou pour mon cas...

Je vous remercie,

Isabelle