probleme avec simpson

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

merci beaucoup sa fait plus de 2h que je cherche le pb et je n'est pas trouvé merci encors