def METHOD_NAME(sds, scales, slope, intercept, ax=None):
    if ax is None:
        fig, ax = plt.subplots()
        fig.suptitle(
            "Standardized Dispersion as a function of Sample-Bin size"
            + ", slope = "
            + str(np.round(slope, 2))
        )
    else:
        fig = None
        ax.set_title(
            "Standardized Dispersion as a function of Sample-Bin size"
            + ", slope = "
            + str(np.round(slope, 2))
        )
    ax.set_ylabel(r"$\log_{10}$(Standardized Dispersion)")
    ax.set_xlabel(r"$\log_{10}$(Bin Size)")
    ax.scatter(np.log10(scales), np.log10(sds), marker="o", zorder=2)
    fit_values = [slope * i + intercept for i in np.log10(scales)]
    ax.plot(
        np.log10(scales),
        fit_values,
        color="#FF9800",
        zorder=1,
        label="Fractal Dimension = " + str(np.round(1 - slope, 2)),
    )
    ax.legend(loc="lower right")