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