code: def linear_response_matrix(self, kwargs_lens=None, kwargs_source=None, kwargs_lens_light=None, kwargs_ps=None):
          
          A = []
          for i in range(self._num_bands):
              if self._compute_bool[i] is True:
                  A_i = self._imageModel_list[i].linear_response_matrix(kwargs_lens, kwargs_source, kwargs_lens_light, kwargs_ps)
                  if A == []:
                      A = A_i
                  else:
                      A = np.append(A, A_i, axis=1)
          return A
docs: computes the linear response matrix (m x n), with n beeing the data size and m being the coefficients
  
          :param kwargs_lens:
          :param kwargs_source:
          :param kwargs_lens_light:
          :param kwargs_ps:
          :return: