start = .45
ending = .55
steps = 50
E = [linspace(start,ending,steps)];
e = linspace(start,ending,steps);
n1 = (4.74./e);
n2= (9.17./e);
n3 = (9.74./e);
n4 = (6.74./e);
n5 = (5.83./e);
n6 = (5.35./e);

o1 = .17^.5
o2 = .17^.5
o3 = .46^.5
o4 = .17^.5
o5 = .19^.5
o6 = .16^.5

K = ((n1-round(n1))).^2+((n2-round(n2))).^2+((n5-round(n5))).^2+((n3-round(n3))).^2+((n4-round(n4))).^2+((n6-round(n6))).^2;
[C,I]=min(K)
e = start + I*((ending-start)/steps)

plot(E,K)
xlabel('e in 10^-9 e.s.u')
ylabel('see explanation below')