void(* evalJ)(const C_INT *, const double *, const double *, const C_INT *, const C_INT *, double *, const C_INT *)
C_INT operator()(evalF f, C_INT *neq, double *y, double *t, double *tout, C_INT *itol, double *rtol, double *atol, C_INT *itask, C_INT *istate, C_INT *iopt, double *rwork, C_INT *lrw, C_INT *iwork, C_INT *liw, evalJ jac, C_INT *jt)