Правильно ли реализован метод вычисления значения интеграла методом трапеций ?

#include <stdio.h>
int main()
{
   float a, b;   // границы отрезка
   float n;    // приращение аргумента
   float s = 0;     // ~ значение интеграла
   int k;       // количество интервалов
   float x;     // аргумент
   float y1,y2; // значение функции в начале и в конце  
   int i; // счетчик
    
   printf("Введите нижнюю границу отрезка -> ");
   scanf("%f", &a);
   
   printf("Введите верхнюю границу отрезка -> ");
   scanf("%f", &b);
   
   printf("Введите приращение аргумента -> ");
   scanf("%f", &n);
   
   x = a;
   
   k = (b - a) / 2*n;  // считаем количество интервалов
   
   for (i = 1; i < k; i++) 
   {
      y1 = (x*x) / (x + 3);   // значение функции в начале интеграла
      x += n;                 // приращение аргумента
      y2 = (x*x) / (x + 3);   // значение функции в конце итеграла
      s += (y1 + y2) * n / 3;
   }
   printf("Значение интеграла: %8.6f", s);
 
return 0;
}