Area under the solution curve

[rafael.guerra]

I think you can use QuadGK as follows:

using QuadGK
julia> AUC, err = quadgk(sol, 0, 1.0)
(0.8641589183483787, 4.908429218630772e-11)