function main()

    h = 1e-2:1e-2:1

    dsin_num = [sin(i)/i for i in h]
    dsin_exact = 1
    dexp_num = [(exp(i)-1)/i for i in h]
    dexp_exact = 1

    err_sin = abs.(dsin_num .- dsin_exact)
    err_exp = abs.(dexp_num .- dexp_exact)
    p1 = plot(h,err_sin,label="Sin Error")
    plot!(p1,h,h,label="Linear")
    plot!(p1,h,h.^2,label="Quadratic")
    plot!(p1,xscale=:log10, yscale=:log10, minorgrid=true)
    p2 = plot(h,err_exp,label="Exp. Error")
    plot!(p2,h,h,label="Linear")
    plot!(p2,h,h.^2,label="Quadratic")
    plot!(p2,xscale=:log10, yscale=:log10, minorgrid=true)

    p = plot(p1,p2, layout=(1,2))
    return p
end

main()