Write a function integrate_trap(f,a,b,N) taking the arguments f (a function), 'a' (a number), 'b' (a number) and 'N' (a positive integer), which estimates the integral of a generic function f on the closed interval [a,b] using the composite trapezoid rule with a total of N points, and returns the value of the approximated integral.