JuliaDiffEq Vern9() solver taking too long

I think you may be measuring the compilation time of the julia code as well. There is a package called BechmarkTools.jl which runs code multiple times to get more accurate data for a benchmark.

You are summing up the time for 1000 runs(=compilation time + 1000 runs). Is fortran doing the same?