function to define a polynomial returns incorrect value when called

Question

I am writing a program to integrate a polynomial using the trapezium rule (https://en.wikipedia.org/wiki/Trapezoidal_rule). Although the program does run, it is not returnin