d  <- read.table("inp_iris.txt",header=TRUE)
x  <- d[,1:4] 
km <- kmeans(x,4,20) #k-means: number of cluster is 4
km$center
cbind(x,km$cluster)
library(cluster)
clusplot(x, km$cluster, color=TRUE, shade=TRUE, labels=4, lines=0)