kill(all)$

van=facebook.com/LearningCalculus$

f(n):=(2+%i)^n+(2-%i)^n$

realpart(f(n));

imagpart(f(n));

L:[]$

for k:0 thru 6 step 1 do L:append(L,[[k,expand(f(k))]]);

L;