code: def _compute_nonlinear_magnitude_term(self, C, mag):
          
          return self._compute_linear_magnitude_term(C, mag) +\
              C["b3"] * ((mag - 7.0) ** 2.)
docs: Computes the non-linear magnitude term