function Y  = shadock( X )

U=rand(1);
if (X==1) then Y=(U<5/6) + 2*(U>5/6 & U<11/12) + 3*(U>11/12);
    else if (X==2) then Y=(U<1/4) + 2*(U>1/4 & U<3/4) + 3*(U>3/4);
        else Y=(U<1/4) + 3* (U>1/4);
        end
    end

endfunction
