qtrap1m002.jl

func(x) = @inline exp(x) - 10.0;


@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:1000000
        a = quad_trap(func, -1.0, 1.0, 10000)
    end
    a
end

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