Hi all, I was wondering why a gsl_odeiv_system expressly needs a jacobian matrix to be defined, while only the bsimp solver appears to be using it. The current formulation makes it sound as if the gsl solvers cannot be used if no jacobian is provided. Greetings, Fleur Kelpin