code: def solve_equilibrium_point(self, analyzer1, analyzer2,
                                  delu_dict={}, delu_default=0, units="nanometers"):
          
  
          # Set up
          wulff1 = analyzer1.wulff_from_chempot(delu_dict=delu_dict,
                                                delu_default=delu_default,
                                                symprec=self.symprec)
          wulff2 = analyzer2.wulff_from_chempot(delu_dict=delu_dict,
                                                delu_default=delu_default,
                                                symprec=self.symprec)
  
          # Now calculate r
          delta_gamma = wulff1.weighted_surface_energy - wulff2.weighted_surface_energy
          delta_E = self.bulk_gform(analyzer1.ucell_entry) - self.bulk_gform(analyzer2.ucell_entry)
          r = ((-3 * delta_gamma) / (delta_E))
  
          return r / 10 if units == "nanometers" else r
docs: Gives the radial size of two particles where equilibrium is reached
              between both particles. NOTE: the solution here is not the same
              as the solution visualized in the plot because solving for r
              requires that both the total surface area and volume of the
              particles are functions of r.
  
          Args:
              analyzer1 (SurfaceEnergyPlotter): Analyzer associated with the
                  first polymorph
              analyzer2 (SurfaceEnergyPlotter): Analyzer associated with the
                  second polymorph
              delu_dict (Dict): Dictionary of the chemical potentials to be set as
                  constant. Note the key should be a sympy Symbol object of the
                  format: Symbol("delu_el") where el is the name of the element.
              delu_default (float): Default value for all unset chemical potentials
              units (str): Can be nanometers or Angstrom
  
          Returns:
              Particle radius in nm