How to write many interaction terms in linear regression model in R?

Question

I have a linear model with 15 independent variables. I want to run the model with all pairwise interaction terms. y~a+b+c+...15 independent variables+ab+bc+...1514