TODO for Optim.jl

more tests for curve_fit()
more documentation

Things to consider

tweak a lighterweight package that does not have Distribution dependency?
simple linear fitting?
use DualNumbers for finite difference (like in Optim.jl)?
add algorithms besides Levenberg-Marquardt?