<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def plot_one_stability_map(self, analyzer, max_r, delu_dict=None, label="", increments=50, delu_default=0, plt=None, from_sphere_area=False, e_units="keV", r_units="nanometers", normalize=False, scale_per_atom=False):
 """ Returns the plot of the formation energy of a particle against its effect radius Args: analyzer (SurfaceEnergyPlotter):
 Analyzer associated with the first polymorph max_r (float):
 The maximum radius of the particle to plot up to. 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. label (str):
 Label of the plot for legend increments (int):
 Number of plot points delu_default (float):
 Default value for all unset chemical potentials plt (pylab):
 Plot from_sphere_area (bool):
 There are two ways to calculate the bulk formation energy. Either by treating the volume and thus surface area of the particle as a perfect sphere, or as a Wulff shape. r_units (str):
 Can be nanometers or Angstrom e_units (str):
 Can be keV or eV normalize (str):
 Whether or not to normalize energy by volume """


        plt = plt if plt else pretty_plot(width=8, height=7)

        wulffshape = analyzer.wulff_from_chempot(delu_dict=delu_dict,
                                                 delu_default=delu_default,
                                                 symprec=self.symprec)

        gform_list, r_list = [], []
        for r in np.linspace(1e-6, max_r, increments):
            gform, r = self.wulff_gform_and_r(wulffshape, analyzer.ucell_entry,
                                              r, from_sphere_area=from_sphere_area,
                                              r_units=r_units, e_units=e_units,
                                              normalize=normalize,
                                              scale_per_atom=scale_per_atom)
            gform_list.append(gform)
            r_list.append(r)

        ru = "nm" if r_units == "nanometers" else "\AA"
        plt.xlabel(r"Particle radius ($%s$)" %(ru))
        eu = "$%s/%s^3$" %(e_units, ru)
        plt.ylabel(r"$G_{form}$ (%s)" %(eu))

        plt.plot(r_list, gform_list, label=label)

        return plt