x <- c(rep(-1,1000),
seq(-1.5,-1,length.out=500),
seq(0,1,length.out=1000),
rep(0, 500),
rep(.75, 500),
rep(1, 500),
seq(0,1,length.out=1000))
y <- c(seq(-1,1, length.out=1000),
rep(-1,500),
rep(0,1000),
seq(-1,1,length.out=500),
seq(-1, 0, length.out=500),
seq(0, 1, length.out=500),
rep(1,1000))
jr <- rbind(x,y)
plot(y~x, xlim=c(-3,2), ylim=c(-3,3))