using python without using any powerful functions other than math module and copy function

Write a function defined as:def Simpson(fcn, a, b, npoints = 21):

Purpose: use Simpsons 1/3 rule to estimate the integral of fcn(x), between the limits of a and b.

fcn: the function we want to integrate

a and b : the lower and upper limits of integration

npoints : The number of integration points used in the range a to b (inclusive). npoints must be an ODD number.

If npoints is not ODD, then add 1 to make it odd!

return value: the estimate of the integral

Write and call a main() function that uses your Simpson function to estimate and print the integral of:

: with a=1, b=3, npoints=10

with a=2, b=3, npoints= 23

2 3cos(1) cos(2.c) *23