Get WCSS (Within-Cluster Sum of Square) for optimal number of clusters in kmeans

I want to use kmeans for a clustering project, and I need to get the optimal number of clusters with the elbow method, is there a way to get the WCSS in the implementation of Clustering.jl to build something like this

the dataset could be like the typical iris

using RDatasets, Clustering, Plots
iris = dataset("datasets", "iris"); # load the data