WikiのようなProjective Varietyのような絵(こちら)が書きたかったが… my.pr <- function(xy,x0,y0,y1){ k <- 1/(xy[,2]-y1+1) new.y <- y0 * k new.x <- x0 + (xy[,1]-x0) * k return(cbind(new.x,new.y)) } x <- y <- seq(from=-10,to=10,length=1000) xyz…