function outHex = ByteMult(n1Hex, n2Hex)
% outHex = ByteMult(n1Hex, n2Hex)
% inputs:  two length 2 hex strings (bytes)/ upper or lower case format
% output:  an (upper case) length 2 hex string that is the product of the
% two inputs viewed in the field GF(256).

%Convert inputted hex symbols to binary vectors representing polynomials
n1BinVec = int2baseb(hex2dec(n1Hex),2,8);
n2BinVec = int2baseb(hex2dec(n2Hex),2,8);
% multiply these polynomials in Z_2[X]
n1n2BinVecPre = ZpPolyMult(n1BinVec,n2BinVec,2);

%Apply division algorithm to reduce the polynomial mod X^4+X^3+X+1 to remainder
[qx rx] = ZpDivAlg(n1n2BinVecPre,[1 0 0 0 1 1 0 1 1],2);

outBinVec = rx;
%convert to hex by first converting to integer
outInt = baseb2int(outBinVec,2);
outHex = dec2hex(outInt);