построить среднее значение и стандартное значение после зацикливания по нескольким плотностям

Question

У меня есть эта функция

def plot_density(X, bandwidth, ax):

    x_axis = np.linspace(0, 1, 101)
    dens_list = []
    for x in X:
        kernel = stats.norm(x, bandwidth).pdf(x_axis)
        dens_list.append(kernel)

    ax.set_xlim(0,1)

    density = np.mean(density_list, axis = 0)

    ax.plot(x_axis, density ,color = 'red')
    sns.rugplot(X)

Я хочу l oop for i in range (4):, затем получить среднее и стандартное значение повторяющихся плотностей и построить среднее значение и стандартное значение, как я могу это сделать