calculate_silhouettes function

Calculate silhouette scores