function y = lin2d(x) % % eigenvalues of A are -0.8 and +0.8 % A = [ -0.2240 -0.7680 -0.7680 0.2240 ]; b = [1;1]; y = A*x + b;