x<-c(1, 1, 1.5, 2, 2)
y<-c(1, 2, 1.5, 1, 2)

plot(x, y, xlim = c(0.5, 2.5), ylim = c(0.5, 2.5), pch=c(15, 17, 1, 18, 19))