Duck typing and consistent return types (Vector{typeof(y0)}) for all solvers; see #7. Tests adjusted accordingly.
Changed order of arguments to F,y0,tspan and use varargs for tspan.
Add keyword jacobian for ode4s*.