x1 <- c(4,-1,3)
x2 <- c(1,3,5)
origin <- c( 0, 0, 0)
m <- t(cbind(x1,x2,origin))
library("scatterplot3d") # load
scatterplot3d(m,
type="p",
color = c("red","blue","black"),
pch = 1,
xlim = c(-6,6),
ylim = c(-6,6),
zlim = c(-6,6))

I'm afraid I don't know what a triplot is but this is a scatterplot3d