% Zobrazeni P(theta)

ThetaT=0:0.01:pi;
Pg1=Pmax1*sin(ThetaT);
Pg2=Pmax2*sin(ThetaT);
Pg3=Pmax3*sin(ThetaT);

plot(ThetaT, Pg1, ThetaT, Pg2, ThetaT, Pg3);
xlabel('Theta (rad)'); 
ylabel('P (MW)');
legend('Bezporuchovy stav','Poruchovy stav','Stav po vypnuti poruchy');
title('Zavislost dodavaneho cinneho vykonu na zateznem uhlu'); 
grid on;
