Optimize performance

[mjungo]

Performance of solve_ivp only marginally better than basic finite differences

• understand and optimize hyperparameters of solve_ivp
• explore Numba
• explore JitCode
• np.einsum instead of np.matmul?
• implement Jacobian (issue Implement analytical Jacobian #14)
• ...