code: def _compute_magnitude(self, rup, C):
          
          return C['b1'] + (C['b2'] * rup.mag) + (C['b3'] * (rup.mag ** 2))
docs: Compute the first term of the equation described on p. 199:
  
          ``b1 + b2 * M + b3 * M**2``