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