ot=0:0.01:2*pi;

beta=pi/2; % reaktivno breme
% beta=0; % ohmsko breme
P1=sin(ot).*sin(ot+beta);
P2=sin(ot-2*pi/3).*sin(ot-2*pi/3+beta);
P3=sin(ot+2*pi/3).*sin(ot+2*pi/3+beta);
%P3=sin(ot+2*pi/3).*sin(ot+2*pi/3);
P=P1+P2+P3;
plot(ot, P1,':g',ot,P2,':r',ot,P3,':b',ot,P,'-k','LineWidth',2)
xlabel('\omega.t / rad')
ylabel('Moc / W')