function [ A ] = SpSymmRand( n, tau )
%
nel = 0;
R = [1:tau];
C = [1:tau];
V = [1:tau];
for j=1:n
    R(j)=j;
    C(j)=j;
    V(j)=rand(1,1);
end
for j=n+1:2:tau
    Rj = min(int32(1+rand(1,1)*n),n);
    Cj = min(int32(1+rand(1,1)*n),n);
    Vj = rand(1,1);
    R(j)=Rj;
    C(j)=Cj;
    V(j)=Vj;
    R(j+1)=Cj;
    C(j+1)=Rj;
    V(j+1)=Vj;
end
A=sparse(R,C,V,n,n);
%spy(A);pause
end

