A = [1 -1; 1 3];



k = 10;



u0 = [1 0]';

uk = A^k*u0;

uk2 = 2^k*[1-k/2; k/2];





Ak2 = 2^(k-1)*[2-k -k; k 2+k];