IRIS Clustering

Fix and complete the program below so that it clusters the iris data. Print the clustering accuracy (what percentage of the data is clustered correctly)
import numpy as np
import matplotlib.pyplot as plt
from sklearn import datasets
from sklearn.cluster import kmeans
iris = datasets.load_iris()
# iris
_, ax = plt.subplots()
scatter = ax.scatter(iris.data[:, 0], iris.data[:, 1], c=iris.target)
ax.set(xlabel=iris.feature_names[0], ylabel=iris.feature_names[1])
_ = ax.legend(
scatter.legend_elements()[0], iris.target_names, loc="lower right", title="Classes"
)
plt.scatter(iris.data[:, 0], iris.data[:, 1], c=iris.target)
# plt.scatter(centroids[:, 0], centroids[:, 1], c='black', s=200, alpha=0.5)
# plt.show()
