Bugfix on initialization of jacobian matrix

we've been saved only because this->n_reactions() > this->n_species()