function x = trap(f, a, b, n)
    x = (b-a)/n * sum(f(a:(b-a)/n:b)) - (b-a)/(2*n) * (f(a) + f(b));