Probleme avec simpson

Résolu
zeroduito Messages postés 2 Date d'inscription mercredi 12 novembre 2008 Statut Membre Dernière intervention 12 novembre 2008 - 12 nov. 2008 à 21:12
zeroduito Messages postés 2 Date d'inscription mercredi 12 novembre 2008 Statut Membre Dernière intervention 12 novembre 2008 - 12 nov. 2008 à 22:00
salut a tout le monde
j'ai un probleme avec calcule d'integrale (simpson).
j'ai essaye de créer une fonction pour calculer l'intergrale mais  le programe n'affiche rien

#include <stdlib.h>
#include <stdio.h>
#include <math.h>

 double f(double x)
 {
 return exp(-x*x/2.0);
 }

double simpson(double n, double a, double b)
 {
 double  h, res, som1, som2,  i;

 h = (b - a) / n;

 som1=0.0;
 for(i=1;i<=(n/2)-1 ;i+2)
                        {

                         som1=som1+ f(a+(2*i*h));

                        }

 som2=0.0;
 for(i=1;i<=(n/2) ;i+2)
                        {
                         som2=som2+ f(a+(2*(i-1)*h));

                        }

res=(h/3.0)*(f(a)+f(b)+(2*som1)+(4*som2));

 return (res);
 }

int main()
{
double n=8.0, a=8.0, b=54.0, integ;
 simpson(n,a,b);
integ=simpson(n,a,b);
printf("\n %f",integ);

}
merci d'anavce
A voir également:

2 réponses

cs_darunia Messages postés 354 Date d'inscription mercredi 18 décembre 2002 Statut Membre Dernière intervention 24 mars 2011 2
12 nov. 2008 à 21:50
J'ai pas regardé en détail, mais 
for(i= 1;i<=(n/2)-1 ;i+2) ça risque pas de t'incrémenter i

 for(i =1;i<=(n/2)-1 ;i=i+2) serait mieux

D@runia
3
zeroduito Messages postés 2 Date d'inscription mercredi 12 novembre 2008 Statut Membre Dernière intervention 12 novembre 2008
12 nov. 2008 à 22:00
merci beaucoup sa fait plus de 2h que je cherche le pb et je n'est pas trouvé merci encors
0
Rejoignez-nous