Duck typing and consistent return types (Vector{typeof(y0)}) for all solvers. Changed order of arguments to F,y0,tspan. Added keyword jacobian for ode4s*. Tests adjusted accordingly.
Fixed deprecation warning for + in fdjacobian and ode4.