void Rossler::jac(const_it_t &state, state_t & out, double time) { assert(out.size() == dim*dim); auto it = out.begin(); it[0] = 0; it[1] = -1; it[2] = -1; it += dim; it[0] = 1; it[1] = a; it[2] = 0; it += dim; it[0] = state[2]; it[1] = 0; it[2] = state[0]-c; }