def METHOD_NAME(
    savefolder,
    varvalues,
    status=None,
    label1="var1",
    label2="var2",
    xlab="variable",
    ylab="Normalized frequency",
    linewidth=2,
    framealpha=0.3,
    markersize=5,
    fontsize=18,
):
    det_mean, det_std = det_multi_ensemble(savefolder, varvalues, status)
    figure()
    grid("on")
    plot(
        varvalues,
        det_mean,
        "b-o",
        linewidth=linewidth,
        markersize=markersize,
        label=label1,
    )
    plot(
        varvalues,
        det_std,
        "r-o",
        linewidth=linewidth,
        markersize=markersize,
        label=label2,
    )
    xlabel(xlab, fontsize=fontsize)
    ylabel(ylab, fontsize=fontsize)
    #    xlim(0,xmax)
    tick_params(axis="both", which="major", labelsize=fontsize)
    legend(fancybox=True, framealpha=framealpha, loc="best", fontsize=fontsize)
    return