Simulating a LTI system in R

Question

I want to simulate a LTI system with input

    X[n] = delta[n] + delta [n+1]

And the system is this:

    y[n] = y[n-1] + 2x[n] + 3