#NEWLINE# def outermost_radius( self , rho ) : #NEWLINE# #INDENT# '\n        Find the outermost radius at which the density of the envelope has\n        fallen to `rho`.\n\n        Parameters\n        ----------\n        rho : float\n            The density for which to determine the radius\n\n        Returns\n        -------\n        r : float\n            The radius at which the density has fallen to ``rho``\n        ' #NEWLINE# return ( self . r_0 * ( ( rho / rho . rho_0 ) ** ( 1.0 / self . power ) ) )