from sklearn.cluster import KMeans
model = KMeans(n_clusters=40, init='k-means++', max_iter=300, n_init=50, random_state=3)
model.fit(data)