В R как бы представлять полиномиальные выражения и делать полиномиальную математику с числовыми векторными объектами? Например:

x1 <- c(2,1)  # 2 + x
x2 <- c(-1,3)  # -1 + 3*x

И хотите:

x1 * x2 # to return -2 + 5*x + 3*x^2 

Примечание: Сегодня утром я ответил на вопрос, а затем, по-видимому, постер удалил его (заставив меня задуматься, не было ли это домашнее задание). Поэтому я перекладываю вопрос по памяти.