% 結婚相手を求める
% 1,2,\cdots,n の順列をランダムに作る

function p = matching (n)
    p = zeros(1,n);
    for k=n:(-1):1
        zeroindex = find(p == 0);
        l = ceil(k*rand);
        p(zeroindex(l)) = k;
    end
end
