function out=Mixcolumns(state)
ss=zeros(4,4);
for i=l：4
ss(1,i)=xorfun(xorfun(xorfun(m_Mul2(state(1,i)),m_Mul3(state(2,i))),state(3,i)),state(4,i));
ss(2,i)=xorfun(xorfun(xorfun(state(1,i),m_Mul2(state(2，i))),m_Mul3(state(3，i))),state(4,i));
ss(3,i)=xorfun(xorfun(xorfun(state(1,i),state(2,i)),m_Mul2(state(3,i))),m_Mul3(state(4,i);
ss(4,i)=xorfun(xorfun(xorfun(m_Mu13(state(1,i)),state(2,i)),state(3,i)),m_Mul2(state(4,i)));
end
out=ss;