code: def linear_response_matrix(self, kwargs_lens=None, kwargs_source=None, kwargs_lens_light=None, kwargs_ps=None):
          
          A = self._response_matrix(self.ImageNumerics.ra_grid_ray_shooting, self.ImageNumerics.dec_grid_ray_shooting,
                                        kwargs_lens, kwargs_source, kwargs_lens_light, kwargs_ps, self.ImageNumerics.mask)
          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: