code: def get_jac(self):
          
          if self._jac is True:
              if self.sparse is True:
                  self._jac, self._colptrs, self._rowvals = self.be.sparse_jacobian_csc(self.exprs, self.dep)
              elif self.band is not None:  # Banded
                  self._jac = self.be.banded_jacobian(self.exprs, self.dep, *self.band)
              else:
                  f = self.be.Matrix(1, self.ny, self.exprs)
                  self._jac = f.jacobian(self.be.Matrix(1, self.ny, self.dep))
          elif self._jac is False:
              return False
  
          return self._jac
docs: Derives the jacobian from ``self.exprs`` and ``self.dep``.