e1 = [1 0]';
e2 = [0 1]';
ax1 = [ 1 1 ]'; ax1 = ax1/norm(ax1);
ax2 = [ 1 -1 ]'; ax2 = 2*ax2/norm(ax2);
theta = 0:(pi/64):(2*pi);
plot( e1(1)*cos(theta)+e2(1)*sin(theta), ...
e1(2)*cos(theta)+e2(2)*sin(theta), ...
'r-')
hold on
plot( [0 e1(1)], [0 e1(2)], 'r--d' )
plot( [0 e2(1)], [0 e2(2)], 'r--x' )
plot( ax1(1)*cos(theta)+ax2(1)*sin(theta), ...
ax1(2)*cos(theta)+ax2(2)*sin(theta), ...
'b-')
plot( [0 ax1(1)], [0 ax1(2)], 'b--d' )
plot( [0 ax2(1)], [0 ax2(2)], 'b--x' )
hold off
grid on
axis('square')