call11(B)

function out = call11(B)

A = 4.7;      %%%%% Hyperfine constant in GHz (As in AI.J)
mu_bohrB = B; %%%%% First-order Zeeman Energy/Planck Constant at 1000G in GHz = 1.39 GHz.
g_J = 2;

H = zeros(8);
H(1,1) = 0.75*A + mu_bohrB;
H(2,2) = 0.25*A + mu_bohrB;
H(3,3) = -0.25*A + mu_bohrB;
H(4,4) = -0.75*A + mu_bohrB;
H(5,5) = -0.75*A - mu_bohrB;
H(6,6) = -0.25*A - mu_bohrB;
H(7,7) = 0.25*A - mu_bohrB;
H(8,8) = 0.75*A -mu_bohrB;
H(2,5) = ((3^0.5)/2)*A;
H(3,6) = A;
H(4,7) = ((3^0.5)/2)*A;
H(5,2) = ((3^0.5)/2)*A;
H(6,3) = A;
H(7,4) = ((3^0.5)/2)*A;

[v,e] = eig(H);
A = diag(e);
v;
[vs,vi] = sort(A);
[x, vr] = sort(vi);
vr;
pos=0;

%for i=1:1:8,
   %if vr(i)==8,
      %pos = i;
   %end
%end
%out = (v(1,pos))^2;   %%%%% Hyperfine state |2,2>


for i=1:1:8,
   if vr(i)==1,
      pos = i;
   end
end
pos;
out = (v(2,pos))^2;   %%%%% Hyperfine state |2,1>