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

double a=0, b=1,c=1.5,d=3;

void Sims(float);
float y(float);
int main()
{
	printf("Metod Simpsona:\n");
    printf("h=0.1\n");
	Sims(50);
    printf("h=0.2\n");
	Sims(25);
    printf("h=0.5\n");
	Sims(10);
	getch();
}

void Sims(float n)
{
  float E=0.01,h;
  float sum=y(a)-y(b);  
  h=(b-a)/(2*n);
  for(int i=1;i<=(2*n-1);i++)
  {
    sum=sum+(3+E)*y(a+(i*h));
    E=(-1)*E;
  }
  sum=h/3*sum;  
  printf ("Rezul\'tat = %f\n",sum);
}

float y(float x)
{
  return (((c*x)/2)+(cos(d*x)));
}