A=[10 20 30; 40 50 60] As = reshape(A,prod(size(A)),1) u = [1 2]' v = [2 3 4]' K = kron(v,u) disp("res1 = u'*A*v ") res1 = u'*A*v disp("res2 = kron(v,u)'*As ") res2 = kron(v,u)'*As