find_etaoids

Given the normalizing matrix A, we compute the etaoids
(in terminology of E. Forest)
 these can be related to the dispersion

0001 function [H1 H2 H3]=find_etaoids(A)
0002 %Given the normalizing matrix A, we compute the etaoids
0003 %(in terminology of E. Forest)
0004 % these can be related to the dispersion
0005 
0006 I1 = ...
0007 [1 0 0 0 0 0;...
0008  0 1 0 0 0 0;...
0009  0 0 0 0 0 0;...
0010  0 0 0 0 0 0;...
0011  0 0 0 0 0 0;...
0012  0 0 0 0 0 0];
0013 
0014 I2 = ...
0015 [0 0 0 0 0 0;...
0016  0 0 0 0 0 0;...
0017  0 0 1 0 0 0;...
0018  0 0 0 1 0 0;...
0019  0 0 0 0 0 0;...
0020  0 0 0 0 0 0];
0021 
0022 I3 = ...
0023 [0 0 0 0 0 0;...
0024  0 0 0 0 0 0;...
0025  0 0 0 0 0 0;...
0026  0 0 0 0 0 0;...
0027  0 0 0 0 1 0;...
0028  0 0 0 0 0 1];
0029 
0030 
0031 H1 = A*I1*inv(A);
0032 H2 = A*I2*inv(A);
0033 H3 = A*I3*inv(A);