samer@4: function z=qforma(A,u), z=sum(u.*(A*u));
samer@4: % qforma - Quadratic form (array version) qforma(A,x)=diag(x'*A*x)
samer@4: %
samer@4: % qforma :: [[N,N]], [[N,T]] -> [[1,T]].