clear
clf
alpha=0:1:180;
N=length(alpha);
Vm=220*sqrt(2);
for i=1:N
    alph=alpha(i);
    alphrad=(alph/180)*pi;
    sinalpi=(sin(alphrad))/pi;
    Vrms(i)=(Vm/2)*sqrt(1-(alph/180)+sinalpi);
    Veff=Vrms(i);
    Daya(i)=(Veff^2)/100;
    Dayamax=Daya(1);
    DayaPersen(i)=(Daya(i)/Dayamax)*100;
end
plot(alpha,Vrms,alpha,Daya,alpha,DayaPersen,'r'), grid on
Hasil=[alpha' Vrms' DayaPersen']