code: def _rho02alpha(self, rho0, Rs):
  
          
  
          theta_Rs = rho0 * (4 * Rs ** 2 * (1 + np.log(1. / 2.)))
          return theta_Rs
docs: convert rho0 to angle at Rs
  
          :param rho0:
          :param Rs:
          :return: