label = kmeans.labels_
sil_coeff = silhouette_score(X, label, metric='euclidean')
print('For n_clusters={}, The Silhouette Coefficient is {}'.format(n_cluster,
    sil_coeff))
Scikit Learn - K-Means - Elbow - criterion