#include <iostream>
using std::cout;
using std::cin;
using std::endl;

#include "TrapeziodalIntegration.h"
int main()
{

  cout << "Calculation of an integral with the Trapezium method" << endl;
  double startX, endX ;
  cout << "Low value of the integration interval: ";
  cin >> startX ;
  cout<< "High value of the integration interval: ";
  cin >> endX ;
  cout << "Number of sub-intervals: ";
  int nSubIntervals;
  cin >> nSubIntervals;

  double area = TrapezoidalArea(startX, endX, nSubIntervals);
    
  cout << "The integral of the function in the interval "<<startX<<" "<<endX<<"is :" << area <<endl;

  return 0;

}