function [keyTable, wmMatrix] = mylsb(srcMatrix)
[ms, ns] = size(srcMatrix);
wmLen = ms * ns;
randomWM = randint(1,wmLen,[0,1]);
%randomWM = zeros(1, wmLen);
%for i = 1:wmLen
modifyQ = [1,2,4,8];
keyTable = zeros(ms, ns);
wmMatrix = srcMatrix;
for i = 1:ms
for j = 1:ns
index = floor(rand*4) + 1;
if (randomWM(1, (i-1)*ms+j) == 1) & ( bitand( srcMatrix(i,j), index) == 0 )
wmMatrix(i,j) = srcMatrix(i,j) + modifyQ(index);
end
if (randomWM(1, (i-1)*ms+j) == 0) & ( bitand( srcMatrix(i,j), index) ~= 0 )
wmMatrix(i,j) = srcMatrix(i,j) - modifyQ(index);
end
keyTable(i,j) = index;
end
end