plist = RecurrenceTable[{p[n + 1] == 4 p[n] - p[n - 1], p[1] == 2,
p[0] == 1}, p, {n, 1, 10}]
qlist = RecurrenceTable[{q[n + 1] == 4 q[n] - q[n - 1], q[1] == 1,
q[0] == 0}, q, {n, 1, 10}]
Thread[{#1, #2, #3} & @@ {Range[10], plist, qlist}]
RecurrenceTable[{p[n + 1] == 2 p[n] + 3 q[n],
q[n + 1] == p[n] + 2 q[n], p[1] == 2, q[1] == 1},
{p, q}, {n, 1, 10}] |