set logscale x
set logscale y
set format y "%.e"
set format x "%.e"
plot [40:] [1e-15:] "rk4.dat" using 1:2 pt 7 ti "taylor2",\
    "" using 1:3 pt 11 ps 1.2 ti "trapezoid",\
    "" using 1:4 pt 9 ps 1.2 ti "rk2",\
    "" using 1:5 pt 5 ti "rk4",\
    1/x**2,1/x**4