#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)));
}