function v=perm(n)
v=1;
for i=2:n; 
    j=ceil(i*rand);
    v=[v j];
    v([j i])=[i v(j)];
end
