def plot_posterior_distribution(ax, N_lam, N_turb, obj_to_rasterize, simplified=False):
    a = (N_lam + (1.0 / 2))
    b = (N_turb + (1.0 / 2))
    ax.plot(x, beta.pdf(x, a, b), label=(((('$N = ' + str((N_lam + N_turb))) + '$\n$l = ') + str(N_lam)) + '$'), linewidth=(4 if simplified else 2))
    ax.plot([beta.mean(a, b)], [beta.pdf(beta.mean(a, b), a, b)], 'ro', markersize=(12 if simplified else 8))
    (_, y_max) = ax.get_ylim()
    ax.set_xticks([0.0, 0.2, 0.4, 0.6, 0.8, 1.0])
    if simplified:
        ax.set_xticklabels([])
        ax.set_yticklabels([])
    else:
        obj = ax.fill_between([beta.ppf(0.1, a, b), beta.ppf(0.9, a, b)], [0.0, 0.0], [y_max, y_max], alpha=0.3)
        obj_to_rasterize.append(obj)
        ax.legend(fontsize=12)
        ax.set_xticks([0.0, 0.2, 0.4, 0.6, 0.8, 1.0])
        ax.set_xlabel('$p$')
def plot_posterior_distribution(ax, N_lam, N_turb, obj_to_rasterize, simplified=False):
    a = (N_lam + (1.0 / 2))
    b = (N_turb + (1.0 / 2))
    ax.plot(x, beta.pdf(x, a, b), label=(((('$N = ' + str((N_lam + N_turb))) + '$\n$l = ') + str(N_lam)) + '$'), linewidth=(4 if simplified else 2))
    ax.plot([beta.mean(a, b)], [beta.pdf(beta.mean(a, b), a, b)], 'ro', markersize=(12 if simplified else 8))
    (_, y_max) = ax.get_ylim()
    ax.set_xticks([0.0, 0.2, 0.4, 0.6, 0.8, 1.0])
    if simplified:
        ax.set_xticklabels([])
        ax.set_yticklabels([])
    else:
        obj = ax.fill_between([beta.ppf(0.1, a, b), beta.ppf(0.9, a, b)], [0.0, 0.0], [y_max, y_max], alpha=0.3)
        obj_to_rasterize.append(obj)
        ax.legend(fontsize=12)
        ax.set_xticks([0.0, 0.2, 0.4, 0.6, 0.8, 1.0])
        ax.set_xlabel('$p$')