docs: def dim_dsigmai(self, pars):
          r"""
          :math:Add formula
          """
          self._set_parameters(pars)
          result = np.sum(- self.m * self.otc * np.sin(self.ang) / self.denom,
                          axis=1)
  
          return result