def solve_equilibrium_point(self, analyzer1, analyzer2,
                                delu_dict={}, delu_default=0, units="nanometers"):
        """
        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
        """

        # 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