added numerical evaluation of Jacobian matrices using numdifftools (untested)