function [out]=GF8_multiply(a,b);
c=double(a);
d=uint16(0);
for i=1:8
    d=bitxor(uint16(mod(floor((c/(2^(i-1)))),2))*b*2^(i-1),d); %輸出乘法結(jié)果,大于8位
end
f=double(d);
for j=1:7                         %輸出乘法結(jié)果,
    mius(j)=floor(f/(2^(j+7)));
    f=f-mius(j)*(2^(j+7));
end     
miu=uint16(mius);
for k=1:7
    out_temp=bitxor(d,miu(k))