modal assurance criterion

[Daniel Tourde]

Hi !

Has any body written an octave script which implements the "modal
assurance criterion" ?

I would like to trace the variation of the eigenvalues of a matrix D(k)
in function of k. The modal assurance criteriondoes it. The idea is to
compute the scalar product of the new eigenvectors with the old
eigenvectors :

[value_1,egv_1]=eig(D(k));
[value_2,egv_2]=eig(D(k+dk));

However, there are a lot of subcases that can enter into account, I
would like to know if someone has already written a "well proven"
procedure.


Thanks in advance

                        Daniel Tourde
-- 
*********************************************************************
Daniel TOURDE                                   E-mail : address@hidden
The Aeronautical Research Institute of Sweden   Tel : +46 8 634 13 44
P.O. Box 11021 S-161 11 BROMMA, Sweden          Fax : +46 8  25 34 81
*********************************************************************