function out=Subbytes(state)
H=[1 0 0 0 1 1 1 1;1 1 0 0 0 1 1 1;1 1 1 0 0 0 1 1;1 1 1 1 0 0 0 1;...
   1 1 1 1 1 0 0 0;0 1 1 1 1 1 0 0;0 0 1 1 1 1 1 0;0 0 0 1 1 1 1 1];
p=[1 1 0 0 0 1 1 0]';
for i=0:3
    for j=0:3
        (state(i,j))'=H*(state(i,j))'+P;
    end
end