find_betaoids

[H1 H2 H3]=find_betaoids(A)
Given the normalizing matrix A, we compute the betaoids
(in Forest's terminology)
 these can be related to the invariants

0001 function [B1 B2 B3]=find_betaoids(A)
0002 %[H1 H2 H3]=find_betaoids(A)
0003 %Given the normalizing matrix A, we compute the betaoids
0004 %(in Forest's terminology)
0005 % these can be related to the invariants
0006 
0007 I1 = ...
0008 [1 0 0 0 0 0;...
0009  0 1 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  0 0 0 0 0 0];
0014 
0015 I2 = ...
0016 [0 0 0 0 0 0;...
0017  0 0 0 0 0 0;...
0018  0 0 1 0 0 0;...
0019  0 0 0 1 0 0;...
0020  0 0 0 0 0 0;...
0021  0 0 0 0 0 0];
0022 
0023 I3 = ...
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 0 0;...
0028  0 0 0 0 1 0;...
0029  0 0 0 0 0 1];
0030 
0031 J=jmat(3);
0032 
0033 B1 = A*J*I1*inv(A);
0034 B2 = A*J*I2*inv(A);
0035 B3 = A*J*I3*inv(A);