
eta1=1;
phi1=0;
eta2=1;
phi2=0;

s1=eta1*exp(i*phi1);
s2=eta2*exp(i*phi2);
a=-s1-s2;
b=-s1+s2;
c=s1-s2;
d=s1+s2;

figure
%plot(real(a),imag(a),'bo', real(b),imag(b),'ro',real(c),imag(c),'ro',real(d),imag(d),'bo')
plot(real(a),imag(a),'bx', real(b),imag(b),'rx',real(d),imag(d),'bx')
axis([-3 3 -2 2])
line([0,0], ylim, 'Color', 'k', 'LineWidth', 0.001); % Draw line for Y axis.
line(xlim, [0,0], 'Color', 'k', 'LineWidth', 0.001); % Draw line for X axis.
line([4/3,4/3], ylim, 'Color', 'r', 'LineWidth', 0.001); 
line([-4/3,-4/3], ylim, 'Color', 'r', 'LineWidth', 0.001); % Draw line for decision regions
title('Decision regions')
xlabel('Re');
ylabel('Im');