<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_all_stability_map(self, max_r, increments=50, delu_dict=None, delu_default=0, plt=None, labels=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 particles of different polymorphs against its effect radius Args: max_r (float):
 The maximum radius of the particle to plot up to. increments (int):
 Number of plot points 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 plt (pylab):
 Plot labels (list):
 List of labels for each plot, corresponds to the list of se_analyzers 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. """


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

        for i, analyzer in enumerate(self.se_analyzers):
            label = labels[i] if labels else ""
            plt = self.plot_one_stability_map(analyzer, max_r, delu_dict,
                                              label=label, plt=plt,
                                              increments=increments,
                                              delu_default=delu_default,
                                              from_sphere_area=from_sphere_area,
                                              e_units=e_units, r_units=r_units,
                                              normalize=normalize,
                                              scale_per_atom=scale_per_atom)

        return plt