function A=SymMatrix(n),

e=ones(n,1);
A=spdiags([-e 2*e -e],-1:1,n,n);
Q=eye(n)-2*e*e'/(e'*e);
A=Q*A*Q;

