code: def incidence(self):
          
          self.C = \
              spmatrix(self.u, range(self.n), self.a1, (self.n, self.nb), 'd') -\
              spmatrix(self.u, range(self.n), self.a2, (self.n, self.nb), 'd')
docs: Build incidence matrix into self.C