code: def dquadmon_from_lambda(lambdav):
      r
      ll = numpy.log(lambdav)
      ai = .194
      bi = .0936
      ci = 0.0474
      di = -4.21 * 10**-3.0
      ei = 1.23 * 10**-4.0
      ln_quad_moment = ai + bi*ll + ci*ll**2.0 + di*ll**3.0 + ei*ll**4.0
      return numpy.exp(ln_quad_moment) - 1
docs: r"""Return the quadrupole moment of a neutron star given its lambda
  
      We use the relations defined here. https://arxiv.org/pdf/1302.4499.pdf.
      Note that the convention we use is that:
  
      .. math::
  
          \mathrm{dquadmon} = \bar{Q} - 1.
  
      Where :math:`\bar{Q}` (dimensionless) is the reduced quadrupole moment.