inc = 36
   ki = seq(1,18203, by=inc)

   for (i in 0:length(ki)) {

     Dat4 = Data4[ki[i+1]:ki[i+2],3:6]

     #CovMat = cov(Dat5)
     #write.table(Dat5, file="CovMatrix.csv", sep = ",", col.names = F,
     append = T)
     #e = eigen(CovMat)
     #EValue = CovMat$values
     #write.table(EValue, file = "EigenValues.csv", sep = ",", col.names =
     F, append = T)

     library(tidyverse)
     library(cluster)
     library(factoextra)
     library(gridExtra)

     name2 = "F:/"
     ext = ".pdf"
     path4 = paste(name2,i,ext)
     pdf(file=path4)

     library(fpc)
     m<- kmeansruns(Dat4,krange = 1:8,criterion="asw", plot=F)
     cust=m$bestk

     km = eclust(Dat4, FUNcluster = c("kmeans"), k = cust, graph = FALSE,
     hc_metric = "euclidean")
     f = fviz_cluster(km, geom = "point", ellipse.type = "norm",
     ellipse.level = 0.68) + ggtitle("TrainingSet")

     CenterClust = km$centers
     #Cente = c(k$centers)
     #print(Cente)
     write.table(CenterClust, file = "KMCC.csv", sep = ",", col.names = F,
     append = T)

     grid.arrange(f, nrow = 2)

     dev.off()

   }