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); %輸出乘法結果,大于8位
end
f=double(d);
for j=1:7                      
    mius(j)=floor(f/(2^(15-j)));
    f=f-mius(j)*(2^(15-j));
end     
miu=uint16(mius);
e=uint16(f);

out_temp=bitxor(e,miu(1)*19);
out_temp=bitxor(out_temp,miu(2)*135);
out_temp=bitxor(out_temp,miu(3)*205);
out_temp=bitxor(out_temp,miu(4)*232);
out_temp=bitxor(out_temp,miu(5)*116);
out_temp=bitxor(out_temp,miu(6)*58);
out_temp=bitxor(out_temp,miu(7)*29);

out=double(out_temp);