how to calculate Galois field number system ?

hi guys i need a code to calculate the following Galois field number system for generating LFSR sequence, i want to modify this code x_gf = gf(x,m,prim_poly) to work for random binary number specified by:

x= dec2bin(rand(10,1)*2^32)

and it should calculate the matrix as in this figure below:the last array is the initial value.