for i, h in enumerate([0.01, 0.1, 1, 5]):
    my_pdf = gaussian_kde(osservazioni, h)
axarr[i // 2, i % 2].plot(x, my_pdf(x), 'r')
axarr[i // 2, i % 2].set_title('Bandwidth: {0}'.format(h))
axarr[i // 2, i % 2].hist(osservazioni, normed=1, alpha=0.3)
Python Gaussian Kernel density calculate score for new values