function y = rs_symbolmultiple( x1, x2)
%本函數實現GF(2^5)域中的乘法
% x1, x2 是輸入，y是輸出
x1 = convert10to2(x1);
x2 = convert10to2(x2);
dx = zeros(5,9);
for i=1:5
    dx(i,i:i+length(x2)-1)=x2;
end
y = x1*dx;
y = mod(y,2);
dp = zeros( 4 ,9);
for i=1:4
    dp(i,i:i+5)=[1 0 1 0 0 1];
end
for i=4:-1:1
    if y(5+i)==1
        y = y-dp(i,1:9);
        y = mod(y,2);
        y = y.*y;
    end
end
y =convert2to10(y);