A = [2,1;1,3]
L1 = 3.618034;
X1 = [0.525731; 0.850651];
L2 = 1.381966;
X2 = [0.850651; -0.525731];
scf(0);
clf()
t = 0:%pi/50:2*%pi;
plot(cos(t),sin(t),'b');
a =gca(); a.isoview="on";
// changed from dividing by L to dividing by sqrt(abs(L))
xarrows([0; X1(1,1)/sqrt(abs(L1))],  [0; X1(2,1)/sqrt(abs(L1))], 1, 5);
xarrows([0; X2(1,1)/sqrt(abs(L2))],  [0; X2(2,1)/sqrt(abs(L2))], 1, 3);
xstring(X1(1,1)/L1, X1(2,1)/L1,"Largest L-value, E-vector/sqrt(L-value)");
xstring(X2(1,1)/L1, X2(2,1)/L1,"Smallest L-value, E-vector/sqrt(L-value)");
// dividing by L1 again was a mistake, but I liked the way it looked.
xtitle("Pictorial Representation of E-vectors and L-values");
xs2eps(0,'arrows.eps');