function RK=roundKeyiterations (X,P, N,r,Sbox,m)
% DIfferent size of input block
%m=1;
RK=zeros(r,4*m);
for itr=1:r
    for itm=1:m
   y=skewtententierperturbed(X(itr),P(itr),N)
  %Y(i)=double(typecast(y,'uint8'))
  op=double(typecast(uint32(y),'uint8'))
    opd(1)=bitxor(bitxor(op(1),op(2)),op(3));
    opd(2)=bitxor(bitxor(op(1),op(2)),op(4));
    opd(3)=bitxor(bitxor(op(1),op(4)),op(3));
    opd(4)=bitxor(bitxor(op(4),op(2)),op(3));
    opd=Sbox(opd+1);
  RK(itr,(itm-1)*4+1:itm*4)=opd(:);
   X(itr)=y;
    end
end