def METHOD_NAME(microstates, lifetimes, out, ax_prop=None, ax_distrib=None):
    states = np.unique(microstates)
    # Plot
    if ax_prop is None and ax_distrib is None:
        fig, axes = plt.subplots(ncols=2)
        ax_prop = axes[0]
        ax_distrib = axes[1]
    else:
        fig = None
    for s in states:
        ax_prop.bar(s, out[str(s) + "_Proportion"])
        ax_distrib.plot(lifetimes[s], label=str(s))
    plt.legend()
    ax_prop.set_xticks(range(len(states)))
    ax_prop.set_title("Proportion")
    ax_distrib.set_title("Lifetime Distribution")
    return fig