code: def dim_dc(self, pars):
          r
          self._set_parameters(pars)
          # term 1
          num1a = self.m * np.sin(self.ang) * np.log(self.w * self.tau)\
              * self.otc
          num1b = self.m * self.otc * np.pi / 2 * np.cos(np.pi / 2)
          term1 = self.sigma0 * (-num1a - num1b) / self.denom
  
          # term 2
          num2a = -self.m * self.otc * np.cos(self.ang)
          num2b = -2 * np.log(self.w * self.tau) * self.otc * np.cos(self.ang)
          num2c = 2 * self.otc * np.pi / 2 * np.cos(self.ang)
          num2d = 2 * np.log(self.w * self.tau) * self.otc2
          numerator = num2a * (num2b + num2c) + num2d
          term2 = self.sigma0 * numerator / (self.denom ** 2)
  
          result = term1 + term2
  
          return result
docs: r"""
          :math:Add formula