samer@4: function z=iqform(A,u), z=u'*(A\u);
samer@4: % iqform - Inverse quadratic form: iqform(A,x)=x'*inv(A)*x
samer@4: %
samer@4: % iqform :: [[N,N]], [[N,T]] -> [[T,T]].