function fas2simplot(Va,Vb,Vc)
close('all')
a=exp(1j*2*pi/3);
V0=(Va+Vb+Vc)/3; if abs(V0)<1e-10, V0=0; end
V1=(Va+a*Vb+a^2*Vc)/3; if abs(V1)<1e-10, V1=0; end
V2=(Va+a^2*Vb+a*Vc)/3; if abs(V2)<1e-10, V2=0; end
vec=[Va,Vb,Vc];
vec0=V0*[1,1,1];
vec1=V1*[1,a^2,a];
vec2=V2*[1,a,a^2];
factor=1.1;

subplot(2,4,1)
quiver([0,0,0],[0,0,0],real(vec),imag(vec),0,'linewidth',2)
hold
m=max(abs([real(vec),imag(vec)]));
m=factor*m;
axis([-m m -m m])
title('Vector')

subplot(2,4,2)
quiver([0],[0],real(V0),imag(V0),0,'linewidth',2,'color',[0,0.75,0.75])
hold
m=max(abs([real(vec0),imag(vec0)]));
m=factor*m;
if m~=0
    axis([-m m -m m])
end
title('Sec 0')

subplot(2,4,5)
quiver([0,0,0],[0,0,0],real(vec1),imag(vec1),0,'linewidth',2,'color','r')
hold
m=max(abs([real(vec1),imag(vec1)]));
m=factor*m;
axis([-m m -m m])
title('Sec 1')

subplot(2,4,6)
quiver([0,0,0],[0,0,0],real(vec2),imag(vec2),0,'linewidth',2,'color',[0,0.5,0])
hold
m=max(abs([real(vec2),imag(vec2)]));
m=factor*m;
if m~=0
    axis([-m m -m m])
end
title('Sec 2')

subplot(1,2,2)
quiver([0,0,0],[0,0,0],real(vec),imag(vec),0,'linewidth',4)
hold
% quiver([0,0,0],[0,0,0],real(vec2),imag(vec2),0,'linewidth',2)
% quiver(real(vec2),imag(vec2),real(vec1),imag(vec1),0,'linewidth',2)
% quiver(real(vec1+vec2),imag(vec1+vec2),real(vec0),imag(vec0),0,'linewidth',2)
quiver([0,0,0],[0,0,0],real(vec1),imag(vec1),0,'linewidth',2)
quiver(real(vec1),imag(vec1),real(vec2),imag(vec2),0,'linewidth',2)
quiver(real(vec1+vec2),imag(vec1+vec2),real(vec0),imag(vec0),0,'linewidth',2)

m=max(abs([real(vec),imag(vec),real(vec2),imag(vec2),real(vec1+vec2),imag(vec1+vec2)]));
m=factor*m;
axis([-m m -m m])
grid

%% para animar la sinusoidal
% clear all
% figure(1)
% for i=1:1000
%     T(i)=cos(i);
%     plot(T)
%     drawnow
%     pause(0.01)
% end