alpha = [0:1:180];
N = length(alpha);
Vm = 4*220*sqrt(2);
for I = 1:N
    Ea(I) = (Vm/pi)*(1 + cos((pi/180)*alpha(I)));
end
plot(alpha,Ea), grid on, xlabel('alpha [derajat]'), ylabel('Ea[Volt]')