ダメ出し：par(new=TRUE) は使わない

# P 君の最初の位置と、A さん B さんそれぞれの位置
P <- c(0, 0)
A <- c(0, 10)
B <- c(6, 8)

library(gplots)
niter <- 1000
a <- 1/2
b <- 1/2
pos <- matrix(NA, nr=niter, nc=2)
pos[1, ] <- a * (A - P)
for(i in 2:niter){
    if(odd(i)){
        pos[i, ] <- (pos[i - 1,] - P) + a * (A - pos[i - 1, ])
    }else{
        pos[i, ] <- (pos[i - 1,] - P) + b * (B - pos[i - 1, ])
    }
}
pos <- rbind(P, pos)

tri <- rbind(P, A, B)
plot(tri, xlim=range(tri[, 1]), ylim=range(tri[, 2]), xlab="", ylab="")
# par(new=TRUE)
# plot(pos, type="l", xlim=range(tri[,1]), ylim=range(tri[,2]), xlab="", ylab="")
lines(pos)
abline(0, pos[niter - 1, 2] / pos[niter - 1, 1], lty=2)
abline(0, pos[niter, 2] / pos[niter,1], lty=2)
polygon(tri[, 1], tri[, 2])