// Evaluate the integral with the trapezium method
// input: range of integration and number of subinterval 
// output: the integral
// the fuction to be itegrated is declared in myFunction.h
// e defined in  myFunction.cpp

double TrapezoidalArea(double, double, int);