#include "stdio.h"
#include "math.h"

float f(float par)
{
return 4*par*par*par - 6*par -12;   //4*x^3 - 6*x -12
}


int main()
{

float a,b,erg=0,n,h;
int i;
char nochmal;

printf("simpson-näherung");

printf("\n\n\n\nlinke integrationsgrenze: ");
scanf("%f",&a);
printf("\nrechte integrationsgrenze: ");
scanf("%f",&b);
printf("\niterationszahl: ");
scanf("%f",&n);

h=(b-a)/2/n;

erg += f(a);   //f_0

for(i=1; i<=2*n; i+=2)  
{
erg += 4*f(a+i*h);   //f_1, f_3, f_5, ...., f_2n-1
}

for(i=2; i<2*n; i+=2)    
{
erg += 2*f(a+i*h);	//f_2, f_4, ...., f_2n-2
}

erg += f(a+2*n*h);	//f_2n

erg = erg*h/3;

printf("\n\nlösung: %f",erg);


}