qtrap1m001.jl

@fastmath function quad_trap(f, a, b, n)
    h = (b - a) / n
    res = 0.5 * (f(a) + f(b))
    for i in 1:n-1
        res += f(a + i*h)
    end
    res * h
end


function g()
    a = 0.0
    for _ in 1:100000
        a = quad_trap(x -> exp(x) - 10.0, -1.0, 1.0, 10000)
    end
    a
end

println(g())
flush(stdout)
@timev g()