samer@4: function z=qform(A,u), z=u'*A*u;
samer@4: % qform - Quadratic form: qform(A,x)=x'*A*x
samer@4: %
samer@4: % qform :: [[N,N]], [[N,T]] -> [[T,T]].