v0.1.8 Optimized Jacobian approximation + complete docstring documentation

Highlights of this relase:

• set central differences (3-point mode) as default jacobian approximation jacmode=3
• automatic jaceps based on machine epsilon if jaceps=None
• complete docstring documentation for all functions and classes