def nu_max(self, *args):
return <NUM_LIT>* (self.mass(*args) /<EOL>(self.radius(*args)**<NUM_LIT:2> * np.sqrt(self.Teff(*args)/<NUM_LIT>)))<EOL>
Returns asteroseismic nu_max in uHz

        reference: https://arxiv.org/pdf/1312.3853v1.pdf, Eq (3)