code: def dim_dm(self, pars):
          r
          self._set_parameters(pars)
          num1 = self.otc * np.sin(self.ang)
          result = -self.sigmai * num1 / self.denom
  
          return result
docs: r"""
          :math:Add formula