Найти определенный интеграл методом Трапеций

#define _USE_MATH_DEFINES
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
// интегрируемая функция:
double f(double x)
{
    return sin(x*exp(-x*x));
}
double Integral(double ax,double bx,double eps)
{
    int k=100; // начальное к-во разбиений
    // значения интегралов на текущей и предыдущей интерациях:
    double val1=0,val2=0.; 
    do
    {
        val1=val2; // предыдущее значение равно текущему
        double s=0; // переменная для суммирования 
        double x=ax; // начинаем с ax
        double hx=(bx-ax)/k; // шаг интегрирования по x
        for(int i=0;i<k;i++)// повторяем к раз
        {
            double t=hx*(f(x)+f(x+hx))/2.;
            s+=t; // суммируем площадь очередной трапеции
            x+=hx;
        }
        val2 = s; // вычисленное новое значение
        k*=2; // удвоение количества разбиений
    }while (fabs(val1-val2)>eps); // повторяем до достижения точности eps
    return val2; // возвращаем результат интегрирования
}
int main() 
{
    //system("chcp 1251 > NUL");
    double eps=1e-3; // точность интегрирования
    double xb=0; // начало интервала
    double xe=M_PI; // конец интервала
    double res=Integral(xb,xe,eps); // вычисление нового значения
    printf("Вычисленное значение=%lf\n",res);
    //getchar();
    return 0;
}