disp(['',kron(a=kron(a=mod(1./hilb(3),2),a),a)*56+32])