cs=crystalSymmetry('cubic');
ss=crystalSymmetry('triclinic');
d=rotation('Euler', 135*degree, 90*degree, 354.736*degree,'Bunge'); %d is the burgers relationship
phi_medio=54.74;
theta_medio=45;
psi_medio=0;
rot=rotation('Euler',phi_medio*degree, theta_medio*degree, psi_medio*degree, 'Bunge');
ori_beta=orientation(rot,cs,ss); % orientation reconstructed beta
s=symmetrise(ori_beta);
for i=1:24
var(i)=s(2*i-1)*d;
end
cs=crystalSymmetry('hexagonal');
ss=crystalSymmetry('triclinic');
% 12 theoretical variants from reconstructed beta
or(1)= orientation(var(1), cs,ss);
or(2)= orientation(var(3), cs,ss);
or(3)= orientation(var(5), cs,ss);
or(4)= orientation(var(7), cs,ss);
or(5)= orientation(var(8), cs,ss);
or(6)= orientation(var(9), cs,ss);
or(7)= orientation(var(11), cs,ss);
or(8)= orientation(var(12), cs,ss);
or(9)= orientation(var(13), cs,ss);
or(10)= orientation(var(14), cs,ss);
or(11)= orientation(var(18), cs,ss);
or(12)= orientation(var(19), cs,ss);