function A = myexpmatrix(n,r)
u = (1-2*r)*ones(n,1);
v = r*ones(n-1,1);
A = diag(u) + diag(v,1) + diag(v,-1);