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