Hello, the Jacobian is optional, so you can safely omit it. However, specifying the Jacobian (or at least its structure) can speed up the integration of the odepack solvers, e.g. for 1D, 2D and 3D systems. The Runge-Kutta solvers don't use the Jacobian. More, see help files ?lsoda ?lsode etc. Thomas